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PREFACE 


The study of electrical transients is a branch of the study of 
vibrations, and has always been of great mathematical interest. 
With recent developments in engineering it has now become 
of great practical importance. Certain branches of engineering, 
as for example communication engineering, are essentially the 
study of controlled transients. 

This book deals especially with electrical engineering pro- 
blems. The formal mathematical treatment has been discarded, 
and the Heaviside operational methods, now widely used in 
technical literature, adopted. I do not profess to have given 
a rigorous treatment of operational methods, but reference to 
tl^e writings of Bromwich, Jeffreys, Carson, and others will 
f^adily supply this deficiency. 

j!iS6tually, the work may be regarded as consisting of two 
inain sections. The first section. Chapters I-VI, gives the theory 
of lumped circuits, wherein the expansion yields negative 
powers of the operator p. The second section deals with smooth 
circuits, or repeated lumped circuits, wherein expansion gives 
fractional positive powers of p. Problems of the latter section 
are analogous to those of weighted strings. On writing down 
a differential equation for such a string the tacit assumption 
is made that dyidx is continuous, but this need not be the case. 
To avoid this difficulty the modern method is to view the system 
as the limiting conditions of a discrete system. This method 
has been adopted. 

In Chapter X I have discussed Fourier series and Fourier 
integrals in order to give further insight into operational 
methods and to show that the Fourier integral yields a long 
hst of operators and their equivalencies. * . ,! 

Some of the iaethods*6f dealing with variable circuits have 
been discussed in Chapters XI and XII. This is a very wide 
subject, and considerably more mathematical research is needed 
to put it on a firmer (i.e. more rigorous) basis. 

Throughout the work the technical press has been freely 
drawn upon for practical illustrative matter. In several investi- 
gations, simplifying assumptions have been made in order to 
bring the work within the compass of this book. Indebtedness 
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to the writings of others is evident from the many references 
in the text. 

I am also indebted to my colleague, Assistant Professor 
W. O. Richmond, M.A.Sc., for reading the manuscript and 
for several helpful criticisms, and also to the publishers for 
their painstaking work. 

W. B. C. 

Vancouver 
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TRANSIENTS IN 
ELECTRIC CIRCUITS 

CHAPTER I 

GENERAL THEORY 

1. Definitions. Transient conditions refer to what happens 
between steady-state conditions when the initial stability is 
disturbed. They are usually of short duration and of decreasing 
amplitude with respect to time, space, or both. In engineering 
practice, studies of these conditions are of increasing import- 
ance; for example, communication engineering is in reality 
the study of controlled transients. 

It is usual to segregate transients as free or forced. In the free 
condition, the transient is occasioned by switching or some 
such change of circuit conditions. The duration of the disturb- 
ance will depend on circuit constants. In the forced condition 
the transient is set up by induction, radiation, etc., and the 
period of the disturbance will depend on the source creating 
the disturbance. In the solution of differential equations the 
particular integral yields the forced vibration, while the 
complementary integral gives the free vibration. 

2. Representation. The mathematical representation of a 
transient term will now be explained. Every engineer is 
familiar with the fact that a harmonically varying quantity 
may be represented by A sin (ot, and that such a quantity may 
be represented as a vector of length A rotating in a plane about 
the origin with angular velocity 27r/, where /is the frequency. 

For transient conditions the quantity varies with respect to 
time or space exponentially. Then the length of the vector 
instead of being constant (= A) would vary as Ae~°‘‘ or Ae~ “®, 
and the full expression would be Ae~'^ sin cot or sin cox. 

The locus of the extremity of the vector would now be a 
centre-seeking spiral. 

3. Derivation of Network Equations. The basic laws of the 
differential equations in electrical work, from which the unknown 
may be found in terms of the known quantities, are Kirchhoflf’s 

1 
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laws. These are usually considered for steady-state conditions, 
but apply equally well to transient conditions. These 
laws are — 

1. The sum of the e.m.f.s in any mesh or closed circuit is 
equal to the sum of the potential drops in that circuit. 

2. The sum of the currents at any junction in a network is 
zero. Current flowing towards the junction is reckoned as 
positive, and that leaving is considered as negative. 

There are also certain theorems which are of use in setting up 
the equations in certain instances, as follows — 

(i) Superposition Theorem. The current that flows in a 
neWork, due to the simultaneous action of a number of e.m.f.s 



Fig. 1 


appUed in any manner, is the sum of the currents due to each 
of the component e.m.f.s acting separately. 

(ii) Reciprocity Theorem. If an e.m.f., located at a point 
in a circuit, produces a certain current at any other point in 
the circuit, then the same e.m.f. acting at the second point would 
produce the same current at the first point. In other words, 
the receiving point and the sending point may be interchanged. 

(iii) Thevenin’s Theorem.* The current in any circuit is 
equal to the open-circuit voltage divided by the sum of the 
external and internal impedances. For the circuit shown in 
Fig. 1, consisting of impedances Z and Z^, where Z includes 
the impedance of the source, let the voltage across ab with Z^ 
absent be e^. Then the current with present is given by 

* — ^t/{^ + ^i) 

(iv) Compensation Theorem. If a network is modified by 
making a change AZ in the impedance of the network, the 
current increment A/, produced at any point in the network, 
is equal to the current that would be produced at that i)oint 
by a compensating e.m.f. of value / . AZ acting in series with 
the modified branch, where I is the original current in the 
network. 


* Comptes rendtis, 1883, 97, 159. 
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4. Differential Equations of a Network. In all our subsequent 
discussions we assume that the current in every branch is 
directly proportional to the voltage, i.e. the parameters or 
coefficients of the equation are constants. In a circuit con- 
sisting of a resistance, inductance, and capacity in series we have 
that the voltage drop in the resistance is Ri by Ohm’s law, 
in the inductance is L . dijdt by Lenz’s law, and in the capacity 

is J ijC , dt. So that, if V is the applied voltage, we have 
V = Ri L . dijdt -f ji/C . dt 

= {R L . djdt -f Jl/C* . dt)i 

— Zi^i 

where denotes the self-impedance of the circuit. When 
circuits are coupled as shown in Fig. 2, there is a mutual effect 


Za Zx 



Fig. 2 


between the currents flowing in the two circuits. For the 
current directions as shown, we have the following equations 

— A^C + ^2 ~ ^ 

whence 

A(^a + Z,) (Z, + Z,) - I,Zf = (Z, + Z,)B, - Z,E, 

i T ^ {Zb + ^c)^l — ■^c-^2 

^ {Z„-i-Z,)(Zb-\-Z,)-Zf 

Denoting by Z^^ the self-impedance of circuit 1, viz. 
{Za -j- Zg), by Z 22 the self-impedance of circuit 2, viz. 
{Zc -i- Zb), and by (= ^ 21 ) mutual impedance of 
circuit 1 with respect to circuit 2, viz. Z^, we have, on 
substitution, 

Z^Bi Z12B2 



4 


TBANSIENTS IN ELECTRIC CIRCUITS 


For a complex network we get 

4 “ ^ 12-^2 “ 1 ” ^ 13-^3 •••■!“ ^in^n 
^2 = ■^ 21^1 + ^ 22^2 ^ 23-^3 • • • + ^2n^n 

and so on. 

From this array of n linear equations we obtain the equation 
for I„ in a concise form by use of determinants as follows. 

Let D be the determinant of the impedance operators, viz. 

■^11 ‘^12 • • • ^in 

^21 ^22 • • • ^ 2 n 

^nl ^n2 • • • ^nn 

Then the minor Mat, is formed from the above by deleting 
the ath column and the 6th row. 

We then obtain the following set of differential equations — 

— El . MiijD + E.^ . MiJD En ‘ M^jD 

I 2 = El . M 21 ID 4“ E 2 . M 22 ID 4- • . . 4* • ^2nl^ 

and so on. 

The general differential equation is given by 



« = 1 


where s = 1, 2, . . . n. 

5. Heaviside Notation. Following Heaviside, we introduce the 
symbol p = d/dt ; the inverse operation of integration we denote 

by 1/p = Prior to Heaviside, this notation had been 

used by both Caque and Boole (1864-5). For repeated differen- 
tiation we have p" = d”fdt^, and for repeated integration we 

have 1/p" = . . . w times. This latter statement may be 

demonstrated by using Caque’s method of solution. Thus, 
suppose we have 

Mi -f- Lpi — E 

Then pi = (E — Ri)/L 

i “ {^/E |* Edt — {MjL) r idt 
Jo Jo 


or 



GBNEBAL THEORY 


5 


On integrating the first term and substituting the whole 
expression for i we obtain 


m/L - Edt - (i?/X) J idt^ 

E ^ ER C C . 


Repeating this process, we get 


._E ER ^ EE^ 

L* 2.2 • 2 ! 2.3 ■ 3 ! 

where a = RjL. 

Another method of dealing with this equation is what 
Heaviside calls algebrizing. On writing the original equation in 
symbolic form we get 

. E . y, E .J. 

^ "5 i rZ. ’ '■^1 r „n i /D#r \ ' ' 


R + Lp 


Lp[l + (R/L) (1/p)] 


Eri_R j_ , J_ 
2.Lp 2. >2 + 2^2 -^3 


. (3) 


Comparing like terms in equations (2) and (3), we obtain that 
(1/p) is equivalent to* t 

(l/p2) = f2/2! 

and iyp”) F= t"/n\ . . . . • (4) 

Incidentally, we also obtain that 

(2) = (l/«)(l-£-“0 . . . (6) 

Also, it may be demonstrated that the operator p obeys the 
usual laws of differentiation. Thus if A is a constant, and 
X, y variables with respect to time, we have 

p{Ax) = Ap{x) 
p{x + p) = p{x) + p(p) 

p”p'"(a:) = p(« + »») (a;) 


* The sign == means “is equivalent to” wherever it occurs in the text. 
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A word of caution must be inserted for expressions like 
(i), which is not equivalent to pe"®' (1). Thep operates 
only on subsequent terms, and not on prior terms. 

In the above discussion we have used a symbol introduced 
by Heaviside, namely (i). This denotes his unit function. He 
defined this function as having a !',ero value before t — 0 and a 



constant value of unity thereafter (Fig. 3). So it represents 
a square wave of infinite period. For the beginner it is best 
to consider it as a battery of unit e.m.f. switched on to the 
line at f = 0. The mathematical basis of the idea will be dealt 
with later. 

6. Solution ol the Differential Equation. From Section 4 we 
have the differential equation 



«= 1 


Dropping the subscripts, and, in conformity with usual 
practice, writing Y{p) and Z{p) in place of M and D, we have 

Y(P) 


E 


Z{p) 


Now Z{p) will be a polynomial in p and may be written as 
Z{p) = p” + + bp'*~^ -h . . . 

and Y (p) wiU be a polynomial in p also, but of degree less than 
that oi Z{p). In view of the simplicity and straightforward 
nature of Heaviside’s method of solution of such equations, we 
will proceed with his methods and neglect the formal mathe- 
matical methods which are obtainable elsewhere. Heaviside 
now proceeds as follows. 

(i) By algebrization of the Y{p)IZ{p) term in a power series. 
Frequently the result is difficult to convert to a known function. 
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This method has already been illustrated in Caqu6’s method of 
Section 6. 

(ii) By the expansion method, described in his Electromagnetic 
Theory, Vol. II, p. 127, where he gives the following expansion 

“The Z in the operational solution is an operator, a function 
of p, the time differentiator. But Z is entirely algebraical. 
Thus Z(0) is the algebraical function obtained by putting p = 0 
in Z. Then, in the summation, dZfdp is the ordinary differential 
coefficient of Z with respect to p as a quantity. Lastly, the 
summation ranges over all the roots of the algebraical equation 
Z = 0.” 

Heaviside merely sketched the proof of this expansion. 
However, several rigorous proofs now exist (see Bibliography), 
and we now give a method due to Goto* which is based on a 
method duo to Press.! 

The fundamental equation is given by 

i = E[Y(p)IZ{py\ 

By a well-known theorem on partial fractions! we have that 
f[x)=f{x)l{x~-a)-^f{x)l{x-b)+ . . . 
where f{x) — {x — a) (x — b) . . . 

So with X = awe get 

[/'(«)L = o= (a-c) . . [/(a; )/(.'« -«)]*_„ 

Also 

As Y{p) and Z{p) are polynomials in p, we may write 

Y{p)IZ{p) = A/{p-pi) + B/lp-pz) + ... . (6) 

where p^, p^, p^, . . . are the roots oi Z — 0, and A, B, C . . . 
are constants. To determine these constants we proceed as 
follows. 

Cross-multiplying equation (6), we obtain 

Y(p) = A{p-p^) ip-Ps) . . . +B(p-pi) 

ip-Ps) ■ ^ ‘ +Gip-Pi) {P-P%) ■ ■ ■ 

* Researches of Electrotechnical Laboratory, Tokio, p. 201. 
t Trans. A.I.E.E. (1917), Vol. 36, p. 232. 

X Williamson, Integral Calculus, p. 49. 
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With p — Pi, we get 

Y {Pi) = A{p^- p^) {p^- p^) . . . 


(Vi-Pi) iPi-Pi) • • • 

By the preceding theorem 

A==[YIZ’y^^, 

Similarly B = [F/Z'],_p. 

So that 

riz=ij{p-p,).iriz'],.,, + iHp-p,).iYiz'i.,.+ . . . 

We consider the case when the apphed voltage is given by 
jE7e“‘. Then 

i=^E{lnp-p^).[Y|Z'],.,^ 

+ i/{p-Pz)4YIZ'],,„+ . . . }e“‘ . (7) 

Now, by power expansion 

y(p — ?).£“*= (1/p) (1 + qjp + +...)«“* 


Noting that 

(1/p) . e“‘ = J" e’^dt = (1/a) (e®* — 1) 

(l/p2) .e“‘ = J‘ (1/a) (e“*- \)dt =(l/a*) (e“'- 1)- (</a) 

and generally 

(1/p”) . £“* = (l/a")£“‘— 1/a" - </a"-i 

— < 72 !a ”-2 . . . — f"-V(w— l)!a 

So that 

1/(P — ?) • e"* = (Vp) [e“‘(l + ?/« + • • •) 

— (?/«) (1 + ?/« + • • •)- (ff^Va) 

(1 + g/a + gVa® ...)•••] 

= (1/p) [(e“‘ - g/a - g*</a . . .) • 1/(1 - ?/«)] 
= (1/p) (£“‘-g£«ya).l/(l-?/a) 


Carrying out the integration, we get 
1/(P ~ ?) • = («“* ~ £**)/(« “ ?) 
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Substituting this expression in equation (7), we get 
i = I!{[l/(x- p,)] (e-‘ - e^‘‘) [ r/Z% . + [l/(« - p,)] 

= £Je^{[l/(a-p,)] [r/Z%^„ + [l/(«-Pa)] 

[r/Z%^^,+ . . .}-^{[e^‘V(a-Pi)][F/Z']^„^. 
-j-[aP‘‘/(o^-p,)][Y/Z’]^.„-i- . . .} 

Now 

[i/(oc-Pi)][rm,^^, + [i/(oc-p,)][r/z%„„ 

+ ... = [Y/Zl^. 

So that we get 

i = ^ 2” . (8) 


The value of a is determined by the form of the applied 
voltage. For a d.c. voltage, a = 0 , and then equation ( 8 ) 
becomes 


i = ^ {[ Y/Z], „ 0 + [ Y/p„Z% „ } 

which we will call the^rsi kind of expansion. 
For an a.c. voltage, a = jco, and we have 


which we call the second kind of expansion. 


. ( 9 ) 


. ( 10 ) 


For a voltage given by E cos cot, we obtain, on rationalizing 
and discarding imaginaries, 


i = E 


fJL 


+ T 


P« 


P = JO) 


^Pn + 




7. Routine to be Adopted. For a beginner with this method 
of solving equations it would be well to adhere to a fixed 
routine such as the following. 

(i) Deduce the equation using the “p” notation. Suppose 
we get pH + bi = E . {!). This, according to our inter- 

pretation, means a constant voltage E is applied at < = 0 , and 
is maintained indefinitely at that value. 

2— (T.152) 
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(ii) By algebraic manipulation, we have 


E 

ap b 


Al) 


Here note, Y {p) = 1 and Z{p) — p^ ap h. On equating 
Z{p) to zero we get p == [— a i Vi^^~ = PiPi- 

(iii) By differentiation, we have Z'{p) — 2p -{■ a. 

(iv) Forp = 0, we get Z(0) = b and F(0) = 1. 

(v) By collecting terms and substituting in the expansion 
expression 

i = E{[Y/Z],^o + Y[{Y/pZ')er^y^,J 

we get 

i = E[l/b + e*’‘‘M(2pi + a) + + «)] 


(vi) On substituting the values for p^ and p^, the complete 
solution of the differential equation is obtained. It will be 
noted that by this method no integration constants are left to 
be evaluated. This is because of the assumption that at i — 0 
the circuit is dead. This, as will be seen later, does not hold 
for all cases, but by a simple extension we are able to bring 
our solution into line with any circuit conditions. In Chapter II 
we will follow this routine in solving certain circuit problems. 

8. Discussion on Boots of Z{p) — 0. The solution of the 
equation by means of the expansion theorem depends on our 
ability to extract the roots oi Z{p) = Q. A little consideration 
will show the difficulties we are hable to encounter in analytical 
work, but, fortunately, in a numerical example, we are able to 
approximate to the roots by either Horner’s or Graeffe’s method. 
These are explained fully in many texts.* 

If the roots obtained from Z{p) = 0 are all different, then 
the application of the expansion theorem is straightforward. 
But when certain of the roots are equal, then Wagnerf shows 
that for roots p^, Pz ... p„ and p„, the latter occurring m 
times, we have 


i = E 


r(0) 

m 


+ 



Y(Pr) 

PrZ'iPr) 



-}- -f- 


♦ E.g. Whittaker and Robinson: Calculus of Observations (Blackie). 
t Elect, f Ur Elek,, 1916, 159. 
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where . . . are constants determined by terminal 

conditions. 

Also, it should be noted that the proof of the expansion 
theorem is dependent on Y{p) being of degree less than that of 
Z{p). If this condition is not fulfilled, we must divide through 
until our remainder R{p) is of power less than that oi Z{p) as 

Y{p)IZ{p) = B^p‘^ B ^ . . . 4-i2(p)/Z(p) 


Frequently it is possible to get out of this difficulty by 
solving for ijp, say, and then differentiating the result to 
obtain i. 

Prodttcts. Sometimes the impedance operator Z{p)v& given 
as the product of two components, Zi{p) . Z^{p). To evaluate 
such an expression, recourse is had to the superposition theorem 
of Borel (see p. 66), which may be written operationally 
as follows — 

IfZip) = A{t) = djdt J‘ A^- A) . A^{X)dX 

= AS) . ^2(0) + £ AS - A) . ASm^ 

where 4i(<) = llZj^{p) and Aa(<) = \IZ^{p). 

We will illustrate this important theorem by an example. 
We have already shown in equation (5) that 


1 

p + a 


. (i) = l[l-a-“‘] 

a 


Then by differentiation — — — . (1) = £““* 

Suppose now that we are given 

1 a p 

Z(p) ■ p + a " p + a 

Then we have 


and 


1 a 

Ziip) ^ p + a 
1 ^ P 

^2(P) P + a 


(i) = (l-e-“‘) = A(0 


(i) = = AS) 


By substitution we get 

A{t) = [1 - £-“*] + £ [1 - £-«(»- -»)] [- 0Le-^]dX 
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On integrating, we get 

A{t) = [1 — e-“‘] + ede-^ + 


Highly Oscillatory Circuit. Guillemin* shows that for 
a low dissipation circuit we may split our determinantal equation 
into two parts, and, by combining the roots so obtained, get a 
close approximation to the actual roots. For a full discussion 
of this method we refer the reader to the original, but the 
working may be illustrated by the following example taken 
from the original — 

Let 50p* + 8*91 x + 2-03 

X Wp + 10“ 

Now let D{p) — p* + 8'91 x 10®p^ + 10“ = 0 

This gives p — i jSOO, orp = ii^OO. 

Also, let 0{p) = 50p® + 2-03 x lO’p 

On substituting the above values of p, we have 

0{p) = ±i9-36 X 10«, or Q{p) = ± j3-90 x 10» 

Then 


djdp . [D(p)] 


Let 


Then 


or 


d 


d = 


4p3 + 17-82 X 105^ 

G{p) 

D'{p) 

j9-36 X 10» _ 
j6-23 X 108 

j3-90 X 10® _ 


10 


j3-90 X 108 

So that the approximate roots (^ + p) are 
- 10 ± j600. 


15 ± j800 and 


9. Notes on Roots. Frequently in discussions of circuit 
conditions the prime consideration is stability, and so we need 
not solve our equations, for such information (i.e. regarding 
stability) may be obtained from the determinantal equation. 
The underlying theory can be obtained from the theory of 
equations, and we give below a brief summary of some of the 
rules involved. 

* Trans. I.E.E. (1928), Vol. 47, p. 361. 
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General. 1. For /(a;) = 0, there cannot be more positive 
roots than there are changes in sign of f{x). Also, there cannot 
be more negative roots than there are changes in sign of/(— x). 

Thus f(x) = a;® + 5x ^ — x^ lx 2 — 0 

Here there are two changes, so we may expect two positive 
roots. 

Also /(— a:) — — a:® + 5^* + a;® — 7a: + 2 = 0 

Now, there are three changes in sign, and so we may expect 
three negative roots. Altogether there will be nine roots, and 
so the other four must be imaginary. 

2. If all coefficients are real, then the imaginary roots must 
occur in conjugate pairs. 

Quadratic Equation.?. For the equation ax^ bx c = d, 
the roots are given by the well-known expression 

a: = [ — b ^ ^/ - 4ac)]/2a 

If all coefficients are positive, then the roots are negative. 
For roots to be positive, the coefficient b must be negative 
by general rule 1. 

1. For real roots, 6® > 4ac. 

2. For equal roots, 6® — 4ac. 

3. For imaginary roots, 6® < 4ac. 

As these will be of frequent occurrence in the text, we refer 
the reader to later discussion. 

Cubic Equations. Here the general equation is of the form — 
ax^ + 6a;® -f- ca: + d = 0 

On substituting ( — x) for x, we find, in accordance with 
Rule 1, that there are three changes in sign, and so there may 
be three negative roots. Usually in solving these equations 
[x— (6/3a)] is substituted for x, and an equation of the form 
a:® -f- 5'a: + r = 0 results. 

From this equation we are able to state — 

(а) If r®/4 + (7^/27 is positive, then the roots are of the form 

« + / 5 ; + 

(б) If r®/4 -j- ^3/27 is zero, the roots are 2a ; — a ; — a. 

(c) If r®/4 -j- g®/27 is negative, then a® and /3® are imaginary 
quantities. So let a == y + = 7 

Then the roots are 2y; — y— 6^/3; — y + i.e. real 
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quantities. The conditions for non-oscillation may be shown 
to be 6^ — 3 ac > 0 and — 4 ac® > — 18 abed -f 

Biquadratic Equations. Let us consider an equation of 
the form 

m* + Ciin^ -f CgW® C^m 4- <^4 = 0 

The roots of this equation, if C^, C^, C^, C4 are real and 
positive, will be real or complex conjugate. They will occur in 
pairs, two reals and two complex. Now, it may be shown that 
if Pi, P2, Ps, Pi, are the roots, then 

(i?i + Pt) (i>i + Pz) {Pi + Pi) {Pz + Pz) {Pz + Pi) 

= - CiC^ - C32 = X (say) 

(a) If the roots are 4 : and Og 4 j^z> obtain 

2(ai -f aj) = — C*! • • • (12) 

and 

X. ~ 4 <*j<z 2 [(% “f" ^2)^ 4 (^1 4 ^2)*] [(®i 4 ^2)* 4 (^1 ^2)*] 

( 13 ) 

From equation (12), and for a stable damped system, a^, 
must be negative, and so we have X positive. 

(b) If the roots are 4 j^i and ag 4 ^ 2 > we have 

2 (ai 4 « 2 ) = -C . . . ( 14 ) 

and X = 4aia2 {[(«i 4 «2)’^ 4 ^1® — 4 } . ( 16 ) 

Here we have the same criterion, viz. if a^, a2 are both 
negative, then X is positive. 

(c) If the roots are Pi 4 i?i> :izjsi, then we have the 
following relationships between roots and coefficients — 

Pi 4 = - GJ2 

Pi** 4 4 rj^ 4 «!** 4 ^PiTi = 

r\{Pi 4 gi**) 4 Pi(ri® 4 V) = - <^ 3/2 
(Pi** 4 gi**) (V 4 V) = Gi 

If we make p^ zero, we then get^'g^ and 4 so that 

ri = - C'i/2 

gi** 4 V 4 «!** = Ca 
^igi** = - c'3/2 
gi**(?-i® 4 «!**) = Gi 
Whence we get gi®(^2 “ gi*) = Gi 



GiC^C^-G^^-C^^Ci = 0 


. ( 16 ) 
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This is an important relationship in systems, since for 
stability and must be negative. Now increase to, say, 
Pi = 0. Then equation (16) gives the limiting characteristics 
for stability. 
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CHAPTER II 


APPLICATION OF EXPANSION THEOREM 
TO LUMPED CIRCUITS 


{ A) Making the Circuit with D.C. Applied Volts 


1. We now apply the expansion theorem to simple lumped 
circuits, i.e. with concentrated characteristics. We concentrate 
on the j)roblem of the making of the circuit by the application 
of a d.c. battery, of voltage E, to the circuit. At the time 
i < 0, the circuit is dead, and at time t = 0 a, voltage E is 
applied. We are supposing that an ideal switch to give instan- 
taneous application of voltage, without any leakage current, 
is available. 

2. The RL Circuit. For a circuit consisting of resistance R 
and inductance L, at any instant we have the following equa- 
tion + Lpi — E 


For a suddenly applied voltage, we write, in accord with 
previous discussion — 

Ri -|- Lpi = E{1) 


So that 


i = E 




where Z(p) = R A- Rp and Y{p) = 1. 


Following our routine, we get — 

R -j- Lp = 0 ; p ~ — RJL 
Z'ip) - L 
Z(S>) = R, 7(0) = 1 


So, on substituting in the expansion theorem, we have — 

• =®[5 + LFiP) 

where e is the base of Napierian logarithms. 

As t aj)proaches infinity, the current approaches EfR. This 
is the steady-state current, and is usually denoted by Iq. Our 
equation becomes — - 

i = ^[1 — 


where a = R/L. 
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3. Deductions from Above Equation, (a) On differentiating 
with respect to t, we obtain 

di/dt = 

and at i = 0 we get di/dt — oc/q 

The initial rate of growth is proportional to a, so a coil of large 
resistance and small inductance will reach its steady state 
more quickly than a coil of small resistance and large induct- 
ance. After a lapse of time 1/a = L/R, the current is /o[l — 1/e], 
i.e. 0‘623/o. This value of time is known as the time constant 
of the circuit. 

(6) The voltage drop and power consumption in the com- 
ponents of the circuit are dealt with as follows. Across the 
resistance the voltage drop is Ri, or E[1 — e““‘], and across 
the inductance it is L{di/dt), or Ee~'^K By adding these com- 
ponents we see that tlie sum of the voltage drops is equal to 
the applied voltage at any instant. 

The energy dissipated as heat at any instant is 
Ri'i = {E^IR)[l - 2e-“‘ + e-2«‘] 
and the energy stored in the magnetic field at any instant is 
iLidi/dt) = (^Vi2)£-“‘[1 - e-“‘] 

The power consumed between t — 0 and t = T is given by 

I* RiHt -f f Li{dildt)dt watts 
Jo Jo 

Substituting the value of i and simplifying, we get — 
{EyR)[T+{l/a)e-'^-{lla)] 


When T is large, this expression becomes {EyR)[T — 1/a) 
and as a limiting case it becomes E^T/R. 

4. The Condenser Circuit [C'J?]. In considering the charging 
of a condenser we pass over the interesting case of charging a 
perfect condenser and consider the case of a condenser C with 
a resistance R in the circuit. Initially, let the charge on the 
condenser be zero, so that the voltage becomes 
Ri + {l/pC)i = E(l) 

1 T, Cp 


i.e. 


E 


R + {l/Gp) 


{1) = E 


RCp + 1 


(I) 


Here Y{p) — Cp, and Z{p) — RCp 1, so 7(0) = 0, 
Z(0) = 1, and Z'{p) = RG. 
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From Z{‘p) = 0 we have p = — l/CR. From the expansion 
theorem we get 

The charge ? = f idt = — e““‘], where a = IJCR. 

Jo 

The values of the time constant, the voltage across the circuit 
components, the energy, and the power consumption may be 
determined in the same way as shown above in Section 3. 

5. The RLC Series Circuit. For a resistance, inductance, 
and capacity in series, the equation becomes 
Ri + Lpi + (l/pC)t = £J (1) 


so that i 
Now, 


E 


R + Lp {IfpC) 
Y{p) = Cp and Z{p) 


{1)^E 


Cp 


LCp^ + RCp + 1 
LCp^ + RCp + 1 


U) 


also F(0) = 0 and Z{0) = 1 
The roots obtained from Z{p) — 0 are 



It is usual practice to denote R/2L by a and l/LC by coq*, 

where coq is the natural angular velocity. Also 

is denoted by co^. Then — a.^— cwq^. Thus the roots are 
p = — a ± wi- 
By differentiation we get 

Z'ip) = [2LCp + RC]^^,^, = ± ViR^C^- ^LC) - ± 2LCco, 
On substituting in the expansion formula we get 

i = FO To 4- ^2 

^ L Pi- 2LC(Oi Pa . 2LCcOi 

= (EI2Lcoi) 


J\4tL^ LCj 


Inserting the values of p^ and pa, we get 

i = (F/2ft)iL) e""**] . . (19) 

Further development depends on the value of coi. From 
Section 9, Chapter I, we see that there are three cases — 

(i) If (Oj^ is positive, i.e. RyiL^ > l/LC, we get 
i = {E/(OxL) e“ sinh coit 


(20) 
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The value of i approaches its final value slowly, and the condi- 
tion is known as overdamped. 

(ii) If is negative, i.e. If LG > Ry4:L^, then we must 
substitute jwi for co^ in the above equation, and so obtain 

i = (E I cd^L) e~'^ sin co^t . . (21) 

In this case the current is oscillatory, its period is equal to the 
free period of the circuit, and the circuit is said to be under - 
damped. 

(iii) If (Oy = 0, i.e. — IfLC, then sinh coit approaches 

(Ojt, and thus the limit of (sinh (o^t)l(x>i is t. 

This gives i={EIL)te~^ . . . (22) 

The current here approaches its final value without oscillation 
and at a faster rate than under (i). This condition is known 
as critical damping, and the value of R obtained from 

R^jil? = IjLC 

is 2 (L/C) and is known as the critical resistance of the circuit. 


{B) Making the Circuit with A.C. AppUed Volts 


6. The RL Circuit. Since an alternating voltage may be 
represented as the real part of Ee^’^* = £/[cos cot -f- j sin cot], 
we now apply equation (10) of Section 6, Chapter I, namely — 

Lzo'o)) ,4iO«i-Pr)Z'(3>,),.,.J 

The values of Y{p), Z{p), Z'{p), and the roots are identical 
with those obtained in Section 2, Chapter II ; so by substitu- 
tion in the above formula we get 

R-fj^^‘^~ [jco+{R/LyiL I 
’cos cot j sin cot "j 

R + jcoL] 


t = jE( 


E 


"c 


R + jcoL 

Rationalizing and discarding imaginaries from the above 
equation, we have 

~R cos cot -(- coL sin cot R 


E 


[■ 


R2 + 0,2^2 

Letting tan ^ = coL/R, we get 

i = [E/^/{R^ -4- 0 )*^*)] [cos {cot-<f>) 


ij2 4. ^2^2 


g-{RIL)t 


] 


cos^.e-W^)‘] . (23) 
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The steady-state current will be {EfZ) cos {cot — <!>). It will be 
noted that the time constant of the circuit is not affected by 
the type of voltage applied to the circuit. 

We next consider a voltage applied to the circuit, 

i.e. the voltage is Es^^ when the time is zero. Then we have 


" + 6 ) 

_R -f- jcoL 


R + j(oL 


Discarding imaginaries, we get 
i = [E/V{R^ -f co^L^)] 

[cos {(ot+ d-<f>)- cos {0-cf>). e- (24) 


When 0 = 0 it is seen that this expression reduces to that 
given above. The transient term depends on the value of the 
angle 6. When 6 = cf), we get the maximum transient current. 

We may regard the transient term as raising or lowering the 
time axis. On this distorted axis the steady-state condition is 
superimposed. So, if we produce the resultant steady -state 
current back to meet the initial current, it will be found that 
sometimes they meet without any phase distortion but at 
others do not do so. In the former case there would be no 
transient. 

7. The CR Circuit. The voltage equation may be written as 
Rpq + qlC = + (1) 


Completion of the solution along the lines of the last paragraph 
gives 

q = [EC/Vi^ + co^R^C^)] 

[cos {cot -1- d — <f))— cos {d — (f>) . e~ 


where tan ^ = coCR. By differentiation we get 
i = - [EG/^{1 -I- (om^G^)] 

[co sin («< -f 0 - - (l/GR) cos ( 0 - .^)e - (26) 


8. The RLG Circuit. For an applied voltage Ee^"‘^ we get 
i[R + Lp+ (l/pG)] = Ee^<^^{l) 

* "" R + Lp-\- {l/pG) 

Adopting the notation of Section 5, Chapter II, we have 
p = — a d: and Z'{p) = dr ^LGco^, and on substituting we 

get 


i = EG 


(o^LG jcoRG -j- 1 dr {ja — p„) . 2LGcoi 
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On simplifying and discarding imaginaries, we have that 


EG 


I 


CO cos {cot + (^) 


‘ at 


+ 


+ (1 - co^LC f] ^ CO^LG 
[(a® — (o^Y + sinh co^t — ^occoyco^ cosh co^n 

(a^ — co ^ — co^)^ + J ” 

where to.n^ = {\ — (o^LG)lcoEG . . (26a) 

The development of this equation depends on the relative 
values of and 1/LG. 

(i) When E^j^tL^ > l/LG, then co/^ is positive and we are 
not able to simplify the above expression. 

(ii) When Ey4:L^ < l/LG, then co^ is negative and we have 

cos {cot + 


(27) 

(27a) 


^ + {ooL- l/coGf] 

g-a« ^0^2 ^^2J pQQ I 

co^L \/[(a^ + ctq^— co^)^ + 4a^a)^]j 
where tan cp = [(a^ + + a>^(a^ — co^)]l2(xcoyCo^ 

(iii) When coj, = 0, the circuit is critically damped. We get, 
as coi approaches 0, that (sinh co^t^co^ approaches t and cosh co^t 
approaches 1, so equation (26) becomes 

cos {cot + (f>) 


E 


^[E^ + {coL- 1/coG)^] 


+ 




2oico^ 


+ 0)2 


(a^ + 


]) 


(28) 


It is of interest to see the effects of impressed oscillations. 
Substituting a = i2/2L and coq® — l/LG, equation (27) becomes 


. E 

% = V 


CO cos {cot + </>) 


+ ((yj, 2 _ ^ 2 ) 2 J 

where is a constant. 


.de “ cos {co^ ip) 


It is evident that — 

(а) The system gives up its natural period (a)^) and assumes 
that of the impressed force {co) with lapse of time. 

(б) When the impressed and natural frequencies are almost 
equal and the damping coefficient a is small, the two com- 
ponents will assist each other during one interv/il and neutralize 
during the next. The system would thus produce a series of 

beats.” 

(c) When co = co^ and a/co is very small, then coq* —<«* = «• 
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and tan <f> = <x./2(o, which approaches zero. The amplitude may 
be written as 

, Eo) E 

* ~ L^/[4:(o^a.^ + (coo^— w2)2] ~ LV'(4a>V + a*) 

Eco E(o E 

olLV {ico^ + a^) 2xLa)V (1 + 

which is very large. 

9. Numerical and Logarithmic Decrement. Using the expres- 
sion for current as given in Section 6, Chapter II, viz. 

i = {Ef(OiL) e “ sin cojt, 
and differentiating, we get 

difdt — {Ejoi^L) e ~ [coj cos (o^t — a sin 

On equating dijdt — 0 , we obtain that tan = (cojot.) or 
generally tan {coit 4- 2»7r) = (ft>i/a) where n is an integer. So 
at t„ = [tan~^ (<^ 1 /“) — 27i7r]/cai the current is a maximum. 
By substituting this value of time, we get that 

imax — {ElcDiL)e~°^”sin [tan~^ (coja)— 2n7r] 

The first maximum will be obtained when n = 0 — 
hmax ~ {Elo)iL)e~'^« sin tan~^{o)Jcf.) — [E{L\/{c>? -f coi2)]e“ 
For w = 1 we get the second maximum as 
Hmax = [E|L^/{<x^ + a)i2)]e-«t. 

Or generally we get 

inmax = [E/L\/{a.^ -f (Ui2)]e-»‘n 

The interval between successive maxima is 27r/a)i and the 
difference of the currents is 

in-in^i = + O] [e - - e - 

= [E/L^im^ V)] [1 - >e-^n 

or = 1 _ e + «2(«M) (29) 

This expression is known as the numerical decrement. 

For the logarithmic decrement we take 

log, {i„lin + i) = log. . e“<2"/«<>.) = a(277'/(Oi) 

= R/2Lfi = d (say) 


where /i = (oJ 2 n. 


( 30 ) 
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10. Natural Period and (Jritical Resistances. We have defined 
the natural period of the RLG circuit as 

4:L^) 


II L. 

2tt~ 2^r^/\LC 
So for an undamped circuit we obtain 


and thus we get 


fo’ 

n 


■■ (l/27r) ^/{IILG) 


(31) 


4:7T^ \ / 

= fo- 

.*./o//i= V(H-^747r») 

This latter formula is well known to radio engineers. 

We have seen that the value of the resistance which makes 
the term under the root sign vanish is given by 2\/ {LfC). 
This is the critical resistance of the circuit. 

If the value of the applied frequency in an RLC circuit 
given in equation (26) be altered, then the maximum steady- 
state current is obtained when co^R'^C^ -j- (1 — co^LC)^ is a 
minimum, i.e. when 2ioR^C^ -f 4:(o^L^C^ — icoLC = 0. 


So we get a value of 

R=± Vii^L/CI)- = ± ^/[(l - (o^LC)ILC] (32) 

These may be defined as frequency critical resistances. 

11. General Tuning or Resonance Conditions. The equation 
for forced periodic motion may be written as 

apW + bpO cd = F cos cot 

Then, in accord with Section 8, we obtain that the steady 
state deflection is given by 

6 = F cos {(ot — <!>)/ \/[<o^b^ + (c — aco*)®] 
where cot ^ = (c — aco^)Jcob. 

Let the impressed frequency be continuously varied. Then 
the amplitude of 6 will be a maximum when 

cy® = (c/a) — ^(bfa)^ = Wf^ (say) = resonant angular velocity . 

As the natural angular velocity of the system is given by 
we have Wf^ = Wq® — l{bla)^. 

By substitution, the ampUtude is given by 
2aFlb'\/ (4ac — 6®) = Fjbcon 
where (w„® = (4ac — 6®)/4a® = <o^ — l(bfa)^. 
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Now, the critical value of 6 = 6, 
amplitude of the deflection as 

Let y = CO fo}Q and d 

e = F/cy'[(l 


. r.. in V 62 co 2 I 


1 

b/K 


2\/{ca), so we have the 
b^co^ 


bc^cOf, 

We then have 

y2)a 4^2y2J 


(33) 



When <w ->■ 0, i.e. applied frequency approaches zero, we get 
= F/C, so that 

amplitude 6 1 

amplitude at zero frequency ~ ~ 7^)^ id^^] ’ ' 

Also tan <^ = 257/(1 — y2) . . . (36) 

As these results are very general, we have plotted the results 
in Figs. 4 and 6 for various degrees of damping (5) against 7. 
It is seen that the peak of the oscillations tends toward co/coq == 1 
as the damping is reduced. Inspection of the phase-angle 
curves of Fig. 5 shows that for low frequencies the resulting 
motion and deflecting forces are almost in phase. But as the 
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frequency is increased the motion lags more and more behind 
the deflecting force. At (w/ccq = 1 they are in quadrature. 

12. Plotting of Results. In analytical work mathematicians 
have made much use of non-dimensional quantities. We use 
this idea for plotting the results of equations such as have 
been given in previous sections. By proper selection of a time 
base we are able to use one curve for all problems of a type. 
Thus for the LR circuit we have 

Ri -f- Lpi — E 

If t = {LfR)x, then dt = {LfR)dx. Then, on substituting, 
we get 

dx di 

Ri + L-j-y- = Ri-\-RM = E 
' dt dx 

where A = d/dx or 

= {E/R)[l~c-^] . (36) 

Thus for such a circuit we 
only need to plot the curve 
(1 — e ~ ®) against x, as in Fig. 
6. Such a curve we will desig- 
nate as the master curve. To 
use it on any circuit, all that is 
necessary is to multiply the 
ordinates by E/R and the value 
of X by L/R. 

equations of the second power, 

we have 


* ~ i2(l -f A) 



When this idea is applied to 


or 


p2y + (a -f -i- ccpy = E 

y — — n) 

y (p + a) (p -f /?) ^ ' 


By partial fractions we get 

_ L_l. 

^ j6— aLp + a P-f|Sj 

Thus, considering an RLC circuit as 

p^q + {R/L)pq -f {llLC)q = EjL 
R^ 


{!) 


we now 


, „ R ^ K R^ I \ 

have «. /? - 2 L ^ Lc) 
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In the above equation q has two components, and q^. 
Letting e = {EIL)I{§ — a), we have 


?1 = 

1 

P. 1 1 

1 

p + a 

U /5-«J 

“ ® p + a 




1 

?2 





So that by carrying out the method illustrated for a first 
power equation we get that 

a< = 

= .... (37) 

and = — 

where — y^. So that 

• (38) 

Thus we now plot (1 — e”®) and (1 — e"*') and multiply the 
ordinates by EjLot.{^— a) and EfL^{^— a) respectively. By 
subtracting ordinates at like times, we obtain the values of q. 



Thus in Fig. 7 let OA and OB be identical times. Then the 
value of q = AC — BD. Curves of this nature are useful where 
a considerable amount of checking is necessary. In order to 
plot the master curves, all that is required is a set of tables of 
exponential functions. 
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Practical Examples 

We will now discuss certain applications which may be solv( 
by the matter already developed. Actually we have use 
electrical circuits in our discussion, but we will show that tl 
principles hold for other subjects as well. 

13. Sweep Circuit for Cathode-ray Oscilloscope. For oi 
first example we take that of obtaining a variable time bas 
or sweep, in the cathode-ray oscillograph. In this instrumei 
the current to be examined is applied to one set of coils; 1 
another set of coils, at right angles to the former set, a “sweep 
current is applied. This current causes the beam to be swej 
across the viewing screen, and thus produces a wave forr 
One method is to use a circuit containing a reactance. Fro: 
previous work we know that the current grows according to 

i = iE/R)[l- s-"^] 

= {EjR) [at - aHy2 ! + ...] . . (31 

As the time is here of the order of microseconds, it is seen th) 

i ^ {EIR)at — Et/L . . • (4< 

i.e. the sweep deflecting force varies inversely as L. 

For a particle of charge e moving at a speed v, the equivalei 
coil = ev. If JB is the flux density set up by the sweep magm 
coil, the force on a particle is evB per unit length. But for* 
— mass X acceleration, and I = \at^. 

.'. I = \at^ — \{evBlm)t^ 
and t — \/{2lmlevB) = iLjE 

or i = {LJE)H^{evBI2m) cm. . . • (41 

where I is the width of the screen swept out. 


Electrical 

Mechanical 

Name 

Symbol 

Symbol 

Name 

Time 

t 

t 

Time 

Charge 

q 

6 

Deflection 

Inductance 

L 

M 

Mass 

Current 

qt-^ 
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14. Motion of Indicatii^ Instruments. For allied mechanical 
problems we are able to make use of the principle of duality. 
The table given on page 28 has been drawn up showing corre- 
sponding quantities in electrical and mechanical units. We 
assume that inductance and mass and also charge and deflec- 
tion are corresponding quantities. From a table such as this, 
once the differential equation has been set up, we are able to 
discuss the problem fully by referring to the corresponding 
electrical equation. 

In the case of an instrument movement, a force F acts on 
the moving coil. This force is set up by the interaction of the 
current in the coil and the fleld of the permanent magnet. 
It will be resisted by a force due to inertia, a force due to the 
damping vanes, and by elastic forces set up in the instrument 
springs. For small movements we may take the damping force 
as proportional to the velocity, and the shear or control force 
as proportional to the deflection. 

If A is the deflection, then dAJdt = pA is the velocity and 
p^A the acceleration. Now, the force equation may bo written 

mp^A + ppA tA = F 

where p — resisting force per unit velocity and t = control 
force per unit deflection. 

F 

••• A = — 2-^ ^ • (i) 

mp^ -\- pp T 

for an impulse force. 

Now, Z{p) — mp^ pp T and Z {0) — r, and the roots are 
given by 



Z'(p) = 2mp -+- p = dz 2m/9 


so that A 


= 

_T 


+ 


1 /£»’*« 

2m^ \Pi~ P2 


)] 


Inserting the values of Pi and p 2 , we obtain that 

.-ai 

A = F ■■ 


T 2m|8(a^- 




[a(e/>«-£-^‘)+y3(£^‘ + £-^‘)]| (42) 


As for the corresponding electrical circuit, further development 
depends on the value of /3. 
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(i) If /3 is positive, then the above equation gives the final 
result. The instrument is said to be overdamped, and the 
pointer moves slowly to its final reading without overshooting. 

(ii) If/? is zero, then 

^ = + 

The instrument approaches its final reading more quickly 
than in case (i), and is now said to be critically damped. Usually 
the inherent damping is not sufficient to produce this result, 
and external resistance is incorporated. 

(iii) When is negative, then the root is of the form i 
and we obtain 

A = if’ [- - ^2) cos ipt - <^)J . (44) 

where tan <f> = a//?. 

Here the movement is oscillatory, the pointer swinging about 
its final reading. In this condition the instrument is said to 
be underdamped. 

15. Oscillograph Vibrators. The general theory of oscillo- 
graphs may be obtained elsewhere. Here we will be concerned 
with the interesting case of obtaining a correct trace, i.e. the 
relative amplitudes and phase relationships of the various 
harmonics to the fundamental must be preserved. The dis- 
turbing force varies harmonically and is given by F sin o)t for 
the fundamental, so in our notation the deflection due to the 
fundamental is given by 

F&inoit 

^ = 2~1 1 • (^) 

mp^ + /op -f T ' 

By the expansion theorem we obtain that 

( Jait g- at 

f f 

T — co^m + jap 2mp {ja — Pi ) {ja — p^ ) 

— P2)e^*— O — Pi)e"^‘]j 

where p^ and pg f^e roots of mp^ -f- pp -f t == 0, viz. 
— ad;/?. 

For recurrent phenomena we may allow the vibrator to 
remain in circuit for some time without registering, and so we 
are interested in the particular integral only, viz. — 

pjmt 

A = F . : — 

T — ma^ 4- jap 
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On rationalizing we obtain 
A = F 


sin {(ot — <f>) 


(46) 


+ (t — mco®)^] 

where tan <f> = pa)/{r — mo)^). 

In oscillographic work the disturbing force will be of the 

n **00 

form I„ sin {neat — /3„), where c is a constant converting 

naa 1 

the current into a force and n is the number of the harmonic. 
Under this condition the integral becomes 


"Y ( sin (ncot- /?„ — (/>„) ) 


• (46) 


and tan^n = pneoKr — mn^co^). 

Inspection shows that if m and p were zero, the trace would 
be an exact facsimile of the wave. As this condition is impos- 
sible of attainment in practice, the damping of the instrument 
should be as near critical as possible, or p^ = ^mr. 

Substituting, we have 


A = T -^nSin {ncot- p„ — <l>n) 
^ ^ ^ T + n^co^m 

n 1 ' 

where tan = 2no)\/{Tm)l{T — 

But when m = p = 0, the deflection is 


Ao = 



I„ sin {ncot — jg„) 

T 


. (47) 


. (48) 


The difference between Aj and A^ gives the error existing in 
the resultant trace. The ratio of the amphtudes is 

Ai T 1 ^ 

Aq ~ T + ~ 1 {n^co^mlT) ~~ 


The natural or free period of the instrument is given by 
Tq = 2-n\/{mlr) 

Also n^coy^TT = n^P, where / is the fundamental frequency of 
the phenomena under investigation. Let 1// = T = the period, 
then 


" “ 1 + nyTJTf 


( 60 ) 
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As this ratio is less than unity, there will be an error depending 
on the harmonic — the higher the harmonic, the greater the 
error in the amplitude. Also the jihase angle now becomes 

an 9n- _ n^TJT) ’ ^ ’ 

and it is seen that there is a small phase shift among the com- 
ponents. 

In the recording of transient phenomena we cannot allow 
the time delay necessary for the decay of the transient deflec- 
tion, so it is essential that free vibrations are non-existent or at 
least damped out as rapidly as possible. As the critically 
damped or aperiodic instrument gives the best results under 
steady-state conditions, we will consider the response of this 
instrument to transient conditions. In general, the force due 
to the action of current in the vibrator will be represented as 

n =0 0 

F — 2 FnS ~ sin nwt 

n = j 

and on substituting this in our force equation we get 


A = 






where tan 
and 





CO 


F„e ~ sin (ncot — (/>„) 

+ n^co\2/^-K)^ 

nco{K — 2^) 

+ //(//- A) 

= rjm, K — pl2m . . . (62) 


The first term of the above solution represents free oscilla- 
tions, which will die away unless A ^ 0. Here we require A 
to be large, or the value of co^ — AY4 = 0, i.e. aperiodic motion. 
So then we have that the first term will be given by 


FT (i) g( a+W (y £< ““ 

m Lc>2 + (/S - a)2 2f~ ~ co^ -f ^ - cuf 2p 

When fi approaches zero, we have 


Ay 


F 

m 


CO 

' co^ -f a2 




(53) 


16 . To Increase the Rai^e of an Oscillograph. As an oscillo- 
graph requires to be critically damped to give a true record 
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over a range of frequencies, the range of frequencies is restricted 
to the natural period of the vibrator. An interesting method 
of increasing the range of an instrument by means of a resonant 
shunt has been devised by Irwin* (Oscillograph) and Wynn 
Williams. t The circuit is as shown in Fig. 8. 



Fio. 8 


Now, as the current through the oscillograph is infinitesimal, 
we may write 

[(L, + L)p2 + + R)p + 1/C" + llG]q - V 

as the voltage across the whole circuit. 

Across the oscillograph strip we get 

[Lp^ + Bp + 1/C]q — e 

But e is proportional to the oscillograph deflection x, and so 
we write that 

[Lp2 + i2p + 1/G]q = X 

Thus we obtain 

V{Lp^ -\-Bp+ 1/C) 

"" ~ (Cl + L)p^ + (R^ + B)p + 1/Ci + 1/C 
Let be neghgible and write a = Rf'AL, = i?i/2C, 
= ^LG, and = 1/CCi. 


We have 


V{p^ + 2ap + (V) 

-j- 2(« 4- ai)p + (coq^ -f 


From this equation we see that the damping is increased to 
(a + a^) and the natural period to \/{(>}^ + a>i^)- If we assume 
aj = 4a and then it is seen that the damping 

coefficient becomes 5a and the natural period = 5a>o. Thus 
we would obtain an amplitude scale which is the same as if 
the galvanometer had a period of 5 coq and a damping constant 
of 6a. 


♦ Oscillographs. Pitman 192.5. 
t Phil. Mag. 11926), Vol. 60, p. 1 
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1. General Solutions. The operational notation leads to neater 
expressions when we are dealing with series and parallel circuits. 
To illustrate the method we will consider RL circuits. 

(i) For the series arrangement we have that 

^[^1 + 

and i\R^ + AgP] == 

On adding these expressions, we have 

JS = Cl + Cg = + i?2 d" (^1 “1“ -^2)p] 

.-. i = [E/{Ri + JSa )] [ 1 - e - <«> + . . (55 ) 

The time constant is given by {L^ + L^)I{R^ + iJg)- 

(ii) For the parallel arrangement we have 

(i?! + Lip)ii — E 
and (Ra + L^)h = ^ 

But — I 

I = {E/R^) [1 - £ - <«‘«/^‘>] + {EIR 2 ) [1 - £ ■ 

= E{{R^ + R^)IR^R^~ (l/i?i)£-<«-'/^>) 

- (l/i?2) . (56) 

As the general method is 
obvious, we consider now 
certain practical ax)plications. 

2. High-speed Circuit- 
breaker. To reduce the 
R 2 L 2 damage due to “flashing 

Fig. 9 over ” of rotary converters 

on electric traction systems, 
a circuit-breaker with a high speed of operation has been 
developed. For details of construction reference should be 
made to McNairy.* A simphfied diagram of connections is 
shown in Fig. 9, where R^, are the trip-coil constants, R^, 
those of the holding-in coil, and R, L those of the load. 

♦ MoNairy: Trans. A.I.E.E. (1926), Vol. 45, p. 962. 
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We now have that 


+ + hp] ■ • (67) 

and also i^-\- = I . . . (68) 

As we may neglect the impedances of the coils in comparison 
with that of the load under fault conditions, we get 

^ + 

From equations (57) and (58) we get 

h[-®i + -^2 + (-^1 + J^i)p'\ — [-^2 + L2p]I 
or i /(jgg + L^P) 

_ ^ ^2 + L 2 P 

R + Lp + iJg + (-£'1 + -£' 2 )^^ 

We now use Borel’s theorem (page 56) to evaluate this 
operator. Denoting RJL by a and (JSi + ^2)/(^i ■+■ ^ 2 ) ^ 

we get 

(p + «Mj> + ^)] ' 

But ri/(p + a)] . (7) = (1/a) [1 

and [pUp + o^)] • (7) ?= (see page 11) 

So the expression for the admittance is 

A{t) = AS ) . ^ 2 ( 0 ) + £ AS - . AS ^) . dX 

For the first term of equation (59) it becomes 
A{t) = l/(p + 1/a) {p + i5) 

- AS- X).AS^)dX 

= £ (l/a)[l-£-«<‘-'‘»]£-'’^dA 

On integrating we get 

(1/a) [(1 - e -/'‘)/^ - (e - _ e - «0/(a - 
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For the second term of equation (69) we get 
A(t) =pl{p + a) (p + /?) 

_ r e “ e ~ dA 
Jo 

Combining the results, we have 

L[L^A-L,]lcc\ /3 “ u-fi ) 


E rR^ 1-e-^' (R^ la 

L^ + lXr' ~\R l) 



Now, ii is the trip-coil current, and from equation (61) the 
rate of initial rise is given by 

[§L, = 

The current will attain its maximum when dijdt — 0 ; and 
after further working we get 

£(« - = a(RA - L^R)j{a.RA - ^L^R) 

So that the time lapse for the trip-coil current to attain its 
maximum is given by 

t = [l/(a - . Ign. [^{RA - L^R)l{aiRA - ^L^R)] (62) 

3. Impulse Generator. In routine testing of insulators, it is 
usual practice to apply an impulse voltage of wave shape 
approximating to that of a lightning surge. The circuit adopted 
is that due to Marx.* It consists of a bank of condensers which 
are connected in parallel for charging, and are discharged in 
series. A bank consisting of forty-eight condensers, each 

* Marx: E.T.Z. (1924). 
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charged to 75 kV, would give 3 600 kV on discharge. Due to 
internal impedance, a pressure of the order of 3 000 kV would 
be available. An elementary connection diagram is shown in 
Fig. 10, where is the lumped generator capacity and the 
capacity of the insulator under test and the capacity to 
ground of the whole equipment. 

Let the voltage across the insulator be e, and the voltage 
on the condenser bo E, Sut t = 0. Also, let the current 



Fig. 10 


distribution be as shown in the diagram. Let the volts across 
the condenser be e^. 

Then = (i? + Lp)i + c 

Also Cc = (1/C'iP) • (I — i) and (1/p)/ = Q 

Thus {R + Lp)i + e — • (/ — i) 

or {R Lp -\- 1/C'ip)i + e = QjCj^ = E 

But i e[(l + R^G^p)IR^ . (/) 


SO that 


p /? 2 /(l + B2G2P) ,,, 

^ R + Lp + \IG,p + R^I{1 + Rfy>) 


This expression may be put in the following form — 

jp PI^^2 l■l\. 

P^ + {RIL+ llR^C^)p^ + {RILC^R^ d- 1/L6'i ' ^ ^ ^ ^ 

+ \ILC^)p + IILC^C^R^ 

= E . — p--, {]) 

p^ + ap^ + op + c 

where A = l/LC^, 

a = {R/L + IjR^C^), 

b = RILC^R^ + IJECy + IILC^ 

and c — IjLC^C^Rz 
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As these coefficients are all positive, there can be no positive 
roots. For a third-degree equation, one root will be real and 
negative, and the other two imaginary and negative. So we 
may write equation (64) as 

Ap 


^ ^(i> + a)[(p-j-A)2 + a>2] 

Factorizing, we get 


■{ 1 ) 


AE 


(A — a)^ -1- |_p -f- a 


J 


P 


cop 


CO ' (p -4- co^ 

P^ + ^P , 


2 


(P + -f 0)2^ 
By the expansion theorem we have (see Appendix) 


e = 


AE 


(A - a)2 + 0)2 


. — at 


e - sin cot 


CO 


+ £“■** cos cot 


] ■ 


(65) 


The last two terms are oscillatory and should be eliminated, 
i.e. 0)2 must be 0 or negative. 

With o) = 0 we have 


^ p ' Z -f 0)2 (A— a)te ^ + £ . (65 a) 

As the value of is very large, we may consider the circuit 

for all practical purposes as 
RLC^G^, and from Chapter II 
we have that the critical resis- 
tance is given by 

R = 2[Vi^(l/Ci -f l/Ca)] 
This value of resistance must 
not be exceeded, otherwise the 
discharge will be oscillatory. 
In practice a wave is inspected 
by means of an oscillograph for ripples, and these should not 
exceed ± 5 per cent of the maximum amplitude of the wave. 

4. Mutually Coupled Circuits in Free Oscillation. Systems 
may be so coupled that a disturbance in one member is 
reflected into a second member. The reflection occurs electro- 
magnetically as in a transformer, or mechanically as in the 


Ri Rz 




SEEIES AND PARALLEL CIRCUITS 


39 


case of a turbo-alternator (Chapter VI). Let us consider the 
case shown in Fig. 11. At a time < = 0 let an impulse oscillation 
be set up in the primary circuit. Then we have the voltage 
equations as 

Z-ipii -f Ryi^ -f (l/Cip)!! -I- =0 • (66) 

Zfgpia + Rzh + iy(^2P)h + ^Ph = 6 . (67) 

Eliminating between these equations, we get 

i,[(L,L^- M^)p^ + (i?,La + R,L,)p^ 

-f {LJC, -f R,R, + LJC^)p^ 

+ + ^2l^l)P + ^ 

Adopting the notation — RJiL^, ag = RJ^L^, o)-^ — l/LiCi, 
0 )^ = l/LgCa, and K = we get 

- K.^)P*‘ + 2(ai + aa)^® -f (w^^ -f 4aaa -j- (a^^)p'^ 

+ 2(aift)a® + a2f«i®)p -f Wi^tUa®] = 0 (68) 

and an identical equation for ig. 

The determinantal equation is of the form 

op* -f- bp^ -f cp^ + dp + / = 0 . . (69) 

For a highly oscillatory circuit with low dissipation, the roots 
may be extracted by the method of Guillemin (Chapter I) 
as follows — 

D(p) = op* -f cp* + / 

and so p^ = [— c ± y' (c* — 4a/)]/2a 

0{p) — hp^ + dp 
— d = {bp^ + d)l{4:ap^ -f- 2c) 

The roots are given by 

be T bVic^ — 4a/) - 2ad IF c ^ y (c^ — 4a/) 1 

± 2y (c^ — 4a/) ^ •'/y L 2a J ' ^ 

As the roots are conjugate, denote them by — ± and 

— ya ± j^ 2 - Now, with conjugate roots the integration con- 
stants are also conjugate, so let them be given in the form 
+ and The solution is thus in the form 

ii = Ae-y** -f- 

+ Be- >'•< [£^<'*** + •” + £- + ®>] 

— cos (jffi< + A) JBi£“’'*‘ cos (jSa^ -1-5) . (71) 
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For the natural angular velocity we have that 
= [cT 'v/(c"-4a/)]/2a 

(coj^ _j_ 4a^a2 + T V[ + 4ai«2 — 


~ 2(1 - K^) 

and, on neglecting a^aa, we get 




a _ (V + 4(1 - K^)(Oj_^(Oi^] 


2(1 - K^) 


(72) 

(73) 


From this equation we see that — 

(i) The response angular velocity is not the same as for 
each circuit acting independently. 

(ii) When a>i = coa = !□, then the response angular velocity is 

P ± K). Thus the amount of difference depends on the 

value of K. For a graphical method illustrating the frequency 
variation due to change of K, reference should be made to 
Plummer.* When K == 1, /?2 “ when 

K — 0 then — O, i.e. the oscillations are of one 

frequency only. 

(iii) The expression for the damping factor becomes 

y = (6^32 + d)/( 4 a )32 + 2c) 

(ai + + (iXiCCg^ + «2«^1^) 

~ 2(1 — + (Wi* + 4aia2 + oi^) 

At resonance Wi = cog = Q ; then, on neglecting a^ag, we get 

«! + ag r /32 + 1 

2 l{i- 

With /3 = f2/\/(l ± K)y then y approximates to 
i(«i + «2)/(l i 

So, corresponding to the two values of there are two values 
of the damping coefficient y. 

We now have that 


' ii = 4ie->'‘*cos(/3i« + Ai) + ^iC-y-'ces (^gi + dj) j 
and tg = >4ge~^‘‘co3 {^xt + Ag) + J?ge“’'*‘cos (/Sgi + ^g)) 

where AxA^BxB^ are arbitrary constants. These may be shown 
not to be all independent. 


♦ Plummer: Phil Mag. (1927), Vol. 34, p. 612. 
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On neglecting R^, R^ and substituting in equation (66), 
we get 

^ ("i®— ^ 1 ^) sin (fiit + Ai) + ^ (a>i® — (52®) sin + ^i) 

M M ^ 

- ^2/3i 2 sin -b Aa) - sin = 0 (76) 

■L/iPl ^lP2 

Equating terms of like frequency, we get 

^i(ct)i®— /^i^) = {MfLi)A 2 f^j^ and A^ = Ag = A 
(i/) = (ikf/Ljj^jjSa® s-nd = (Jg = <5 
So that we have 

cos + A) + cos (/Sgi + ^) 
ia = LJM l]^i cos (^^^ + A) 

+ j^(^) “■ cos + ^)) 

When damping is considered, then the equations may be 
written as 



— A^e cos {Bit "f" A) “h B^e cos (B^t -f- d) 

andia = ^ViLJL^) {[{(aJBif— l]^ie-’'-‘cos (Bit + A) 
— [1 — cos (B^t + 5) 


(77) 


The general character of the oscillations may be deduced from 
equations (77). It is seen that there are two damped harmonic 
oscillations in each circuit, which leads to a periodic rise and 
fall of the amplitude corresponding to the acoustic beat between 
two musical notes. 

Further, as while /?2 > it is seen that the low 

frequency oscillations are in ])hase in the two circuits, while 
the higher frequencies are tt radians out of phase. Thus the 
beats will be out of phase with the energy surging back and 
forth between the two circuits at a beat frequency of 

If we assume at < = 0, — 0, and that the charges 

on the condensers are and = 0, it is easily shown that 

for the case of tuned circuits that 


a 


V(1 +^) 




• Yit 


Q. 


sm 


V(1 + K) 


^ i i 


Q 


V{l-K) 


4-- (T.152) 
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and 




a 


\/(i + 


?0« 


■y,t 


Q 


Sin 


-\/(l + 
<1^ 


t — \ 


»'** sin 


D. 


V(i-^) 

Q 

t 




So the effect of the couphng is to lower the frequency of the 
gravest tone and to raise it for the other. 

5. Coupled Circuits in Forced Oscillation. When there is a 
voltage Eq sin (ot applied to the primary side, the circuit 
equations are 

+ i?!*! + (l/^iP)h + ^ 

L^pit 4 - + Mpix — 0 


We will here investigate the steady state only, i.e. we assume 
that the circuit has been made long enough for steady-state 
conditions to obtain. The currents and will be sinusoidal, 
and so we get and 


Hence [Ri + j(<oLi — l/o)Ci)]Ii + jcoMI^ = E 
and [iij + j((oL^ — llioG^)]!^ + jcoMIx = 0 
or Z^Ii + — E 

and Z2/2 + j(oMIi = 0 

c J.I. ± T ^ _j r —jcoME 

So that Ix = _p corny - z^Zx + coWV^a]* 


As the actual current is the imaginary part of and I^, so 
setting 

= Zx-\- co^MyZ^ 

Then/i= {EIZx')e-j^ 
or ix — (EJZx) sin {oit — ^1) 
and /a = - (joMEjZ^'Zx' )e 
or 1*2 = {cdMEIZ^'Zx) sin {cot — — <^2 ~ 

where = V[(J?i + + (Zj - 

tan <l>x = (Xx- (o^M^XJZ 2 ^)I(Rx + (om^R^IZ^^) 

Z, = V(JS2* + X,^) 

tan (f>2 = X2IR2 

6. Testing Transformers. Now let us consider the case of a 
transformer with a bushing of conductance g and capacity C 
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as burden. The circuit is as in Fig. 12. On neglecting the 
primary resistance and transferring all quantities to the 
secondary side, we get that* 

— ® • • (78) 

But *2 = (gr -1- Cp)e . . . (79) 



and, on dropping subscripts, we get 

E = {R alyp) (g + Cp)e + e . . (80) 

where e is the voltage across the bushing and a = 1 — 

Now E will be of the form E cos {oit — 0), where 6 is the 
“switching-in ” angle, so that with aL — L 

E cos {cot — 0) 


LC\jf>^ H- {gjC + RIL)p -f (1 + Rg)l{LC)] 


. (7) 


(81) 


Let the roots of the determinantal equation be given by 
(— a ± ^5). Then, by the expansion theorem, we have, for the 
steady state, that 

E cos {cat — 0 — (f)) 

Cq = 


coCV {/?2[1 -f- gL/CBf + [(1 -f Rg)l(oC - coLf} 

. , ^ , R[l+ gL/GR] 

where ton ^ “ (1 + ' 

The transient term will be given by 


(82) 


■"2 


where Z'{p) 
Let 
and 


j^ijco— p„)Z'{p„) 

: [2p ± {g/G R/L)]LG = ± 2LG^. 
a = gl2G + i?/2L 

= ^/i{g|2G + Rl2Lf- (1 + Rg)ILG] 


( 83 ) 


* Russell; A.O. Theory, Vol. II, p. 260. C.U.P. 1914. 
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Then we get 

cyTna J 




-Pt 


2LCP Ijo) + a — jw + a + /5J 
Ee~'^ e~ r {jco + a) sinh ^ cosh /Si* 




2LCp L (a2-^2_ ^2) ^_^-(2aa>) J 

On rationaUzing and discarding imaginaries, we get 
Ee-’** 


(84) 


«< 


LCp 

cos (6 + y) . (a sinh + /S cosh j8<) + c) sin (0 + v) • sinh/3< 


[ 

2coa r 

where tan v = f, [ 


V[(“^ ~ ~ o)^)^ + 4a W] 

1+gL/CR 


(1 + Rg)IC(o — (oL 


] 


] 


(85) 


— ^ tan <f> 


When the resistance R is small, then the value of will 
approximate to jl\/ (LC), so that tan ^ = 0 or ^ = 0. With a 
high natural frequency we may neglect terms containing w//3 
and a//S, so that we get 

e = ^o[cos {(ot — 0) — e ~ cos 0 . cos . (86) 

The maximum voltage will occur when the switch is closed 
at (ot = 6 — 0, i.e. at maximum applied volts. Further, due 
to the value of circuit constants, a voltage higher than requisite 
may be applied unwittingly to the bushing. 

7. Arc Shunted by Resistance. As the current in an a.o. arc 
must follow the voltage fluctuations, the arc will be extinguished 
and restruck periodically. The phenomenon will depend on the 
rate of building up of the voltage across the arc space. We 
consider a circuit as shown in Fig. 13. Let the current through 
the inductance be the current through the capacity be 
and the total current be I. For a voltage applied to the circuit 
we get, when the arc current is zero — 


E = RI+ {llC-p)i, = RI + q,IC 
I = iL + ic = H + Mo 
pi = pi^ + p^q. 


. (87) 


Combining these equations, we get that 

^ = (lo/C) + R[ii, + 
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On dilferentiation this becomes 

0 = p{<lclG) + Rpi,. + Rp% 

But Lpi^ = qJC 

Rp\ -f {l/C)pq, 4- (RqdLC) = 0 

L 


R 


Fig. 13 

The roots of this equation are given by 

p^- 1/2GR ± V(l/402i?2 - l/LC) 

so that if IjAG^R^ < l/LC, or R is greater than the critical 
resistance, we get as solution — 

q^ = cos (ot B sin m) . . (88) 

where oi = ^/(\|LG- 

Thus the voltage across the arc is 

Ca = E — qJG = E — B~ [A cos cot B sin cotyG (89) 

The constants may be determined from terminal conditions 
at < = 0, as, for example, e,, = E^ and — — /j, where I is the 
current at which the arc fails. So that E^ — E — AjG. 

The condenser current is — pq^, so we have 
-I^ = B(o- AI2GR 

A = G[E - LJ and L - [- /j + (L - Li)/2i2]/tG 
and ea = E — s~ [{E— E^) cos cot + (1/wC) 

[{E — Ei)l2R— 1 1 ] sin cot . . . (90) 

From similar considerations it may be shown that if 
1/LG < 1/4L^C®, i.e. if the resistance is less than the critical 
value, the voltage across the arc is 

e, = L- Li) cosh^i + (1/Ci8)(^^^- a) 

sinhjSfJ . . . (91) 
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On differentiating these expressions, the rate of building up 
of the voltage across the arc space may be determined. We have 

deJdt = [IJCl^, . . . ( 92 ) 

Thus it is seen that the initial rate of voltage recovery is 
inversely proportional to the capacity of the system and 
directly proportional to the current at which the arc fails. 



CHAPTER IV* 

BOUNDARY CONDITIONS 


1. Mathematical Theorem: Subsidiary Equation. So far we 

have considered cases where the circuit was dead at a time 
< == 0. As this condition does not always hold, we will now 
give a subsidiary equation which in combination with the 
standard form of the expansion theorem gives the solution 
under these conditions. 


Let the general equation be given by 
Ti/i\ , d”~^y , d’^~^y , 


+ ««y 


dtn ' ' ^dt^- 

= jp”y + p"" V d- «2 + . • • + • (93) 

At < = 0 let p''~^y = c^, p^~^y — c^, and so on. Also let 
p^y = u{t). 

Then by integration we obtain 

P’'~ V = ^u{t)dt + Cl 
Repeating this process, we get 

p”~^y = + Cit + ^2 

= p-^u{t) + Cit + Cg 

So y = p~”u{t) + Cit”'~ ^l(n— 1) ! + c^V'~'^l{n — 2) ! 

+ . . . + c„ . . . . . (94) 


Substituting in equation (93), we have 
F{t) = u{t) + aip~^u{t) + a2P~^u{t) . . . 

+ + <*2(^1^ + C2) . . . 

= u(t)[l + aip-^ + . . . + a„2>“”] + CiUi 

+ 02(^1^ *^2) •••■!" V(^ ^)- 

+ C2<'*-2/(w- 2)! . . .) .... (96) 


By equation (94), we have 

u{t) = p^[y — Cif'-''^l(n— \)\ — c^V*-^l{n— 2)\ . . . — c„] 

So by equation (96) we get 

F{t) = [p” + + . . . + a„] [y—CiV'-^l{n— 1)! — . . . 

— c«] + + . . . + a„[Cit”-^/{n— 1)1 + . . .] (96) 

* Reproduced by kind permission of Editors of Philosophical Magazrne, 
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Let /(c) = 1)! + — 2)! + . . . 

Then !>"/(«) = • • • + 

and aiP"~y(c) = {a^c-y + + • . .) 


Thus equation (96) reduces to 

F{t) = [p» + • . • + (Inly 

- ClP - C2{p^ + «l2>) + • • . 

• _ F{t) + c^p + C2(j)^ + «iff) + • • • 

p” + a^p ^- 1 + a^p^- 2 . . . + a„ 


We will demonstrate the use of this equation by a few 
examples. 

2. Breaking Circuit. Letusconsider an i2L circuit. Initially, 
the current will be / at i = 0, and the voltage equation at 
break is given by Ri + Lpi = 0 or pi + {RIL)i — 0. From 
equation (97) we have that F{t) = 0 and = I. So that we 
get 


P 


P + RIL 


[ 1 ) 


On expansion this becomes 

i = Ie- 


For the RLG circuit we have the following voltage equation — 

Ri + Lpi + (1/C'p)i = 0 

or p^q + (lilL)pq + {^/LC)q = 0 

At f — 0, q — Q, and we have that F{t) — 0, = R/L, 

ttg = If LG, Cl — 0, and Cg = Q as the constants of equation (97). 
Thus equation (97) now becomes 
. ^ Q{p^ + Rp/L) 

^ RpjL + 1/LG ^ ’ 

By the expansion theorem this becomes 

^ _ p r CPi + , P2 + RIL)e^>n 

^ ^ I 2pi + R/L 2p, + R/L J 

The values of the roots are 

p = — R/2L ± V — ^LG) = — a db Wq, say. 
Hence q — Qe~‘‘‘[a, sinh co^t + coq cosh cO(,t]/o}Q 
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For oscillating conditions we have coq = jo), say, and 

g' = [Q . . GOs{a)t — i>)]/co V (a® + o*) 

where tan (f> = a/co. 

And pq = i z= — (Qs'^ sin cot)lcOfj{x^ + coq®) 

3. Field Discharge Resistance. When the field circuit of a 
machine is broken, the energy stored in the magnetic field 
must be dissipated so that the sparking at the switch contacts 
is minimized. Usually this is 
done by opening the circuit 
through a discharge resistance 
(Fig. 14) which is inserted in 
the circuit before the main circuit 
is broken. So we have that 

{R + + Lpi — 0 

i.e. pi {R R 2 )ilL — 0 Fig. 14 

At f = 0 let i = 7. So, neglecting the presence of iron, we get 



P 


P -^ 2 )/-^ 


{ 1 )= + 


(98) 


The energy stored in the magnetic field is given by ^LP and 
the energy dissipated in the circuit by resistance is 

(R + R2)Ps-^^ 


So that the net energy in the spark will be 
I\\L- (i?i + i?2)£-2«‘] 

For no sparking we must have a time of switch opening 


1 , L 
2(R + R 2 ) 


sec. 


(99) 


4. Automatic Traffic Signals. A further example of the 
discharge circuit is in traffic-signal systems. An elementary 
diagram is shown in Fig. 15, where A is the relay operated by 
vehicles having the right of way and jB is a demand relay for 
vehicles wanting the road. With the relay A on the upper 
contact the condenser is charging, so that 


E 

+ P^ 




_ e - tICR, 

R, 


. ( 100 ) 
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and the voltage across the condenser is 


f* (ilC)dt 

Jo 


(101) 


When the relay B closes, the voltage across the neon tube is 

E-E[l- . . (102) 

which is insufficient to cause the neon tube to fire. After a 
definite time the relay A is moved to the lower contact by a 

— . master timer and the condenser 

JL I discharges through R^, giving 

O H = . (103) 

^ The voltage across the condenser 
‘-T-' is illCp)i 2 — and across 

^ / y the neon tube is E[1 — 

I , When this voltage is equal to the 

critical volts e^, the tube fires and 
"T" 5 Ji2 time of discharge is 



Fig. 15 


(7J?2 Ign. [E/{E 


The discharge of the neon tube 
actuates a relay, which in turn operates the signal light. The 
size of the condenser is the hmiting factor in the time equation, 
and should be such that it will not allow _ _ 

the neon tube to fire before the last ' TfHlf y 

vehicle is clear of the intersection.* 

6. Opening of a Circuit through a \ * j. 

Condenser. Frequently condensers are f-Sj / 

placed across the contacts of a switch '''p I . T 

to reduce sparking. In Fig. 16 is shown I | 

an LR circuit with a condenser across jq 

the switch contacts. Initially, the cur- 
rent with the switch closed is steady and equal to — EJR. 

On opening the switch the condenser is inserted. So, neglect- 
ing arc phenomena at switch contacts, we get 

i 

~ LC[p^ -f RpJL -f 1/LC] • ^ ^ 

At t = 0, i if,orq — 0 and pq = Thus we get and so 

„= + 

^ p^ + Rp/L + 1/LC - ' 

* Priest: Journ. I.E.E. (1936), Vol. 77, p. 149. 
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• ( 104 ) 
. (105) 


where tan <f> = to/a. 

At < = 0 the current through the condenser = E/R. 

Initially all the current passes through the condenser and 
none through the switch contacts, and thus the voltage drop 
across the contacts is small. 

6. Bestrikiug Voltage subsequent to Circuit Interruption. 

When a circuit-breaker opens, an arc is estabhshed between the 
contacts. On an a.c. circuit, every time the current wave passes 
through zero the arc is extinguished for a period of 100 or so 
microseconds. Should the dielectric strength of the gap be 
re-established at sufficient speed, then the voltage across the 
gap will be insufficient to cause breakdown. The arc would 
then be extinguished. Consideration should also be given to 
the effect of power factor. For example, on short circuit the 
current lags by 90° approximately on the voltage. When the 
current is zero, the appMed voltage will be at its maximum 
value. Thus we have the maximum generator volts maintaining 
the arc. The voltage at the moment the switch breaks may be 
termed the recovery volts. 

Oscillograms of the current under short-circuit conditions 
show high-frequency oscillations. It would appear that the 
ability of a breaker to interrupt a circuit depends on the 
increase in the dielectric strength of the gap occurring at a 
greater rate than the rise of these high-frequency oscillations. 
This voltage is termed the restriking voltage, and may be 
calculated as follows. 

The opening of a switch is equivalent to the introduction into 
circuit of an e.m.f. equal in magnitude but opposite in sign to 
the potential difference appearing at the switch. Let Z{p) be 
the impedance operator of the system and I the current prior 
to the opening of the switch. Then the voltage across the 
switch contacts would be I . Z{p) = e. This is the voltage 
necessary to prevent current flow. To calculate Z{p) we have 
to take into account the circuit constants, including the arc. 
For a power arc the drop is directly proportional to current X 
arc length. Assuming that the switch plunger bar reaches its 
final position instantaneously, then the arc drop is directly 
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proportional to the current. So for one phase tlie circuit may 
be represented as Fig. 17, where RLG are the circuit constants 
(including those of the machine) and r the arc resistance. 
We get 

® ^ ^ ■ LCrp^ + {RCr L^p + R + r ’ 



As e is the voltage across the switch contacts it may be 
termed the restriking voltage. If a — {RCr + L)l2LCr and 
= (72 + r)ILCr, then we have 

i> = — a ± \/{<x^ — too^) — — a ± Wi, say. 


’[ 




From the expansion 

rR . , i2£““‘ 

+ ITT^ smh (o,t + 


R + r'^ co^C 


a)i cosh a)it + a sinh coit \1 
a? + j\ 


(107) 


For an a.c. system we have that 

r r{Lp + -R) 1 /Yx 

^ ^ iLCrp^ + (RCr + L)p + R ^r ' 


The steady-state term is 


IrV 

^'= 2 1“ 


(jR jo>L)e^'“* 


(o^LCr -f- {RCr L)j(o + JS -f- r 
and the transient term is 

p„(R + Lpn)e^’‘* 


] 


^'2 


(Pn^ + (o^){2LCrp„ RCr -f- L) 


V Pn{R + LpJe^»* 
^ ^ 2I/yr(Oi{p„^ -}- o)®) 


(108) 


. ( 109 ) 
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Now, Pn is usually very large compared with co, and so we 
may approximate — 




^ [R + Lp„)e^o* 
^ 2LCro)iP„ 


( 110 ) 


On neglecting circuit resistance and damping, we get 

2LCo)i ' = ~c^ 

Had we considered the circuit as having distributed con- 
stants, then, as will be seen in Chapter VII, there are several 
free periods, and the expression is in the form 


e 


n—« 

e. + 2 

n = l 


InS “ 

sinh (o„< 

C(o„ 


( 112 ) 


Further investigation on these lines may prove of some 
value in the rating of circuit-breakers for a definite system. 
The difficulties respecting the calculation of the natural periods 
of a system are not insuperable. 
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CHAPTER V 

snnrsoiDAL asd damped voltages 
OTHER METHODS 

We have already demonstrated the expansion method of 
dealing with sinusoidal voltages (Chapter II). When we come 
to deal with more complicated voltage expressions, that 
method becomes unwieldy; but we are able to make some 
simplifications. These are — 

1 . To express the voltage function as its equivalent operator. 

2. To use the Heaviside shifting method. 

3. To use the superposition theorem or the Duhamel integral. 
1. Equivalent Operator of the Volt^e Function. This method 

is based on the following mathematical argument — 

cos (at = 

P 

But expressed operationally is ^ • (^) Appendix), 

ande-^"‘ is — p-r- . (2). 

Hence 

coa cat = if— ^ + —^1.(2) = X - ..{1). (113) 
In like manner we obtain 


sin (ot -o-r — 2 . (2) , 

Also cos (co< + jS) = cos cot cos P — sin (at sin /S. 
Hence, on substituting equivalencies, we get 
. P* cos d — cop sin d , 

«» — - ■ w 

For the hyperbolic functions we have 


sinh (ot = 5 . (i) 

and cosh cot = r . 11) 

jps — ^2 ^ / 


(114) 


. (116) 
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2. The Shifting Method. This method was introduced by 
Heaviside,* and he demonstrates it as follows — 

= e“ (dtt/di) — e~ “* att = e “ — a)a . (116) 


It is thus seen that the rule for changing the sequence of p 
and e ~ is that on moving p from the left to the right of e “ 
we must change the p to (p — a) ; on moving to the left we 
change p to p -{■ c>.. Actually we are multiplying /(;p) by either 
g - ap or e + and so delaying or advancing the time of applica- 

tion of the voltage as shown in Fig. 45. Reversing the equation, 
this theorem may be written as 


-e( 


' Z(p) 


{i) = 


Zip + a) 

1 p 


6-“*. (I) 

ii) 


Zip -f a) ■ p -f a 


• (117) 


Thus, for the expression e sin o)< . (I), we have that 


sin tot 


top 

~ p* -f O)® 


so that 

sin (ot ==■ e~^ . — « • (^) 

p® -f <0® ' ' 


(»iP + 

ip + Pf + w* 


i.e. sin tot 


(»{P + /?) P ... 
(p -h /?)® -f o)® > -I- /3 ■ ' ^ 


e-'”. (i) 

. ( 118 ) 


Similarly 

e " COS (ot 


{n 

(p -f /3)® -f o® • ^ ^ 

pip + ^) ... 

(p + |8)® + «>® • ^ 1 


. (119) 

. ( 120 ) 


A further example, taken from line theory, will suffice to 
show the method. The operator is 

e->h^[ip-\r pY-o^-\.il) . .( 121 ) 


Multiplication by e “ . e'’* does not alter the value of the 

operator, so we get 
e~ ^ e~^l” \/\_ip -f- p)* — u®] 

= s-ptg-xlv^ipz_„i)^^il) .( 122 ) 


* O. Heaviside : Electromagnetic Theory^ Vol. II, p. 293. 
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3. Superposition Theorem or Duhamel Integral. This method 
is used by Carslaw in his Theory of Fourier's Series. 
Fundamentally it is a “step by step” method, and may be 
derived in the following manner — 

At a time t let an e.m.f. E{t) be applied to a circuit whose 
admittance is A(t). Then the current i{Jt) = E(t ) . A{t), where 
A{t) is Y(p)IZip). 

If now a voltage E(k) is applied to the circuit at a time 
t = X,we then have 

i{t- X) = E(X) .Ait~ h) . . (123) 

as the additional current. So we would get a total response 
(Fig. 18) of 

E{0).A(t)-\-E{X).A{t-X)-\-E{22.).A{t-2X) + . . . (124) 


jS 



If E{X) — AE is applied at a timj At, then the above 
equation is 

E{0) .A{t) + AE . A{t— At) + AE A{t— 2At) . . . (125) 
or, in the limit, we get 

l{t) = E{0 ) . A(<) + £ A(< - dt) . (dE/dt) dt 

i.e. I{t) = djdt I* E{u) . A{t— u)du . . .( 126 a) 

Jo 

or — dfdt I* E(t— u) . A{u)du . . . ( 126 b) 

Jo 
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By using the well-known partial integration formula 


[-1=. 

f * XdY + C YdX 

Jo Jo 


and noting that, if ^(t) = 

f f(t,A) dX, then 


1 

II 

dt/\ ddt 

-/(Mo)i+/(M.)^ 


we may transform our equation into one of the following 
forms, the exact form to be used in any particular analysis 
depending on the problem and the ease with which the mathe- 
matics may be handled. 

I(t) = E{0) . A{t) + 

« 

A{t — u) . E'{u) du 
'o 

. (1260) 

= A(0).^?(0-f J 

[* E(t — u).A'{u)du 

0 

. (126d) 

= A^(<) . A(0) J 

I* A'{t — u) . E{u) du 

0 

. (126e) 


^( 0 ) .^(0 


u) . A{u) du . (126F) 


4. Applications of the Above Methods. We will now show 
some applications of the above methods. As all equations 
must lead to the same result, we will give only one example 
of each method. 

Method 1. Let us consider the sudden short circuit of an 
alternator. In an elementary way we may represent the 
machine by resistance and inductance of constant values and 
neglect all other effects. When the machine is suddenly short- 
circuited by an ideal switch, we have 

Ri -j- Ijpi = E sin {mt + a) . (I) 
or i = ^ ^ sin (ca< + a) . (i) . . (127) 

Now, using equation (114), we get 

E (op cos a -f p® sin a , 

* "" RA-Lp ' p2 


5— (T.152) 
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The impedance operator is {R Lp){p^ eo*), and it will 
have roots — RJL and i joi respectively. 

On applying the expansion theorem (Form 1) we note that 

dZip) = . z^{f) + Z^'ip) . Z,(p) . (128) 

where Z^ = R Lp and Z^{p) = p^ ft>®. 


On substitution of the roots obtained from Z^ip) = 0, it is 
seen that the first term will cancel out, leaving the second 
term Z 2 {p) . Z^{p). Similarly, on the substitution of those 
obtained from Z^ip) = 0, the first term only will remain. 
(This is an alternative to Borel’s method. Chapter I.) Thus 
we get 


E 


~3l 


cos a — CO* sm a 
2o)^{R + jcoL) 




■ jto* cos a — CO* sin a 
— 2(o^{R — j(oL) ^ 


jcot 


0) COS a — sin a 

^ L[{R/L)^ + CO*] ^ 


(129) 


Let tan <f) = coLjR and Z* = .R* + to^L*. Then we have 


.Rfcos (a — ^) + j sin (a — <^) . 


=i[ 




Ja>t 


cos {oL — <l>) — j sin {a. — (f)) 

_ - 




sin (a — <l>)e 


= {EIZ)[co8 (a — <f>) sin cot — sin (a — <f>) cos cot 

— sin {(X — (f>) . e~ 

= {E/Z) [sin (cot + <^ — a) — sin (a — ^) . e " , (130) 


In this case coL will be very large compared with R, so that 
<f> approaches 90°. The equation then reduces to 

i = (EfcoL) [cos (cot — a) — (131) 


If ^ = a, then we have i = (EfcoL) sin cot, i.e. there is no 
current transient term when the switch is closed at a ~ Tr/2 
or when the supply voltage is a maximum. 

For the maximum transient the switch must be closed when 
cos a is a maximum, i.e. when a = 0, or the supply voltage is 
zero. The subject will be dealt with in more detail in Chapter 
VI. 

Methods 2 and 3. For the case when the applied voltage is 
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of the form Ee ~ ^ sin mt, we apply the shifting method first 
and then use the Duhamel integral. Thus we have 

E 

» = / ? ~r sin (ot . (i) 


R Lj) 




Now, applying the Duhamel integral, we get 

= [l-eG -««]/(/? _AL) 
where a = RfL 

^(0) = 0 
E{u) = E sin (DU 

So that we have 

e» . i(t) = ^(0) . E{t) + E(t-u ) . A'{u) du 
E r‘ , 

~l \ sin co(<— «) dM 

The value of the integral is 


( 133 ) 


So that 


.(y* + (A 


Lie 


a»^ + (A — a)® 

[(A — a) sin cot co cos cot] 


L 1(0^ + (A - a)® + (A - a)®] 

sin {cot + ^)J 

where tan <f> = ft)/(A — a) = coLI{XL — R) 




^/[(o'^L^+{XL- Rf] 

[e ■" sin ^ sin {cot + ^)] 


( 134 ) 


(136) 


(136) 


Method 3. For the application of the third method we 
choose a circuit {RL) to which a damped wave is applied of 
the form Ee ~ 

Now, for a wave given by “ ■** we get 


(137) 
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For a constant e.m.f. the equivalent of l/(iJ + Lp) is 
[1 — e~“*]/i2, so that 

A{t) = [l-£-“<]/i2 

and E{t)^E£-^ 

Thus we have -4(0) = 0 and A'{u) = (a/i2)e“““. 

Then I{t) = A{0)E{t) f E{t— u) . A'{u)du 


or 

7(0 = 

J* (Fa/T? )£-““. e . 

. (138) 

and 

/(<) = 

(Ecf.lB)e-^ J* 



== 

(77a/72) . £ - ^[e - - 1]/(A - a) 


i.e. 


[-fi;/(AL- 7?) [£-“'- £-"“] . 

. (139) 


6. The BLC Circuit with Applied Damped Sine-wave Voltage. 

In radio work, oscillatory e.m.f.s are impressed on RLG circuits. 
If the applied voltage be .£e“^*sin o)t, then the voltage 
equation is 

Lpi 4- -R* + = Ef.~ sin cot . (1) 


or 


EC 


LCp^ + RCp + 1 * 

sin cot . (i) 


, -At 


(140) 


By shifting, we get 


ECe 


■u . 


p— K 


LC(p - A)2 + EG{p - A) + 1 
sin cot . (1) . . (141) 


Dropping the term ECe ~ for the time being, and applying 
the second form of the expansion theorem, we get 

{jco— X)e~ (i>i — A )£ “ *’■* 


y 


LC(jco — A)® + RC{jco 


A) + 1 2LCcoi{jco- 
iPi- A)£-^-* 


2LCcoi{jco — pg) 


Pi) 

• (142) 


where P 1 P 2 = « ± j(^i> a = A — B/2L and coj^ = B^jAL^ — IJliC. 
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So, for the steady-state term — 

— (A cos (at -f CO sin cat) -\- j{(o cos cot — X sin (ot) 
y* ^ - (o^LC - j(o{2XLC -f BO) -f {LCX^ - RCX + 1) 

— sin(cof + 0) + i cos{o)t -f 6) 

= J) ^(^2 ^ ^2) . . . , 

where tan 6 = Xfo). 


(143) 


On rationalizing and dropping imaginary terms we get 

— cos(coi -j- 0) . [2coAiC' -|- EC(o] — sin(co^ -f 0) 

[LCX^ - RCX + 1 - (o^LG] 
+ ft>2) [{LCX^ - LC(o^ - RCX -f 1)2 

+ o> 2(2ALC' -f RC)‘] 

xxr..,. . , (oC{2XL + R) 

With tan <f> - _ ^ 2 ) _ jicx -l 1 

— sin((o< 0 (f>) ci4.a.\ 

^/{X^ + (a^).^/\D^\ ^ 


For the transient term we get, on substituting for p — 

= r '"- . ^+>1 ,-Me] 

2jLC(Oi Liw — jcoi — a jo) + jcoi — « J 

e“* ri(aA — a2 -+-^aco — co^) sin (o^t — co(co -f-^’A)cos to^^ 

jcOiLC [_ — 2Jcoa — (co^ — coj^) 


Vt 


Discarding imaginaries, we have 

(a® -)- coj2 — co2)[a(A — a) — coj®] sincoi^ 
-f- (Oi [A(a2 — co2 (Oi^) — 2co2a] cos (Ojt 


€-<*t 


(OiLC 


[«> 


-f COi^ — to2] -f- 4a2co2 


(146) 


Then the complete solution is given by 

i = ECe-»[y, + y,] . . (146) 

On neglecting damping, a usual procedure with radio 
circuits, we have coq* = lijLC and = A d: j<^o> for the 
steady-state term, after some simplification, we get 

coo* Biji{(ot + 6 + (f>) 

“ VK + ' V[(^* - £0* + V)® + 4A2co2] • ' 
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The transient term becomes 


Vt = 


t.e. = 




2jo}^C 


Lj(a) — <Uo) — 


^ + tyo)- 


d 


£(7 


- (A*— CO* + <Uo*) 

[coq sin a)o< + A cos coq^] — 2A(o* cos coq^ 
(A* — ( 0 * + coq*)* + 4A*co*) 


( 148 ) 


on rationalizing and discarding imaginaries. 

If tan 02 = ^/<yo tan ^2 = 2(oA/(A*— co* + coq*), then 
cos <^2 • sin (cOfli +02) 


r 

y*-Lc[- 


v'(A* + coo*) 
where A = \/[(^* — co* + coq*)* + 4A®co*]. 
So the full expression for the current is 


0) sin ^2 oos (o^t 


]1 


ECe-» 


* = 


t 


'e ^ I (Oq^ b 

A~ ( 7(AM^ sin (co< + 0 + ^ ) + ^ 


cos <j >2 . sin (coo< + 02) 

+ 0>o^) 


CO sin ^2 cos (Offt ]) • (149) 


6. Theory of Ballistic Galvanometer. Usually the action of 
the ballistic galvanometer is explained by considering the 
charge as given instantaneously to the galvanometer. We now 
investigate the movement when the charge is given over a 
period of time, as for example during the Recharge of a con- 
denser, when e == Effi ~ = E^e ~ From the theory of 

instrument movements, we get 

wp*0 + ppB + T0 = F . . (160) 

where p is the damping coefficient due to frictional effects 
and to the coil movement in the magnetic field. This latter 
effect may be expressed as O^r, where 0 is the galvanometer 
constant and r is the resistance of the coil circuit. The 
deflecting force 




Oe 

B 


= 


1 


■ ® = a. . 

mp* + pp + T 

By shifting, we get 
e = Ke-«* 


. e 


-at 


• {1) 


m(p — a)* + p{p — a) + T 


(1) 


. (161) 
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Using the expansion theorem, we get 
P T <<>0 

where a>o = V (p® — ^mr) and Z'(p) = i 2m(Off. 

r 1 £«</ e-(P + «.)l/2m 

0 = Z«-“‘ — 2 ; h— o T — 

Lma^ — pa + T Wq V 2ma — p + 

P - (p - a>«)e/2f)t \ "I 

. (162) 

2ma — p — tt>o / J 

The maximum deflection will occur at some time when 
dd/dt = 0, or 

« 1 j [(p + coo)/2m]e-<^ + ->>‘f^ 

moc^ — pa + T coq \ 2ma + co^ — p 

ma — &)(, — p / 

so that 

i: 


moc^ — pa + T 




2ma«>r 


£ - pt/ 2 m |-£«>,»/ 2 m _ g - «>.J/ 2 inj ( j 53 J 


For a critically damped instrument this reduces to 
OF t 

®"“= 1^:2^ ■ ■ ■ 


(164) 


7. Operational Methods of Evaluating some Integrals. The 
solution by the Duhamel method calls for evaluation of 
integrals of the form 


y = e 


ru =-l 
a. 

Ju - 0 


e““ sin CDU . du 


The integral term may be written down operationally* as 

♦ See Appendix. 

ix -i • (^) 

Shifting e““ to the right, we get 

i (ojp-ct) 

P ’ (l> — a)* + ft)* 

nr 1 0>{p-a.) p O) 

^ ■ (p — a)* + ft)* * ^ — a ‘ ' (p — a)* -f ft)* ■ ' ' 
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Evaluating this term by expansion, we get 

CO ^ ■ 

^ ~ ^ + «>2 2jco((x -i-Jco) 2jco(oc — Jco)^ 


L«^ + ' «» + ft>5 

Let tan <f> = co/oc. 


(a sin cot — CO cos cot) 


sin(to< — <i) + e““‘sin<i 

Then y = — ^ 77^, r\ 

V (“ + «> ) 

If the expression is given as 


/•« «» I 

y = e-“‘ 

Jtl - 0 


e““8in (com -{-d)du 


then, proceeding as before, we get 

1 po) cos 6 4-®* sin d 

p p* + co’* 

By shifting, we get 
(p — a)co cos 6 + (p — «)“ sin 6 p 
p ’ (p — a)^ + co^ " p — OL ' ^ 

(JO cos 5 + (p — a) sin 6 

~ (p— a)® + CO® 

By the expansion theorem we get 


c® + CO® a® + CO® 


(a sin cot — co cos oot) 


4 sin 6 a .”2 

La® + CO® 

If tan (f> — (ojcc, then we have 


—i d — ST—; 5 (« cos co< + CO sin co<) 

CO® a® + CO® ' 


' V(a® + CO®) ^ — — e““‘sin(6 — <^)] 

Thus the operational method is of assistance in evaluating 
certain integrals. For further information on this method 
reference should be made to Heaviside. *** 


♦ O. Hbavibide : Electro -magnetic Theory^ Vol. Ill, p. 234. 



CHAPTER VI 

TRANSIENT EFFECTS IN MACHINES 


General. Machine work provides many illustrations of transient 
phenomena, but complete investigations are impossible, due to 
the complex behaviour of iron in the machine. However, if 
we neglect this effect in the interest of simplification, sufficient 
data are available to deal with the problems. 

1 . Free Period of D.G. Armature : its Motional Capacity Effect. 

In the case of a d.c. motor operating on constant voltage V we 
have that V — E Ri, where E is the generated volts, B the 
resistance, and i the armature current. But we have 


n 


« = '‘^a-60 
= nip 


. 10 - 


the induction 


V 10 ~ ® 

where v = no. of poles and ip = <f)Z - . 

factor. If 6 is the displacement of the armature in geometrical 
radians, then Kddidt is the armature velocity in r.p.m. So we 
may write 

r, dd 


Further, the torque oc Ei ccbipi. 

So the voltage equation becomes 

V = Ri Kipddfdt . . . (165) 

At light load the torque developed = accelerating torque or 
Jd^djdt^ — bipi . . . (166) 

where J is the mass moment of inertia of the armature. 

As the voltage is constant, we get 

Rpi + KippW = pV =0 . . . (167) 

Substituting from equation (156), we have 

Rpi + bKipHJJ = 0 . . . . (168) 

At < = 0, i = V/R, so we obtain 

i = {V/R) 

65 


. ( 169 ) 
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Comparison of this equation with that for a condenser 
circuit, namely i — (VIB) shows that the motional 

capacity effect of the armature is given by J/bKip^. 

When the armature inductance is taken into account, the 
following equations are obtained — 

{R + Lp)i + Kippd = F . . (160) 

by)i = Jp^ . . (161) 

or i[Lp^ + Rp bKipyj] = 0 . . (162) 

When i is oscillatory, then Ry4:L^ < bRyi^JL, and 

The frequency of oscillation is 


1 l/bKy,^ R^\ 
2waJ \ JL ~ 4.Ly • 


. (164) 


Neglecting the value of R as being very small, the frequency 


. 1 


IS 


27r 



and has a nominal value of about J to Vo cycle 


per second. The critical resistance below which the current 
oscillates is given by 


R = 2y/{bKy>^L/J) 


2. Stability of a Motor Operating on Fluctuating Voltage. 

When the applied armature voltage is given by V sin Oi, where 
Q is the frequency of the armature voltage, we have 

Lpi + Ri + Ktppd — V sin 
and Jp^d = byii 

By the above method we get 

LpH -j- Rpi + bK y>H/J = VCl cos Q.t . 

As we are interested in sustained oscillations only, we will 
consider the particular integral, namely — 

cos (Q< + 0) 


(165) 

(166) 

(167) 


V[-B" + (bKrpyja - LQ)^] 
and from equation (166) we have 

_ lixpijj 

so that 

6^ cos {CU + 0) 


( 168 ) 


0» = - F 
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It will be observed that 6 has maximum values as follows — 

1. With B small, when = Kbyi^jJL. 

2. With R s v'[(2 - -s/3) bKy>^L/J], 
when Q = Kbyj^/JL\/3. 

3. When H = 0. 


Also the value of i will be a maximum when fl* = Kbtp^JL, 
i.e. when the frequency of disturbance is equal to the natural 
period of the machine, the value of the current is ( V/R) cos Clt. 

3. Starting Conditions in D.C. Machines. It is sometimes 
found in practice that a motor is slow to pick up speed on 
starting. This is especially the case with compound-wound 
machines. In these it will be shown that the fields may act 
differentially during the starting period. 

(i) Shunt Motor. For a shunt machine we will assume that 
the air-gap flux is constant and that the field distortion is 
corrected, so that the back-e.m.f. depends on speed alone, i.e. 
e = ym. 

The speed will be given by 



where a is the acceleration = torque/inertia = byii jJ. So we 
have |.| 

e = {bxp^jJ) I idt — K I idt 
Jo Jo 

where K — bxp^jJ. 

Then the voltage equation is 

Ri -f- Lpi + (KJp)i = V . . (170) 


where R and L are the armature resistance and 
inductance. Simplifying, and noting that at f = 0 — 

» = /o= V/R=c\ 
dijdt = 0 = cj 


we get* 


W + {R/L)p] 

* ” [pa + (^R/Dp -f {K/L)] • 


i.e. »■ = 4- a)£^ -I- (/5— a)£-'’‘ 

= /q e “ [cosh pt -f- (a//3 ) sinh /3<] 


leakage 


• (ivi) 

. (172) 

. (173) 
. (174) 


♦ Chapter IV, equation (97). 
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When the machine constants are such that By4:L^ < KfL, 
then (8 is imaginary, and we get 

i = /o £ “ [cos sin /8<] . . (176) 

Now, a/^ = Oil *t/ — KJL) < 1, so that the current after 
a time — ttJ 2 is negative under these conditions, and the 
motor speed will “hang up” temporarily. 

(ii) Compound Motor. For a compound -wound motor, 
the circuit may be represented in Fig. 19. Let the number of 
turns on the shunt field be and on the series field be n^. 
Due to the series field, the air-gap flux is ci^ and so the e.m.f. 



induced in the series winding will be x 10 Similarly 

the e.m.f. induced in the shunt field due to its flux acig will be 
p(aa 2 )awil 0 “*. If e is the voltage generated in the armature 
by rotation, we get the following equations for the armature 
and series field — 

El = Biii -f Lipii 4- pct\wi 10"® -|- paa 2 % 10 “® -f- e 

= Biii -f Lipii -f if (pi, + apig) -f e . . . (170) 

For the shunt field — 

= -^2*2 + + P • acit2«%10~® 4- pci,awil0“® 

= fjgia 4- L^VH + Mia^piz + api), . . . (177) 

Now, the value of the armature e.m.f. may be obtained from 
the same considerations as above. Here <f> will be given by 
c(ii -f- ai 2 ) and the torque by 6 ^ 1 ,, so that 

e = Ki<l>n = Ki<f> ^ dt . . (178) 

The complete solution may be obtained from equations 
(176), (177), and (178). But it will be difficult to manipulate. 
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80 we assume that (f> is constant or e = {Klp)i, and the equations 
may be written thus — 

El — {Ri -|- RiP + + Klp)ii + aMpi^ . . (179) 

E 2 = (i?2 + R 2 P + a^Mp)i2 + aMpii . . . (180) 


From these equations, by eliminating and noting that 
pEi = 0 and E^ — E^ = V, we have 


- (M + L)p^ + + JT -I 

,Ap^ + + C'f* + -^-^8 J 


{ 1 ) ■ 


(181) 


where A = L^M -f a^ML^ + LiL^ 

B = MR^ " 1 ~ LiR^ * 1 ” ct/^MRi R^L^ 


and 0 = R 1 R 2 + KL^ + ci^MK 


As the coefficients A, B, and C are positive, the real roots are 
negative (imaginary roots, if they exist, must be in pairs and 
conjugate). Let the roots be given by — a, — and — 6, where 
a, /3, and 5 are positive. Substituting, we obtain 



As R 2 is large, I/R 2 will be small, and the effect of the other 
three terms may be sufficient to cause the field current to be 
negative for small values of time. 

In a similar manner the armature current may be formulated 

r aMp^ 1 

as %i--E ^ Bp^ + Cp-\- AijJ • 


or 


ii = 


E 



aMp^ e~^ 
pdZJdp 


] 


(183) 


For actual values of the currents and ig reference should be 
made to Ludwig’s paper (see Bibliography). 

Various methods have been devised to reduce the interaction 
of the fields, such as reducing the time constant of the shunt 
field and adding resistance to the series field. In certain cases, 
as in welding generators, the interaction between the fields is 
accentuated. On this matter reference should be made to 
papers by Miller.* 


* Mn-LBE: Trans. A.I.E.E. (1934), Vol. 63, p. 1296; also Cbbbdy: Trans. 
A.I.E.E. (1931), Vol. 50, p. 662. 
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4. Synchronous Machines: General Equations for Rotating 
Body 

Free Oscillations. For any rotating body the equation for 
free motion is given by 

{Jig) pH bpd cB — Q . . (184) 

where 6 is the instantaneous value of angular oscillation 
measured in mechanical radians, 

J is the amount of inertia in kg. m.®, 
b is the damping force per unit angular velocity, 
c is the control force per unit angular displacement. 
This equation may be integrated as 

6 = ^e-“<sin (Q< + ^) . . (185) 


where a = damping constant = bgf^J, 
Q = natural angular velocity = 


J 



)■ 


From this equation it is seen that if « is positive, then the 
oscillations die away exponentially and the system is stable. 
But if a is negative, then the amplitude of the oscillations 
increases and the system is unstable. 

Forced Oscillations. When fluctuations of the driving 
torque are considered, as in the case of a rotor driven by a 
reciprocating engine, then the equation becomes 

{Jlg)pH + bpB + c0 = {T„, + T cos a)t) . (J) . (186) 
where is the mean torque, 

T is the amplitude of the oscillating component 
superimposed on Tav, 

CO is the angular velocity of the disturbing force. 
Generally, the oscillating torque curve would be analysed as a 


Fourier’s series and formulated as 


n a o 

1'^ „ cos {ncot + <!>„). 

n- 1 


Here we concentrate on the simple form given above, and 
the additional deflection due to the varying component of 
torque is given by 


6 


COB {cot if>) 

Vlb^co^ + (Jcoyg-c)^] + JCig 

(a® -f- £2®) cos (£2f + V*) 1 

\'\_{a.^ + £2® — to®) + 4a®£2®]j 


. ( 187 ) 
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where a and Q are as defined in the immediately preceding part, 
(a® + (a** + Q2) 


and tan y) — 


2ayi2® 


It is seen that the oscillatory motion consists of two frequencies 
given by and fl/27r, the latter being more pronounced at 
the start and dying away after a time and leaving the forced 
oscillations. When the two are equal, we get a condition of 
mechanical resonance, and the resultant amplitude may be so 
large as to cause damage. When the forced oscillation is repre- 
sented by a series expression, then there are more chances for 
resonance to occur. However, the amplitude is greatest at the 
gravest frequency, and it is usual to base design work on that 
frequency. 

6. Ssmchronous Motor. We now apply the above to the case 
of a 2v-pole synchronous motor which has a uniform reluctance 
for all positions of the armature. Let R and denote the 
armature resistance and synchronous reactance per phase. 
These we consider constant. From the theory of a.c. machines 
we have the following developments for the synchronizing and 
damping torques, which give the values of c and b in the above 
torque equations. 

(i) Synchronizing Torque. This is known to be propor- 
tional to the angle of displacement, and if the machine swings 
0/2 ahead of the mean position, then the net change of power 
per phase may be shown to be 

AP = {E^XJZ^) sin (0/2) cos {612) — El sin (0/2) sin ^ (188) 


So the torque per unit space angle displacement from the mean 
position for an w-phase, 2y-pole machine may be shown to be 

mv® AP mv^ r E^Ls El sin ^ 

= (0/2) . 47r/ ^ ^ 

The first term is due to the circulating current, and the 
second is due to the change in phase angle between E and I. 
With cos ^ = 1 and R neghgible, we have 


(189) 


mv* E^ 
goi ' X, 


(190) 


(ii) The Damping Torque. When the machine moves ahead 
of its mean position it will tend to give out more power and so 
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slow down. If it lags behind, it will tend to motor and so speed 
up. The forces set up will be proportional to speed, and so 
to pd. It may be shown that the power increment is given by 

dP = EdI co3<f> = RjZ^ . . . (191) 

where is the max. value of the flux linkages. Also we have 
dP = torque x space radians per sec — T . dco 


So that torque per radian per second 

dT 4>o^Rv^(R^-(o^L^) 
^ ~ d(o~ 2g {R^ + co2L2)2 


. (192) 


For an w-phase machine, when R is small this reduces to 

b = — E^v^Rmlgoi^L^ . . . (193) 


The negative value of b denotes that the machine is unstable, 
and the free oscillations will increase indefinitely. Insertion of 
damping windings in the pole faces ensures that b is positive. 
This introduces non-symmetry, and the problem becomes 
unwieldy.* 

In the case of forced oscillations discussed in Section 4 above, 
the values of b and c will change with change of speed and of 
displacement respectively. The synchronizing torque will 
consist of two components, one due to steady displacement and 
the other due to the oscillating displacement. To determine the 
final value of the synchronizing torque a graphical method 
due to Rosenbergt is used. Illustrations of this may be found 
in texts on electrical machine design. J 
6. Parallel Operation ol A.G. Generators : Self-excited Vibra- 
tions. A problem which is closely allied to that of synchronous 
motor operation dealt with in the last section is that of parallel 
operation of a.c. generators. Here the successful operation of 
machines in parallel depends on the characteristics of both the 
alternator and the turbine, or in other words we must take into 
account the governor stability as well as the alternator stability. 

For the governor, we adopt the following approximate force 
equation for free vibration — 

mp^x -f- Cipz -f- k^x =0 . . (194) 


• See e.g. Pebscott and Riohaedson: Journ. I.E.E. (1934), Vol. 76, 
p. 497. 

t Rosbnbbeg: E.T.Z. (1902), Vol. 23, p. 460. 

i E.g. Miles Walkbe: Specification and Design of Dynamo-Electric 
Machinery (Longmans). 
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where x is the movement of the governor sleeve, 

m is the equivalent mass of the rotating parts of the 
governor referred to the motion of the sleeve, 

Cl is the damping force constant of the governor due to 
the oil dashpot, 

k is the restoring force per unit displacement of the 
sleeve due to the spring. 

For the alternator we have, under similar conditions, that 
Jp^d + C 2 P 6 -|- = 0 . . (195) 

where 0 is the displacement in mechanical radians of the rotor, 
J is the mass moment of inertia of all the generator and 
turbine rotating parts, 

Cg is the damping constant due to the extra power torque 
developed and to the braking action of the damping 
winding, 

is the synchronizing torque of the generator. 

As we are going to combine both sets we must measure all 
quantities in the same units, say in lb. -sec. units. 

In the complete set the equations (194) and (195) may be 
combined. So we have certain relationships between the 
various quantities. For a small speed change (Aco) let the 
movement of the governor sleeve be x. This movement will 
open the steam valve, and so the engine torque will be increased 
proportionately to x by AT, say. So now we have AT — kx, 
where k is a constant. But due to sleeve movement, we imme- 
diately have a restoring force on the governor amounting to 
— cAo) — — cpd, so that our equations for disturbed motion 
become 

mp^x -j- Cipx + kiX — — cAoi — — cpB . (196) 

and Jp'^0 -f- -f- k^Q = kx . . . . (197) 

and we get 

[(Jp^ + C 2 P + k^) (mp^ -f Cip -f ki) + kcp] 0 = 0 . (198) 
From this equation we have the determinantal equation 

p* -f ttjp^ -f ^dsp -i-at = 0 . . (199) 

where 

— Ci/m -j- cJJ ; — kjm m 

— (^ 1^2 + <^ 2^1 + ck)IJm; = kikJJm 

6— (T.152) 
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Now we may introduce constants derived from the con- 
sideration of each unit working independently as follows — 

(a) Damping coefficients = — cJ2m; aj = — c^2J. 

(b) Natural angular velocities cdj^ = kjm ; Wg® = k^/J. 

It is seen that all coefficients in our determinantal equation 
are positive; further, we obtain by reference to Section 9, 
Chapter I, equation (16) as the criterion of stability, that 
must be positive, or, in the limit, 

aiagOg = ag® -f . . . (201 ) 

We may write down our coefficients as — 

«! = — 2(ai -f ag) ; ag = "i® + 

ttg = — 2(x.2(Oi^ — 2aift>2® + 2y ; a^ = coi^cog® 

where y = ckl2Jm. 

Now, suppose X is the maximum sleeve travel and Aco is the 
maximum change of speed, then we have 

cAco = kiX 

or c — ki{xlAa)) 

Denoting Aco/w by the governor regulation, and noting that 
the normal torque, T^, is given when the sleeve movement is x, 
then we have 

kx = T„ 

Q _ fc kiX _ T„ k^ 

^ ^ ^ 2Jm ‘ a)(Aft)/<w) 2d ’ J(om 

= (ji}^T„l2Joid = (Wi®/2M 

where t = Ja)IT„ is the time to bring the set up to full speed 
from standstill. 

Substituting these relationships in equation (201), we get 

(«i + «2)[*>2*(1 + P*) + ^ajag] [ag -f aj/p® — l/2^<] 

— <Og®(ai +• «a)® — p^(02\«.2 + ai/p® — lj2dt]^ — 0 (202) 

where p = Wi/wg. 

From this equation it may be deduced that — 

(а) When the governor natural frequency exceeds the 
generator natural frequency, then damping of the governor is 
necessary. 

(б) When the governor natural frequency is less than that 
of the alternator, then damping windings are necessary. 

(c) For isosynchronism, damping is necessary in both 
governor and alternator. 
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7. Pulling into Step of a Synchronous Motor. Here we have 
an interesting but insoluble problem. In starting a synchronous 
motor the usual procedure is to supply the stator winding with 
polyphase currents; the short-circuited damping winding on 
field system acts as the rotor winding of an induction motor. 
The motor thus runs up to a subsynchronous speed. In this 
condition the field poles will slip backwards through the 
synchronously rotating stator field. When the field system is 
energized, the two fields will eventually lock together and the 
rotor will then move forward at synchronous speed. 

The phenomenon which occurs between the induction-motor 
and synchronous-motor stages is exceedingly complex, and the 
following points are of importance — 

(а) Due to the damping winding being asymmetrical, the 
induction-motor torque will be uneven. 

( б ) With a sahent-pole machine the reluctance of the air gap 
varies considerably, and will aggravate the non-uniformity of 
the torque under (a). 

(c) When the main field switch is closed, then the field 
current will build up according to the law 1(1 — The 

rate will depend on the time constant of the field, and a period 
of time will elapse before the full synchronizing torque is 
experienced by a motor. 

With these points in mind, we can enumerate the torques 
that will be experienced and must be considered for any 
puUing-into-step operation. 

(а) The amount of the shaft load torque: the greater the 
load, the longer the time that will be required to accelerate it 
and the greater will be the slip of the induction motor. 

( б ) The torque developed as an induction motor and the 
characteristic (torque/slip) curve : under this heading the effect 
of the unbalanced windings and reluctance must be considered. 

(c) The torque developed by the motor when the field poles 
are excited: this will depend on the angle (9) between the 
excitation and armature magnetic axes at the moment of the 
closure of the switch. 

Thus we get an equation of the type 

PiP^d -k -^ 2(1 — b cos 26)p6 -\- Pg sin 26 

-f P 4 sin 0.(1) = Pg (203) 

where P]P ®0 is the inertia power and is equal to 18-6( Wr®) (P/v^) 
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X 10 kW. per electrical degree per sec®, x {d^djdt^) and 
where Wr^ = inertia in lb. ft.®, 

/ = frequency, 

V = number of poles. 


The induction-motor power = P^fSGOfs kW./electrical 
degree/sec., and (1 — 6 cos 26 ) is a term to take into account 
the asymmetrical damper windings. Pg sin 29 is the reluctance 
torque. The value of Pg is given by Doherty and Nickle* as 


S{xg—x,) F® 


2xaX^ ' 1 000 


kW., where Xg and x, are the reactances in the 


direct and quadrature axes respectively. P4 sin 0 . (i) is the 
synchronizing effect of the field current. The unit function sign 
is used to denote that before the field is applied this term is 
zero. 

The solution of these equations cannot be derived by ordinary 
mathematical means, and recourse is had to an integraph.t 
This method has been applied to the solution of various 
problems, t 

From the results obtained in the two papers quoted below 
for the limiting values of slip for successful pulling into step 
we note — 

(a) For a cylindrical rotor 


/r ^4 1 

n VL/(WV®)J 


(204) 


(6) For a salient-pole motor 


8 = 


^ /r 1 

n VLf(Wr^)j 


(205) 


8. Damping Action due to a Shaded Core. In contactor and 
relay work, the iron core is shaded by a close-fitting copper 
ring. In the case of synchronous machines there is the allied 
problem of the short-circuited damping windings. Due to the 
analogy existing, we will here consider the first case, and to 
simplify our problem we assume a laminated core so that eddy 
currents are negligible. 


* Doherty and Nickle: Trans. A.I.E.E. (1926), Vol. 46, p. 912. 
t Bush, Gage, and Stbwabt: “A Continuous Integraph,” J. Franklin 
Inst. (1927), Vol. 203, p. 63. 

t E.g. Edoerton and Fourmaribr: “Pulling into Step of Salient-pole 
Machines,” Trans. A.I.E.E. (1931), Vol. 60, p. 769; Edoerton and Zak: 
“Pulling into Step of a Synchronous Induction Motor.” Journ. I.E.E. (1930), 
Vol. 68, p. 1206 . 



TRANSIENT EFFECTS IN MACHINES 


77 


Consider an exciting winding of n turns on a core arranged 
as in Fig. 20. Then 


V — Ri n d(})ldt . (206) 



1 

with (j) expressed in 10^ megalines. 

L 

epi\/ 



The current set up in the solid copper 
ring will be 

{\IR,)mdt. . (207) 


•trlr 



T 


1 


where R, = resistance of the copper ring. 

Now, the m.m.f. in the core will be 
the vector sum of that due to the two windings : we may write 


(f) — K{ni — i,) 

where K is the flux per ampere-turn 

and Kni — <f> + Ki, — [l -f- {KIRg)p](f> 

So that i = (IJKn) [1 -f- {KIR,)p]<f> 

Substituting in (206), we get 

V — {RIKn)<f> -f {RIR,n)p^ -f- np(f> . (208) 

Denoting KnVjR by the steady-state value of the flux, we 
obtain 

= [1 -f {KIR,)p + {Kn^lR)p]<f> 

Hence we obtain by expansion that 

^ (1 - exp. [- tliKIR, + Kn^lR)]} (209) 

When the copper ring is absent we have 

^ (1 — exp. [— ll{Kn^lR)]} . . (210) 

It is seen thus that the presence of the short-circuit ring 
causes the time lag of the flux to be increased. The increase 
depends on the resistance of the ring and on the square of the 
number of turns in the exciting coil. In any machine equipped 
with short-circuited windings the flux will be tardy in either 
the building-up or the decaying processes. 

9. Short-circuit Current of an Alternator : Symmetrical Fault. 
The calculation of the transient short-circuit current in a.c. 
machines has been the subject of several papers, each author 
eliminating one or more of the assumptions of his predecessor. 
Naturally, as such work proceeds the mathematics become more 
cumbersome. We may start with the elementary discussion 
given in Section 4, Chapter V, wherein we assumed, among 
other things, no variation in the fleld current during the short 
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circuit. We are able to remove this error by a method developed 
by Ku.* Here we assume that the stator windings are sym- 
metrical, the rotor speed constant, and the resistances and 
inductances of the stator winding constant and independent of 
rotor position and of current; also that the mutual effect 
between rotor and stator depends on the relative displacement 
angle; and that the machine is not saturated and the har- 
monics in the distribution of the air-gap flux are negligible. 
This is very close to the case of the machine with a cylindrical 
rotor. 

When a three-phase short circuit occurs on a machine, the 
phase voltages v^, v^, become zero. So we may consider 
that they are annulled by equal and 
opposite voltages at short circuit. The 
current may thus be assumed to be the 
steady-state current at the moment of 
short circuit plus the transient current 
set up by application of these annulhng 
voltages at the terminals of the machine. 
Let the currents, etc., be shown as in 
Fig. 21. Then we have the following equations for suddenly 
applied terminal e.m.f.s 

(f?! -f- Lip)ii + Mpif cos 6 = Vj . (211) 

(i?! -j- LiP)iz + Mpif cos (0 — 277/3) = . (212) 

(i?i -j- + Mpif cos (6 -f- 277/3) = . (213) 

(i ?2 + L^)if + Mp{iy cos B -f- cos (0 — 277/3) 

-f *3 cos (0 -f 277/3)] = 0 . (214) 

Here = synchronous self-inductance of a stator phase plus 
the mutual effect of the other stator phases. 

As we are considering a symmetrical three-phase short 
circuit without neutral connection, we may write 

*1 “f" *2 "I” *3 “ ^ 



+ ^2 + ^3 = 0 

So equation (214) becomes 

(jRa + L^)if 4- cos (0—77/6) -f cos (0 — 77 / 2 )] 

= 0 (216) 

Further, we know that 

= / cos oit = (//2)[e^“‘ e — say. 


* Ku: Trans. A.I.E.E. (1029), Vol. 48, p. 707. 
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and *2 = I cos {cat — 27t/3) = 

= ie-i2W3 ie + i2W8 

= + ta, say. 

On substituting these in equation (211) and equating like 
components of current and voltage, we have 

(i?! + + \M'pifE^ = V . . (216) 

and (Ri + Lip)i + iifpt>-^® = i^ . . (217) 

Equation (216) now becomes 

(■f^2 “I” ^2P)y "I” Y^p^^~ ^ -\- ts ) = 0 (218) 


If we assume that the short circuit occurs when the excitation 
axis coincides with that of phase I, then 6 = cat. On multiplying 
equation (216) by and (217) by e^“‘ and using the 

shifting theorem, we get 

[Bi + +jco)]ie-^ 4- ^M[p + (219) 

[B, + L^lp — j(a)]le^“‘^ -}■ h^lP ~ • (220) 

From equations (218), (219), and (220) the values of » and t are 
obtained, as follows — 


Z(p) 




3 

4 LiL^ 


3 

pip —j(o)v + j pip + 3 (a)v 


and i is given by a similar expression with v substituted for v 
and jca for — jca, i.e. it is a conjugate expression. 

Here o = (1 — and 


Zip) = i»» + 


ri Bi / Bi . 
LotI^AA 


L, 

2i?2 


+ 

\ 1 1B./B,^ 




) 


From Zip) the roots are obtained in the usual way. With the 
normal machine constants, one root is real and the other two 
complex. 

When the voltage v is given by V sin cat = i F/2j) — e~ 

Then v = (F/2j)e’“‘ and v = — (F/2^')fi~^"‘ 
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The values of the currents are now given by 


2j ^ aL 


j(ot . 




V 

2j 


1= 


[( 

3 

ip + ^ilLt — joj) — 2 

rW) [(^ + 1) 

(y + §;+>)- 1 


As these currents are conjugates, we have 
i — i -{-i — 2t (real) 




LiLg 



When the short circuit occurs under load conditions there is an 
angle a between the excitation e.m.f. and the terminal volts. 
On taking account of the time (0) at which the switch is closed, 
the terminal voltage will be fully specified by 
As further analysis is difficult, we will now use the results 
given in Ku’s paper to show the effect of short-circuiting the 
machine at no load but with 0 = — 71 '3°. 

The machine constants are as follows — 

Ri = 0-0232 i ?2 = 0-138 M = 8-34 x lO-^ 

= 5-66 X 10-3 ^ ig.3 X 10-3 ft, = 377 

a 0-046 RJLy = 4-1 RJL^ ^ 7-13 

1/(tAi = 3 840 Vr,^a^ - 188 

Z{p) =--p^ + 248p2 + 144 X 103^ -f 220 X 10® 


The roots are Pi — — 169 ; ~ ~ i ^359 

With p^ — — 169, the contribution is found to be 


"Pi 


1 88 

.•1_ pi(377«~71°-90°) 

2 


X 


3 840 X 6 100 89-5° 

ififlt 

— 246 X 10« 


_ 394 g-109« gjISTTt+lon 


Twice the real part gives the total contribution as 
1 790e-i«3‘ cos (377« + 109°) 

The value of the total current is given as 
i = 1 790£-3«3* cos (377< + 109°) 949£-39-6< cos (736<— 51°) 

-f- 1 038£-39-6« cos (18<- 88-7°) -f 88 cos (Silt + 109°) 
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The last term is the sustained value, and it will be noted that 
the largest term is attenuated the most. The term of double 
frequency is not so strongly damped, and there is a further 
term of about zero frequency. 

The expression for the field current under no-load conditions 
may be shown to be 

3 MV o)p 
^ 2 ' "^) 


In this case it is given by 


if 


= 94 X 10‘ X 188 X 


P 

Z(p) 


With the value of the roots as given we have that the con- 
tribution of Pi is 


. . 

= 94 X 10^ X 188 X - ^.6 ^ X 


== 1 


On evaluating the contributions of the other roots, we have 
if = 1 2106-^^^^ + 1 cos (359i- 160°) 

i.e. a current of practically fundamental frequency super- 
imposed on an ex])onential curve. 

Of recent date several contributions towards the solution 
of unbalanced faults have been made. As these are difficult to 
summarize, we aj^pend some references dealing with this 
subject. 
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CHAPTER VII 

PROPAGATION OF ELECTRIC WAVES ALONG A LINE 


1. Steady-state Solution. A transmission line may be regarded 
as composed of sections ; the greater the number of sections, 
the closer is the representation. In the limit, with sections of 
length dx, the representation is complete, and solutions ob- 
tained under these conditions are designated as smooth line 
solutions. The arrangement of a section is shown in Fig. 22, 



where R is the resistance, L the inductance, C the capacity, 
and O the insulation conductance, all per unit length of the 
line. In an elementary way, the flow of current may be con- 
sidered to occur as follows — 

Due to the impressed voltage at Aa, a current will flow 
in the section ABba. The condenser element will assume a 



charge, and a voltage will be impressed on the next section, 
with a resultant current therein. Repetition of this process 
eventually gives a voltage at the far end of the line. The 
initial diffusion of electricity through the wire is termed the 
building-up process, and the final stage is the steady state. As 
the theory concerning steady-state conditions is well known 
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to all power engineers, we will give a brief summary and 
develop therefrom the theory concerning the process of building 
up. 

For the section given in Fig. 23, let the voltages and currents 
be as shown. Then 

de 

e -j- = e -j- i(Jidx) -j- (Ldx)pi . . . (221) 

di 

and i ^dx = i -\- e{Odx) + {Cdx)pe . . . (222) 

where x is measured from the receiving end of the line. 

These equations reduce to 

~ = ii? + Lpi 
di 

and ^ ~ 


for harmonic quantities the potency of p — jw. So we have 



^ = i[ii + joiL] = iZ . 

. (223) 

and 

~=e[0+j<^ = er . . . 

. (224) 

By differentiation we have 



di 

= ZYe = oL^e 

dX^ dX 

. (226) 

and 

dH de 

- ZYi = 

dX^ dX 

. (226) 

where 

a - ^/{ZY) = V[(^ {G -f joC)] 

(227) 

A solution of these equations is 



E^ = 

. (228) 

and 

h= {cLlZ)[Ae^^- Be-^] 

. (229) 

where A 

and B are constants, determinate from 

terminal 

conditions. 


For example, at a; = 0 let Eg, 1„ be the receiver 
current respectively. On substitution we get 

volts and 


A = and B ^ l[E,- {Zju)!,] 
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So the voltage and current are given by 

E. = \[E, + {ZI^)hY- + \[E^- (Z/«)/Je— . (230) 
= + e-“]/2 + {Zlx)Jg[s^^— e-“*]/2 

= Eg cosh aa; + Ig^/{ZIY) . sinh ax . . (231 ) 

also h = W{^IE){[Eg + hViZm] . 

W(Z/7)].e— } . . . (232) 

= Ig cosh (XX + Eg\/{YIZ) . sinh aa; . . (233) 

If the voltage and current at the transmitter {Ej,Ij.) had been 

used in the determination of the constants, x measured from 
the receiving end, we have 

E^ — E j. cosh ci{l— x) — I rV (ZJY) .sinh ci{l— x)' . (234) 
Ix = iT^osh ol{1—x) — Ej.\/{YJZ) .sinh a.{l — x) . (235) 


When X is measured from the transmitting end we get 

Ex = Ej, cosh aa; — Ij^/{ZIY) . sinh oca; . (236) 

Ix = If cosh aa; — Ef^/{ YfZ) . sinh aa; . (237) 


When the far end of the line is connected to a load, let 
Egjlg = cr, the load impedance. Then we have, for equations 
(230) and (233) 


and 


Ex = Eg 

~ I B 


l^cosh aa; + 
l^cosh aa; + 


y(i) 


sinh XX 


sinh XX 


] 


Denoting a\/{YIZ) by tanh 6, these equations become 


Ex 


, sinh (aa; + 0) 
“ sinh 6 


(238) 


and 


y cosh (xx + 0) 

S3sh0 


(239) 


So then EJI^ = o tanh (al + 6) . coth 6 
= \/{ZJY) . tanh (xl + 0) 
= Zq tanh {xl + 0) 


( 240 ) 
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With these results, the following relationships for various 
terminal conditions may be obtained — 


Line Termination 

<J 

E, 

I, 


Earthed 

0 

^ sinh otx 
^ sinh cd 

J. cosh ctx 
^cosh cd 

\/{ YjZ) . tanh clx 

Open . 

00 

^ cosh ccx 
^ cosh ocl 

sinh ccx 
^sinh cd 

y/( YjZ) . OOth OLX 

Infinite Line . 

^/(Z|Y) 

1 

T 


V{ZIY) = 


2. Roots of a. The value of a is given by 
V[(i2 + i^) ((? + Cp)] 
or for harmonic quantities with p == jcj — 
a = 

= /? ± jy, say 

The values of /3 and y are obtained by squaring both values 
of a and equating the reals and imaginaries. The following 
results are obtained — 

Vim-lYl + BQ-co^LC)]. . (241) 

y=V[W\-\Y\-RO + o>^LC)]. . (242) 

where \Z\ = and \Y\ = 

By differentiating the expression for /3, with respect to co, 
it may be shown that /? is independent of frequency when 
R/L = G/C and a has a value of \/(R0). Under like conditions 
y = co-y/ [LC). 

The physical meaning of these expressions may be illustrated 
by considering an infinite line. With x measured from the 
transmitter, we have 

Now, = cos yx — j sin yx. So we have 

Ex — Ej£~^^[coa yx — j sin ya;] 

Thus an applied voltage Ej, is attenuated to due to 

passage along a line of length x. In addition, it suffers a phase 
change of ya;. Thus, when ya: = 277- or a; = ^rrfy miles, the cur- 
rent and voltage at x are in phase with the sending current 
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and volts ; but when xy = n ox x ~ Trfy, these quantities are 
in phase opposition. A spiral diagram shows this idea clearly. 



Fig. 24 is drawn for a line such that = 0*0024, and y = 0*0176 
radian per mile. The voltages are tabulated in Table I. 

TABLE I 


X 


y 

X 


y 

0 

10 

0 ^ 

260 

6-64 

262 ° 

10 

9‘7 

10 ® 

300 

4*88 

303 ® 

20 

9-38 

20 ® 

360 

4*43 

363 ® 

40 

912 

40 ® 

400 

3-86 

404 ® 

80 

8-27 

80 - 6 ® 

600 

303 

606 ° 

100 

7-87 

101 ® 

600 

2*38 

606 ® 

160 

6-96 

161 ® 

700 

1-88 

707 ® 

200 

6-23 

202 ® 





It is here noticed that at ^Tijy = 27r/0*0176 = 35.5 miles the 
voltage is in phase with the sending voltage. 

It is customary to define p as the attenuation factor and y 
as the wavelength factor. 

Effect of Frequency on p and y. The value of p may be 
written as 

2p^ = ^[{R^ + Co2A2) ((J2 ^ ^2C<2)] ^EO- co^LC 

On differentiating this expression with respect to <w, we obtain 

d(2p^) fr Yl® 

do) ~ ^\r\R‘^ + J 

The value of p increases with increase of frequency, and has 
its minimum value when the frequency is zero. If OR = LO, 
then the value of p is given by \/{RO), which is independent 
of frequency. Likewise it may be shown that for the wave- 
length factor y when CR — LO, then y = o)^/ (LC). The 
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velocity of propogation is now l/\/(LC'), which is also inde- 
pendent of frequency. 

Thus for a line where CR = LG there is no distortion of 
attenuation or of velocity due to frequency change. However, 
phase distortion may still exist. As actual line constants do 
not meet this specification of CR — LG, inductive coils known 
as loading coils are placed in the line. The usual type of fine is 
composed of two parallel copper wires. The capacity and 
inductance per unit length are readily calculable. For lines 
of diameter b and separation d, with d > 6, we have in air that 


Ign. {2dfb) 


and 



Usually R and G are negligibly small, so 
a — i iy — j(*>{LC)^ 

y — and co/y = = 3 X 10® m/sec. 


Also the surge impedance is 



= 276log,o. 



ohm 


For high-frequency work, coaxial cables are now used. They 
consist of concentric conductors separated by porcelain spacers. 
The large area reduces high-frequency resistance loss, and as 
the electric and magnetic fields are confined to the region 
between conductors, interaction is reduced. Here 


C - 


277Kn 


log. i^iK) 


and L 




The velocity of propagation is as before, but now the surge 
impedance reduces to 138 log^Q (r^/rg). 

3. Lines in Series. We may rewrite the line equations as 

Ej, — EgA -(- IjjB 
Ij, = /gA EfG 


where A = cosh cd, £ — y' {ZJY) . sinh aZ, and G = \/{YIZ) 
sinh aZ. 

For two lines of different constants connected in series, where 
the first has a surge impedance Z^ and the second a surge 
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impedance Z^, let the voltage and current at the receiving end 
be Eg and respectively. Then the voltage and current at 
the sending end of that cable will be Ej^ and Ij^ respectively, 
or 

Ej.2 = EgAi 4 " 

and = IsA + EsCi 

The voltage and current at the sending end of the first 
section are 

-®ri ~ “h -^r2-®2 

and In = 

Thus by simple substitution we have 

Eti ~ Ej^AiA^ 4" 4" IAB^A^ 4~ -^yB^ 

and In — I^A^A^ 4- B^C^ 4* E^A^G^ 4“ -^yG^^ 

If now the impedance of the second line is equal to its surge 
impedance, i.e. the line is of infinite length, then 

Et 2 ,I^T 2 — Z-j^= l/I^i 

The impedance of the whole system is 

7 Eiy _ -^2 ~h B^Yi 

^ - In Y,A, 4- <?2 

cosh 0.1 4- Z 2 Fj sinh olI 
Yi cosh al 4- {^ 1 ^ 2 ) sinh cd 
y [Zi cosh aZ 4- Z 2 sinh aZ 
. . j — 2 cosh od 4- Ziy siqh cd 

When Z 2 is almost equal to Z^, say Zi = Z^ Z^, it is seen 
that 

r cosh cd 4- sinh c/1 4- {Z^jZ^) cosh aZl 
~ ® Lcosh cd 4- sinh cd 4- {Z^Z^) sinh azj 

= 7 > (^a/^2) (cosh cd — sinh cd) 

~ ^ L ^ cosh cd 4- sinh cd 4- (Z^Z^) sinh od 

Remembering that ZJZ 2 is small, we have 

7-^7 [y I ^ cosh cd ~ sinh aZ I 
^ 1_ ~^Z 2 cosh cd 4- sinh aZj 
= Z 2 4- Z^e-^ 


7--(T.I52) 



90 


TEANSIENTS IN ELECTBIC CIECUITS 


4. General Solution in Terms of Characteristic Vibrations. In 

order to determine the value of the voltage and current at 
any time and place along a hne, we may use either the expan- 
sion theorem and obtain a result in characteristic vibrations 
or the equivalency of operators and obtain a result in waves. 
The first method, dealt with here, yields results which usually 
converge slowly, so are cumbersome to handle in calculations. 
With the following notation 

= IjLG, a = RJiL, b = QI2G, p = a -f- 6, and a — a — b 

the general transmission equations, with x measured from the 
transmitter, may be written as 


e* = -f 

• (243) 


. (244) 

a = {iMy'iip -f p)2 — a2] , 

. (246) 


Open-circuited Line. In this case we have 
ej: = Eq a.t X = 0 
i = 0 at a: = 1 


as the boundary conditions. The voltage and current equations 
for an applied steady voltage become 

cosh a(l — x) 


ex 


E, 


cosh al 


(I) 


(246) 


E(,/p-\-2b\i sinha(Z— a:) 

~ Lv\p 2a ) cosh ctl 

Considering the voltage equation, we have 

Z{p) = cosh od 
With Z{p) = 0 we get that 

0.1 = ±jn(7T/2)* 

But aV = (p -f p)2 — a® 

•''P = - j = - P±3^n, say 


(247) 


(248) 


♦ n is odd. 
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Now 


dZ dZ do. , . , , P + p 

_ — . — — I sinh at . y- 

dp dx dp av^ 


, Pf sin (W7r/2) 1 
P sin (rt7r/2) 


= ± j^n ■ 


W 7 r /2 


With p = 0 we obtain 


YjO) 

Z{0) 


cosh {[2(Z — a:)'v/ (a6)]/v} 


cosh {[2Zy'(a6)]/v} 

Thus the complete expression for the expansion is 


E, 


cosh j . {l-x) 



n-n 


cos (I — x) . 


± j^n(— P ± j^n) 


n=^l 


P sin (w7r/2) 
V® n‘rTj2 


( 249 ) 


( 260 ) 


( 261 ) 


The summation sign may be simphfied, since 
g-py £}P„t X _ — 2s~'*[p sin p„t + /?« cos fij) 

Wn\-P+j^n~-p-j^n)~ WTJJ) 

— 2 cos {ftj — <f>n) _p, 

^nViP^ + ^n^) ® 

where tan = p/jff„. 

The expression for the voltage becomes 

( cosh [2V(a6) {l-x)/v] 

* ®{ cosh [2(y'a6)Z/v] 


-.t V » sin {n7Tx/2l) cos 
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At X — I we obtain the receiver volts as 

n=QO 


E„ = E, 


'4 


TTV" 


cosh [2{'\/ab)l/v] 


n“l 


The value of the current may be obtained from 
— di/dx = ((? + Cp)e 

This gives 


- ( Ov sinh [2(\/a6) (Z— a;)/t’] , 2v^ , 

^ ^«(2^ cosh [2{^ab)yv] T ^ ' 


n=«l 


cos {niTxJ2l) p 0 cos — <f>x) 

jn \rv(^T¥jr' 


C sin 


in/Snzjj 


(264) 


With Q -- L = 0 the cable is the Heaviside “ideal.” On 
making these substitutions, i.e. a = \/ (EGp), we get 


n=-l 


2 ^ 2 
n=l 


sm . g-n^'tucRp (266) 


2^ 

Rl 


0 |■g-1T•^/4CiW• _j_ g-MliCRf ] 


For very small values of t we get 
i = {2EJBI) 


Characteristic Frequency of Line. The second term of 
equations (252) and (254) contains a component having an 
angular velocity The frequency will be 



1 // 1;% V 

^ V\~~^ 



. (267) 


Consider the case when a = 0, i.e. R/L = O/C. Then we have 
/„ = vnl2l, /i = vl2l = Il2ly/{LC), = dl2ly/{LC), etc. 
This gives the frequency characteristic of the line. 
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Earthed Line. We give one more example of this method, 
that of an earthed line. The boundary conditions are 

E ~ E, at z = 0 
E — 0 at X = I 


So the equations are 






sinh cd 

p 2b \ cosh a(Z — x) 

p 2a) sinh cd 




(268) 

(269) 


As the development is similar to that of Section 3, we 
merely sketch the method. 

Z{p) = 0 gives cd — jutt (w — 1, 2, 3, . . . ) 

— p± ^((72 _ = — p± jpn 

P cos (w7r/2) 


P cosh cd 

Z'{p) = -,ip + p)-^==±jPn 

F(0) _ sinh [2{\/ab) (I — x)lv] 
Z(0) sinh [2{\/ab)l/v] 

Expansion gives 

[2(^ab){l-x)/v] 2 v^7t 


»l7r/2 


(260) 


e, = eS‘^ 

_ SI 


sinh [2{\/ab)llv] 

2 n . ( mTx\co^ {8„t — (hn)'] 

) Vip^ + r 


(261) 


Ir 


„ ((? cosh [2(A/a6) (/— a:)/w] 2v^ 

»inh[2(V<.6)W T‘ 


-Pt 


in/9„<Jj 


2 COS {mTXfl)l'0 GOS {^„t—(f>n) ^ . 

Tn IVi^^TJT) 

The characteristic frequencies now are given by 
/„ = nv7T/27Tl for R/L = OfC 

= ll2WiLC), U - etc. 


(262) 


(263) 


5. Difference Equations: Recurrent or Inactive Networks. 
Certain networks used in practice are built up, and so their 
characteristics are lumped or non-uniform. Such networks 
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must be considered as discrete systems (a smooth system is the 
limit of a discrete system). To handle such systems we resort 
to difference eqmtions. Examples are now given of this method. 

Artificial Line. To simulate actual line conditions an 
artificial fine is built up, so permitting laboratory investigation. 
We will consider the line to be built in T-seetion. For the 
pillar we will use l/(G*+j<w(7 ) = B and on the top we will use equal 
impedances {R + ju>L) = 2A. The general arrangement is as 
shown in Fig. 25. For such a network to simulate the smooth 



line it is essential that the number of sections should be large. 
But when terminal conditions only are required, then it is 
possible to accomplish this in one section, provided that the 
following values are adopted — 

B = Zo/sinh aZ 

where and a are determined from smooth line characteristics. 

Let the network have n sections. From Kirchhoff’s laws, for 
the r** section we have 

if2A + ifiB— i^_iB— = 0 . . (264) 

or ir[2A + 2B] — B[i,_i + = 0 . . (265) 

Let ir = Pe*’’ + Qe-" . . (266) 

Then by substitution we get 

[(2A + 2P) - P(e« + £-«)] (Pe“’- + Qe-^^) = 0 (267) 

2(A + B)IB = e“ + £ “ = 2 cosh a . (268) 

i.e. 1 4- AJB = cosh a is the value of a to satisfy equation 
(266). 

For the determination of P and Q we proceed as follows. In 
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the initial section the voltage is and r = 0 ; for the nth 
section the voltage is and r — n. Then 

E^ = P{A + B— £e“) + Q{A + B- 5e-“) . . (269) 

E^ = P{A + J5 — + QiA + B—BE^)e-”^ (270) 


Whence 


E^ + 

2B sinh a sinh net 


and Q = 


E^ + 

2B sinh a sinh na 


(271) 

(272) 


r E, + E,s-« 1 r E, + E,s-’-‘ 1 

\2B sinh a sinh naj ~ \_2B sinh a sinh naj 
E^ cosh ret + E^ cosh (n — r)et 
B sinh net sinh a 


— ra 

(273) 


When n is very large, this equation reduces to 
j. Ej^ cosh (n — r)et 
“ ~ B sinh net sinh et ’ ' ' 


(274) 


Application op an Impulse. We now consider the appli- 
cation of an impulse voltage E at the transmitter. Then 


I = 


E cosh (n — r)a 
B sinh na sinh a 


. (Z) 


(276) 


In subsequent work the leakage conductance will be neglected 
so that 

2A = Lp -{■ R'> B = yCp 

cosh a = 1 + Gp{Lp -f- i2)/2 . . (276) 

Using the method of Section 3, Chapter VII, we proceed as 
follows — 

Z(p) = sinh na sinh a . . (277) 

For Z{p) — 0 we obtain a = jsTtJn, where s = 1, 2, 3 . . . 
Substituting in quotation (276), we have 
cos {snfn) = 1 -|- Cp(Lp -f R)/2 .... (278) 

R /rJ22(72- SLCil - cos Sn/n)"! 

■•^’ = -^±VL 2^ 

= — p± j^, 

where p = RI2L and /S, = \/[4i;2 sin® {S7rl2n) — /o®]. 
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Now, Z'{p) ~ \n cosh no. sinh a + cosh a sinh na.'\{dxldp) 

= n cosh na sinh a . {daldp) . . . (280) 

Also from (276) we have 

sinh a . {dot.Jdp) = LCp + RC/2 
Z'(p) = n{LCp + RCj2) cosh no. 

= nLC{p + p) cos 57r — (— 1 ynLC(p + p) • (281) 

When p approaches 0, i.e. when a approaches 0 — 
sinh a sinh na. -> na^ 


Now, cosh a = 1 + • • • = 1 + \Gp(Lp + R) 

a? = Gp{Lp + R) . . 

F(0) Cp 1 


So we have 
i 


Z{0) nCp{Lp -y R) nR 


1 ^ '^C'pn cosh (m— r)a . 

^ ~ \_nR / fULCpJp.. + p) cos sttJ 

r 1 6“'’* V ' cosh n — ra . 

“ InR'^nLZ^ (Pn + /o) (- 1)* J 


For the nth section this reduces to 


1 ^ e-W] 

8 = 1 

ri -I 


(282) 

(283) 


(284) 


(286) 


As n is very large, we may write sin {sirl^n) = snl2n. 

Then . . . (286) 

This equation is of the same form as that obtained for a smooth 
line. (See equation (262),) 

6. Filter Circuit. In communication work it is necessary at 
times to accept only one frequency band and to reject all 
others. The network which accomplishes this purpose is known 
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as a filter circuit. The ideal filter circuit would allow all signals 
in a desired frequency band to pass without attenuation and 
would attenuate all others to zero. 

Several circuits may be used as filter circuits, but as their 
actions are, within limits, the same, we select one of the form 
shown in Fig. 26. It will be seen that it is of the same form as 


f □ c 

3 B 


C I 


— I 


Fig. 2G 


the artificial line except that here the circuit is properly 
terminated. 

Let there be n sections in the network. The investigation 
proceeds along lines the same as those in Section 5, and we 
obtain for the 7 th section that 

i^ =- + Qe-^^ . . (287) 

where 2 cosh a — {AJB + 2). 

From terminal conditions we have 

E = ^i^A + B{iQ — i^) for the first section . (288a) 
and 0 == \inA + — i^) for the last section (288b) 

By substitution of (287), we have that 
E cosh (n — r)a 

~ B ' sinh noL sinh a ’ ’ ’ ’ (289) 


Actually we are only interested in the Tith section, where 


E 1 

B ’ sinh noL sinh a 


■{!) 


(290) 


Low-pass Circuit. For a low-pass filter circuit we have 
A = R A- pL and B = IjCp, so that 

2 cosh a == (1? + pL)Cp + 2 

For a steady voltage E , (i), applied at the sending end, the 
current at the receiving end is given by 


s—n 



«=1 


(291) 
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If now an a.o. voltage such as . (I) is applied, we then 
have 

ej).« 




Gp, 

(p/ + co^) Z'(p), 


;] 


The second term becomes 


Ee-'‘^ 

nL 


^±m- 




P ± jP,) + «>^] cos STT 


(292) 


or 


2Eb-^^ V [(p^ + coHp^-^,^)] sin 2pi3, cos p,t 

nL Z, /5.[(p2 + ^.2)2+2o>V-/5.^) + o*]cos«7r 


9 — 1 


It is seen that the last term of the expression gives the 
transient components, the number of which is identical with 
the number of meshes. The attenuation is given as p = E/2L 
and the responsive frequency as 

A = 

— (1/27t)\/[4i;* sin® {svl2n) — p®] 


or 


Neglecting p®, we have 

ft = (v/'”') sin {sTTf2n) 

sin (siTl2n) 

7r\/{LC) 

At 8 — n — 

/„ = l/nViLC) 


. (294) 


the upper frequency 


When 8 = 1 — 


j, sin (7cf2n) 

nViLC) 


= the lower frequency 


i.e. if we increase n we lower the frequency A, which will pass 
through. The value of/„, not being alfected by the number of 
sections, will remain the same. 

Other filter circuits may be considered along these lines. 

7. Potential Distribution along an Insulator String. In a 
suspension-type insulator let and Z^ represent the imped- 
ances between tie pins and between pin and earth respectively. 
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Then we have for n disks the equivalent circuit of Fig. 27. 
As in Section 6, we obtain 

{Z 1 /Z 2 ) + 2 = 2 cosh a 

e, = Pe^ + Qe-^ . . . (296) 

er/^2 = [(fir-l - er) - (Cr “ 

For r = 0 we have P Q =0 

For r = nwe have Pe"“ + Qe~"“ — E 


.'. E = 2P sinh no. 
jj, sinh ra 
*'■ ~ sinhna 

-- li II . . II 

. (296) 

II 

II 

rOi 

=0 = 

ir~2 

II 

r 

*?« 

Fig. 27 


The voltage across the rth disk is 

V g = 6f ^r— 1 


or 


Vr 


2E 

sinh not. 


[sinh \oL cosh a(r — ^ )] 


Vr cosh (r — ^)a 
Vr + 1 ~ cosh (r + J)a 


(297) 


If now Zi and Z 2 are represented by pure capacities, then 
2 cosh a = C' 2 /^i + 2 


i.e. a real number. 

Under these conditions the voltage would taper off exponen- 
tially from the line side to earth. In practice it is found that 
this is not so ; the voltage curve rises on the earth side. This is 
readily explicable when it is realized that Zi should include 
surface leakage due to dirt, moisture, and local surface ioniza- 
tion. Extant literature gives no definite formula to express aU 
these quantities simply, so we take 


“ 1 +j<oC,E 

where J? is a resistance shunting the condenser between pins. 
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Now cosh a = 1 + i — ; — s 

1 + j(oCyR 

so that a is of the form A + j-B and the distribution is no 
longer exponential but may have a definite minimum. This is 
the distribution which holds in practice. 

8. Distribution of Voltage along a Transformer Layer Wind- 
ing. In a transformer winding, non-uniformity of voltage 
distribution has long been recognized, and has resulted in the 




dx 


use of reinforced insulation for the end coils and of a solid ring 
for the end coil. Recent work, however, has thrown some 
doubt on the efficiency of both these methods. Due to the 
complexity of the problem, we make simplifying assumptions 
to bring our discussion within reasonable limits. 

In any single-layer transformer winding the equivalent 
circuit could be represented as in Fig. 28, where we have the 
following quantities — 

(i) Resistance of winding. If r is the resistance per unit 
length of the winding, we have an elemental resistance of 
r . dx. 

(ii) Inductance of winding, due to 9 S (the total flux in the 
core, i.e. the useful flux plus the leakage flux, the latter varying 
along the winding). 

(iii) Capacity of winding to core. As the core is usually 
earthed, this will be the capacity to earth. If C is the total 
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winding capacity to earth, then the elemental capacity will 
be {Cll)dx = Cl . dx. 

(iv) Capacity between turns. This will be a series arrange- 
ment, and so, if C is the total self-capacity, the elemental 
capacity will be Cl/dx = C^jdx. 

(v) Conductances existing between coils, and between coils 
and earth. These we represent as g and 0 respectively. 


( 3 ^ \ 

Equations. Let e and ( e — ^ .dx\ 

be the voltages at A 

\ / 

and B respectively. Then we have 


{0 Gip)e . 

. (298) 


. (299) 

de 

. (300) 

^^ = ^Tx ■ ■ 

de n d(f> 

^ ^ 

. (301) 

d 

^*^dx + ^2 + * 3 ^ 

. (302) 


where n is the number of turns per unit length of winding. 

In equation (301) the value of <f> may be split up into 

(i) Mutual flux linking the entire winding and invariable 
with respect to x. 

(ii) A flux {(f) g) linking part of the windings and variable 
with respect to x. 

(iii) In a multi-layer winding, the variation of the flux 
linkages across a section must be taken into account. For a 
single-layer winding* 


Q-4fn{l.m.t.)n 


i>m + 


w 


n iidx,dxo 
. ' ^ ' 


where W = equivalent air-gap length, 
l.m.t. — length of mean turn. 


♦ Bewlky: Trans. A.I.E.E. (1932), VoL /j 1, p. 299. 



102 


TRANSIENTS IN ELECTRIC CIRCUITS 


On differentiating equation (301), we get* 

^ ^ ^ W X 108 diTdi 

_ L 

~ ^ aa:* Pdx.Bf 


where L = 


0-4:7rn^{l.m.t)P 

If X 108 • 


. ( 303 ) 


Substituting equations (298), (299), and (300) in equation 
(302), we have 

At*! = *4 A(i2 -f- t's) 

= (G* + Cip)e— A{gAe + C^Ae) 

d 

= (C? 4- Gip)e — g^H — G^b^e, where A = ^ (394) 

So we have equation (303) in the following form — 

A*e = i?[(G< 4- C'ii))A*e- (g 4- G^)b!^e\ 

- + Gip)e- (g 4- G^)A^] 

On collecting terms we have 
A^l 4- li(g 4- G,p)e - A\R{Q + G,p) -pp{g + G,p)] 

+ 4- C'ijo)e == 0 . (306) 


From the general equation (305) we are able to deduce the 
voltage distribution in the winding. 

Initial Voltage Distribution. Initially the voltage distri- 
bution will depend solely on the series capacity, so that we have 



CO 

<3 

II 

. (306) 


*8 == CispAe 

. (307) 


*4 = • 

. (308) 

Hence 

Cgi^A^e = Gype 


or 

Q 

b?e = TT e = ni^e . 

. (309) 


* Blvhb and Boyajian: Trans. A.I.E.E. (1919), Vol. 38, p. 577. 
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The solution of this equation is given by 

+ -Be-”*® . . . (310) 

ig = -Be"”*®] • • (311) 

where A and B are constants dependent on terminations. 
If the neutral is earthed through an impedance Z(p), we 
get at a: = 0 e = E 

and a,i X = l e — . Z{p) 

So we have that E = A B 
and iaZ(p) = Z{p) . pC 2 ni[Ae^* Be- — Ae”*^ + -Be-”** 
whence 



Hence 


[Z(p).V(C'iC'2).p +!>-”“ 

X [Z(p)'\/{GiG2) -P cosh ml — sinh mi]-^ 

[Z{pW{G^G2).p + l]B’^^ 

X [Z{p)\/ {CjC2)p cosh ml — sinh ml]-^ 


e^ = E 


^iP)V ■ P *5osh m{l — x) — sinh m{l — x) 


(312) 


Z{p)y/G^G 2 ) . p cosh ml — sinh ml 
when Z{p) = 0, we get the voltage distribution, with the 
neutral dead earthed, as 

sinh m{l — x) ^l — x 


E 


sinh ml 


I 


E 


(313) 


When Z{p) = oo , we get the voltage distribution, when the 
neutral is isolated, as 


„ cosh mil — x) ^ ^ 

E r — i — —B 

cosh ml 


(314) 


It is seen that from equations (313) and (314) that the 
voltage distribution along the winding will be exponential in 
form, ranging from 0 to ^ in the case of the earthed neutral, 
and from ^ to ^ in the case of the isolated neutral. 

Final Distribution. The final distribution of voltage will 
depend on the resistance and conductance elements of the 
winding. Now we have 

h = G'e = A(ii + ta) = Ac + grAc^ 

.’. Ah — p e — nh .... (316) 

1 + gR 
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The solution of the equaton is — 

.e = -f 

* = ^ Me”* — Se-”«] 

K 


( 316 ) 


If the transformer is earthed through an impedance Z{p), 
we get 

^ Th 

^ cosh nij, — x) — sinh n{l — x) 


E^ = E 


n 


( 317 ) 


( 318 ) 


Z{p) ^ cosh nl — sinh nl 

So for a grounded neutral we have 

_ „ sinhn(l-x) 
sinh nl 

As is small, we may approximate : thus 
E, = E\\ - (x/0] 

For the isolated neutral we get 

w - T. coshn(Z-a;) ^ 

^ cosh ni ~ ^ 

Thus the final state of voltage distribution is a straight line 
sloping either from ^ to 0 or 1? to E, depending on the imped- 
ance in the neutral connection. 
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CHAPTER VIII 


BUILDING-UP OF VOLTAGE AND CURRENT IN A LINE : 
REFLECTION PROBLEMS 


General. A mental picture of the building-up of voltage in a 
line may be obtained from Section 4 of this chapter.* In the 
previous chapter we have seen that the expansion theorem 
gives a Fourier’s series which converges very slowly and is 
useless for computation. We now investigate the problem 
using operational equivalences, i.e. Fourier’s integrals. 

1. Infinite Line. In Vol. II, p. 382, Heaviside gave an 
interesting development for the case of the infinite line. From 
Section 4, Chapter VII, we have seen that the current for a line 
earthed at the far end, when G* = 0 is given by 


^ ( n7ix\s\n ^„t'\ 


The step in this summation is wtt/I = w, say. In the limit 
we have Ttfl = dm. The equation now becomes 

sin Pnt 


h = ErC\ 


cos mx 


i»» 


where — p^) 

To evaluate the integral, we note that 
d d dPn 

d^)(co® == ^ ^ d(^) 


dmj 


t sin X i . 




¥ 


^ ^ d sin /?„< ^ ^ 

But ^ — = cos . 

So we have from (322) and (323) that 


sin 


d 


cos = ¥ 


sin 


d{p^) — **’ /3 — 

j d dl 


(320) 

(321) 


(322) 

(323) 


( 324 ) 


♦ Also, see Heaviside ; Electro ^magnetic Theory, Vol. II, pp. 67-77. 
8— (T.132) 105 
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Expansion by Maclaurin’s theorem gives a result of the form 


fip^) ==/(o) + pmo) + + • • • 

where/' (0),/"(0) . . . denote differentiations with respect to p*. 

and so = 

and by ( 324 ) we have 


Here 


mv 


, 1 < sin mvt 

f (0) — — o , etc., 

^ ' p 2 mv 


or 


sin 


} + ^‘^2 + f ‘ 


Y 1 < Y ' si 

!Vp2; + •••_" 


sin mvt 


mv 


or written in full — 
sin sin mvt 




mv 


+ 


*4 


t sin mvt 


I 


+ 

So the integral has the form 
sin mvt ^ 

cos mx dm 

0 mv 

By tables — 


mv 

2! 


dt 


, I ‘ i sm mvt 

tdt\ 

Jo Jo rnv 


dt -f~ . 




sin m{vt + a;) + sin m{vt — x) 


mv 


] 


dm 


I 


Sin mx , 7 T 

dm = 7: 

0 m 2 


So that the integral given above will have the value 


and 


+ 1] if v<> X 

^[1- i] if «< < X 


To avoid the zero or singularity we will start our integration 
at a time t = xfv. We revert to the physical aspect of this 
step later. Now the integral becomes 
sin p„t . 


r 

I cos mx ^ 

Jo Pn 


-dm 





BUILDING-UP OF VOLTAGE AND CURRENT IN A LINE 107 


r r 

Term by term, the integral becomes ^ j> ^ I » 

Thus a series, known as a Bessel series, emerges from this 
integral, or 


r sin/3„<^ ^ rf f ^ 

cos mx —5 — dm = ^ Io\ pi s ) I 

Jo /?« 2v vVj(« >*/.,) 


So that the current in the line becomes 



-J'W'-S)']).. 


*/») 


. (326) 


Physically, no current can exist in a line until the impulse 
has had time to arrive, or < = xfv. Once the impulse arrives, 
the current will build up according to the above expression. 
For example, at a: = 0 the current commences to build up from 
the instant of switch closure as 

but at a distance x along the hne the current is zero until a time 
xfv = 0, say, has elapsed. Thereafter it will grow as 

E,y/{GIL)e-^UJip{t^- em 


2. Operational Development. This result could have been 
developed by direct operational methods as follows — 


For an infinite line 


■-L l(t±3^\ 

* Lv^ \p -\- 2a / 


, . . , • V l/p + 26\,,, 

or at the transmitting end *o = J^J^ p _|_ 2a 

for a steady applied voltage. 


V 

But -F- 


(p + 26) 


■(^) = 




Lv ■\/[{p -f- 2a) {p -f 26)] ■ ' ' Lo ^/{p -f pf- 
Now using the shifting theorem we have 


a^] 


(326) 

(327) 

■il) 


V 


VKp + pf - 


or 


: = 


p 






_F 

I/o' 


1 + 


= e-p* 

Q 


P 


Gp 


ViP^— 
je-'>‘/o((rf) . 


V(P^-<r^) 

= e-'^Ioiat) 


(328) 
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Case 1. With G = 0. Here p = a = a, and we get 

Lv 
V 


V 

and “ Zw 

where 


. (329) 
. (330) 


a = ^ -\/(P® + 2oi)) = I ViP^ + 2<T2)) 


Shifting to the left, we have 

Power expansion yields 


. (331) 


Lv 


■i‘+s 


2 p2 — q.2 ^4 ^p2 — Q.2J2 

+ ^ 


...] 




(p^ — (T^ )® 


Now 


* 3! 
P 


2! ' V* 4! 

• ••][; (p“- '^“)‘]) /.(rf) 

(332) 


5! 

. (i) =/o((rO 


^/(JJ^ — <72) 

ViP^- = p.(l) = 0 


So that all terms in the second brackets are zero, and we are 
left with the first set of terms. 

From the Bessel equation (see equation (454) of Chapter XI) 
given by 

p^Ifiiat) + -^pIo{ot) — aH^{at) — 0 
we have (p® — <t®)/o(<tO == — |/o(o^O = — “7 


at 


and, generally, 

{p^-a^)IUat) = - (2m + 1) {ot) 

On substituting, we get for (332) 

AM) 


This expression may be shown to reduce to 


’• = x; 




>xlv 


. (333) 

. (334) 
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Comparing this result with equation (329) for shows that, 
given the value of current at the transmitter, we may write 
down the current at x by changing the argument of the Bessel 

Heaviside gives several other methods of 


to 


/('-a 


developing this equation. 

Case 2. With all constants present. We now have 




and 


l.=.vv(^C + fjp. 


—VHP + P)*-o*] 

e ^ 

+ pf — o‘ 


•U) 


(336) 

(336) 


Now, by Campbell’s table {loc. cit.) No. 866 — 

P • . /r... ^ .2^ _-^a] • (^) ^ h - -a j J (337) 


So 


V[(?> + 9^) 


+ 


«G»Fo|V'’^/o[or(A2-^*JJdA . (338) 


and Vrr 


V 


- pxlv ^ 


ax 




^alv 




h 

a\ 

1 

G 

Vi 

1 

1 


dX, 


(339) 


I > xlv 


3. Finite Line with Tenninal Impedances (Wave Solution). 
Let a steady voltage E be impressed on the line through an 
impedance jZj, and let the line be closed bA x — I through an 
impedance Z^. On substituting in the general line equations — 

F* = -f Be-**, 

”p -f 26 ' 


and 




z, 


we get [Be-“* — yle“*] = .4e“* + Be-^ . 

^<p) 

A-\-B = E-^{B-A) 


and 


(340) 

(341) 
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With 


] 


= Pi and -fy — = P 2 , it is easily shown that 

(P2- l)g-2“' 


'(») 

^ = 


5 = 


(Pi+1)(P2+1)-(Pi-1)(P2 

(P2 + 1) 


l)e 


-2a) 


E 


{p^ + 1) (P2 4- 1)- (Pl- 1) (P2- 1)6-^' 
On substituting, we get 

P2 ^ -a(2)-a:) 






- z, + z,., 


P 2 + 1 


1 - 


Pl 


1 


P 2 


1 


Denoting 


Pl 


I — Pl + 1 P 2 + 1 

1 


-2a/ 


(342) 


Pl + 1 


^ by — Pl and - ^ 


P 2 + 1 


by — ^ 2 > we have 


E 


fe “® + P 2 ® ***^ *n 

L 1 — /«iP2e~^' J 


(343) 


z, + z,„ 

On expanding we get 
h = y - j-y - [e-“^ + 

"1 “*■ -^(p) 

+ + . . .] 

For an infinite line fed through an impedance Zj at a; = 0 the 
current at x is given operationally by 

E 


(344) 


Zi + Z, 


(v) 


So the first term in (344) is the current at x due to E impressed 
at a: = 0. Similarly the second term gives the current at 
(21— x), and so on. The first term gives the effect of the direct 
wave, while the second gives the effect of the wave reflected 
from the far end, and the third that of the wave reflected from 
the near end. The above equation is thus the sum of direct 
and reflected waves which are effective at times t — xfv, 
t— (21— x)/v, t = (21 -j- x)/v, etc. Thus if we require to know 
the current at a time less than (21 + x)/v, we should take into 
account only the first two terms of equation (344). 
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4. The Voltage after Infinite Time. In practice the usual 
problem is to find the voltage after an infinite time. Let the 
voltage E be impressed on a distortionless line at C. The line 
will attenuate as and have reflection coefficients a and § 
at A and B respectively (Fig. 29). For a voltage wave moving 


ccAI“ 

1-. — a- 


I 

c 


X 

Fig. 29 


BP 


to the right we have that, at x, the wave is v■^^ = 

At B the value of the wave will be and the reflected 

wave is ^ When this wave passes x, its value is 

/3 . . £-««-*) = p . ^^£-9(2«-a-*) 


So at A it will have a value of . \Ee~‘^^^~”'\ and the reflected 
wave will be a.^\Ee~^^^~"'K This is attenuated to 


at X. Now, at x, the value of the voltage will be the sum of 
the voltages, thus 

F* = \E [£-«(*-“) 4- ^e-Q(2l-a-x) _pa^e-«(2<-a+*)^ _ _ ] 




£-«(*-“) [1 + / 5 £- 29 «-®)] 


1 - oL^e 


-2ql 




. (346) 


In a similar manner the waves travelhng to the left 
voltage 

gjg-4(* + 0) |-1 _|_ 

1 - J • 



will give a 
. (346) 


So that the total voltage at A is given by the sum of (345) and 
(346). We thus have 


V 




a£-«‘*] [£«<'-*>+ /3£-«('-*)] 


cat . 


a^e 


-Ql 


] 


. (347) 


A convenient method of following through any voltage 
impulse is by means of a lattice diagram shown in Fig. 30. 
This has been constructed for the pulse travelling to the right 
only. The pulse travelling to the left would need to be super- 
imposed on the diagram. From a knowledge of the length of 
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line the times at which the various pulses arrive could be 
determined, and so the voltage at any point could be found 
by the addition of the pulses at their correct times. This is 
discussed in Section 9 of this chapter. 

From equation (347), by substituting different values of a 
and j8 we are able to write down the values of the voltage at x. 


A _ C B 

^ CL ^ 



For example, with a line open at both ends we have 
a = == + 1, and (347) reduces to 


F = 


E 


sinh qa . cosh q{l — x) 
sinh ql 


. (348) 


When both ends are earthed, then a = jS = — 1 and the 
equation becomes 


F = 


E 


cosh qa . sinh q{l — x) 
sinh ql 


. (349) 


With an earth at A and the line open at B we have that 
a = — 1, /S = 4* li and so 


V = E 


cosh qa . cosh q(l — x) 


(360) 


cosh ql 

These expressions give the steady-state solutions and confirm 
results given in Section 3 of Chapter VII. 

6. EflSciency of Reflectors. In Section 3 we have seen that 
the reflection coefficients at the far and near ends of the line 

^(j>) — -^2 j ^(p) — 


are given by 


Z 


iv) 


'(d) 


+ 


respectively. 
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where Z(j,) — \/{ZIY) is the operational characteristic im- 
pedance of the line. 

For steady-state conditions in a line terminated by an 
impedance at x — I we have the following equations — 

E, = p,[l -f 

The incident voltage is 

1 4- . • • (351) 

and the reflected voltage is 

. . . (362) 


At a: = ? we have Eg — Vi-\- v^. The current equation may 
be written as so that the incident 

current wave is 


— ^ 0^1 ' 

and the reflected is 

*2 ~ — • 

The efficiency of the reflector is 

^ ^0 H 

Zj^ -f- Zq 


. (363) 
. (364) 

. (366) 


We now apply these results to some simple line conditions as 
follows. 

(i) Line Open. In this instance Z^ is infinite, and Zq is 
negligibly small. From equations (351) and (352) we get that 
the incident and the reflected voltages are each equal to ^Eg 
and are both of the same sign. We may write these as Uj = 
^IgZi and V 2 = \I^x- From equations (353) and (364) we 
get that the currents are = \EgYf, and *3 = — ^EgY^, so that 
the voltage at the end of the line is the sum of the incident 
and reflected volts, and the current is zero. 

(ii) Line on Short Circuit. Here Z^ = 0 , so that Fj = 00. 
The value of the incident and reflected voltage waves will be 
given by = \ZfJ[g and = — ^Z^Ig. The current waves are 
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given by ii = \ and ij = ^YffEg. So, for a line on short 
circuit, the voltage at the far end is given by 

En== Vi + Vi = 0 


and the receiver current is 


2*1 — ^sYq 


(iii) Harmonics in Wave. For a wave of composite fre- 
quencies, in the case of a non-dissipative line, the reflection 
coefiicient is 

o ^in -^Qn 

P” - X,, + JCo« 

and thus is independent of frequency. Thus no distortion 
occurs in the reflections. When resistance effects are present, 
then we can no longer consider that = . . . — and 

distortion will occur. The generalized problem becomes very 
complex. 

6. Series Circuit. For an overhead hne and cable in series, 

^ 2 ^ jgr 

by the method of the last section we have 6 — ~ = J ;■ as 

Vi Zj 

the reflection factor into the cable. When is the incident 
volts, the voltage at the junction becomes 


2Z(. "] 

V = v, + vJI = 


This is the voltage transmitted into the cable. The voltage 
reflected into the overhead hne is 


The incident current is 


r z, -) 


and the reflected current is 

i --1* 

*2- Z, 


7 ^0 + Zi 

2Z.Z, 


y 

• 2Z.Z, 


The current transmitted into the cable will be 


*1 ■ - *2 — Y . 


2ZcZi Zg + 
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For a wave travelling from cable into line we have 

vi Zi — Zc 

and V = Vi 

So that the incident voltage becomes 



V ► V 

oCjV {l*ci^)V ^a+oLiWi 

(Uoc,j)v — a+(ti)Vi 



and the voltage reflected back into the cable is 



Similarly the incident current wave is 

< vz, + z, 

“ Ze 2 Z^Z, 


and the reflected current is 



Zi-Z, 
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2*1^ ^ 

So the transmitted current is given by ^ • 

Summarizing, we have the following coefficients- 


From Overhead Line to Cable 
Reflection coefficient 


Zi-Z, 


Zl + Zg 

From Cable to Overhead Line 
Zg-Z, 


; transmission coefficient — 

1 + aj. 


Reflection coefficient a, = 75 — ; — ^ ; transmission coefficient = 

1 + a.. 


We are now able to study the case of a voltage fed into an 
open-circuited line, i.e. when = 1 : as shown in Fig. 31, the 
line losses are neglected. The voltage on the cable builds up to 

2(1 -j- ai)F[l + aa 4 - «2^ + • • •! = 2(1 4- ai)F . - — 

1 CX 2 

As aj = — a. 2 , the final voltage at the junction is 2E. 

7. Series-parallel Circuits. Fig. 32 shows a typical arrange- 
ment consisting of a line of surge impedance Z feeding on to a 



transformer of lumped impedance Zj. From the other side, lines 
of surge impedances Zg and Z 3 are taken. The impedance at A 
looking into the network is 

ZaZa 

The reflection coefficient is 

Z^-Z 
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If Vi is the incident voltage, then the reflected voltage is 


The potential at A will be 

Vi + 






The current through Z will be 


Zft Zf^-\- z 

The potential at B will be equal to the voltage at A minus the 
potential drop in Z^ and so will be 

Zo + Z Z^Z^ + Z^^~ z, + z 
Further, the current through Z^ will be given by 

Z,Z, + Z~ + Z, Zo + zj 

The voltage drop across the whole is given by 


iZi + igZg — 


/ ZiZ2 Z^Z^ -f* Z^Z-y 

Z 4" z\ ZoZn 


To fix ideas let us consider a line of surge impedance Z^ con- 
nected to an inductive coil which is in series with a line of 


Fig. 33 

surge impedance Z^, the junction being earthed through a 
capacity G (Fig, 33). With the above notation we get 

r7 ^ ^2 Z^LCp^ + Lp + Z^ 


Zq — Lp + 


1 “h 


1 + 


and the reflection coefficient at A is 


^ ^0 

^0 + ^1 
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The potential at ^ is v (1 + j®) and the current through Lp 

is ^ (1 + jS). The voltage at B is 
■“0 

(1 + i?) - (1 + /3) i>2> = vx (1 + /3) ( 1 - 

the current through the condenser is 

^. = C7p^;x(l+/3)(l-^) 

So that the current into the line is / — 

It is seen that each individual equation is in operational 
form; for evaluation, the impressed voltage wave form must 
be known. 

8. Effect of Lig h tning Stroke. We will consider a lightning 
stroke on a line terminated by an inductive coil. The consensus 
of opinion seems to be that the voltage surge due to lightning 
may be represented as Ee~^K Here the incident voltage, is 
equal io Ee~'^ . (2), and the reflected voltage is 



By shifting, we get 


L (p — g) — ^0 
L (p — a) + Zq 



The expansion theorem gives 
Ef 

"« = 3r-~z „ + 2.) - 2 v-"'’"] 


The potential at B is given by 

■ f „ [2aL- 

qlJu! /!/q 

and the current through the inductance is 

2v, 2E 

Z, + Lp • Z^+{Lp-a) • 
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If instead of an inductive earth, a capacitative earth had been 
used, then the current through the condenser would have been 


2vi 

1 -f- Zffip 


Cp . ( 1 ) 


or 


2ECe-<^ 
a.Zfp — 1 




9. Overhead Line Terminated by a Cable. A length of cable 
is frequently used to absorb surge voltage and so protect 
machines connected to the system. In the cable system shown 


Zi 

O.H.h. 


A Zz B 

-I I- 


Z3 

Ma(Jtines 


Fio. 34 


in Fig. 34 we will investigate the voltage occurring at the cable 
connection. The effect on the wave form is neglected. 

Neglecting Attenuation and Time Effects. The Vol- 
tage AT B. When an impulse voltage E is applied at A, the 


2Z, 


voltage transmitted to the cable is — r^E. This voltage is 

Z’l + ^2 

reflected at B, and the voltage at B is given by 


2^2 ( 
Zi + Zj\ 


1 + 


Zs — Z^ 
Zg + Z2 


) 


E 


and the reflected part is 

2^2 Z^ Za _ 


On the reflected voltage reaching A it is reflected again, and 
this component is 


2.^2 Z^ Zn Z-i Zo ^ ^Za ,, 

^1 + ^2 ■ ^3 + ^2 ■ ^1 + ^ 2 ^ ^ ^1 + ^ 2 ^ 


where M = 


Zi — Z2 

Z\ + Z2 


Zo Za 


E. 


Earthing both ends of the system, we can concentrate on the 
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wave within the cable. The value of the wave is now 
moving towards B. At B the reflected voltage is 


2 ^a 


M 


^■^2 -tr ^3 -^2 


and the total voltage at B is 


2^2 ( 
+ Z 2 \ 


1 + 


Zo — Zn 

^T^2 


) 


M 


Similarly for the next stage, the reflected voltage at A is 
given by 


2^2 
"f" ^2 


^3 — ■^2 


•^2 
“t" ^2 


M = 


2^2 

^X + ■^2 


Jf 2 


and the voltage at B is 

2 Z, 


M^[ 1 + 


^3 


Zo . 


•^1 + ^3 V "3 

So the total voltage at B is given by 


Zn + Z ^ 




2 Zo 


Zx "I” ^2 




+ if + ilf2 + . . .) 


= E 


( 


1 + 


Zo — Z-t z» — z^ 



m 


where mjE? = Af and m<l. 

When Zi — Z2 = Z^, then E^ = E. 

CoNSiD BEING TiMB EFFECTS. The Voltage AT B. Let the 
time for a wave to travel twice the length of the cable be T : 
reckoning time from the arrival of the first wave at B, we have 
that the initial voltage at B is 


E 


I 

Zx + ^ 2 \ 


1 + 


= EPs" (say) 


The second arrival will be 


The third will be 
and so on. 


EinPe~'*^ 


Em^Pe~^’^ 
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If we sum for n steps, i.e. for nT seconds, the efiFeotive value 
of the first voltage at B (see Fig. 35) is given by 


the second by 
the third by 


E .P . e-»r(n-i) 






Volts 

r-JT 

Q 

N 

_ 1 — 

Tiwe 

m 

h i^-DT J 

P 


Fiq. 35 


and so on. So that the total voltage is 


With a = 0 we get 




As m is less than unity, the value of E may be several times 
that of Eg. For = Z^ = Z^ we get Eg = E. 

For the End A. The initial voltage at A is 




This pulse is partially transmitted and partially reflected. The 
transmitted part is reflected again from B, and on reaching A 
the voltage is given by 


E 


^Za Zo Zn 


( +1+^) 


or 


+ ^2 -^3 “I" ^2\ + '^2> 

2^2 / ^3 — *^2 1 \ 


9— (T.152) 
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2Zo 


The reflected part moves to B as — --fjr m. 

Zii + ^2 

So that the total voltage at A is given by 


E 


2Z, 


^3-^2^ a' 


+ -28 + ^2 
Summing for n steps, we have 

E^ = + Pe-(«-i)«r 


+ P 

E + Pe-(**-i)«T j’ij^ 




■^3 “ 1 “ ^2 
1 — 


] 




where P = 


2Zo 


^1 + ■^2 
With a = 0, we get 


And when 2^ = = this reduces to E^ — 2E. 


10. A Condenser Discharging into a Short Length of Line 
on Open Circuit. Let Z be the impedance of the line and 
l/Gp the operational impedance of the condenser. Then, at 
the condenser terminals, we have the reflection coefficient as 


l/C'p — Z _ g — p _ 
1/Cp -{■ Z a. p~~ 


where a = IfCZ. The transmission coefficient wiU be 

1 t ^ „ 

■^g + p g + p 

The condenser discharge may be formulated as Ee~^, and, 
during the passage of the voltage along the line, let it be 
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attenuated to Ee~^ . e“** = ae. 
marized as follows — 

Condenser 
Charge Discharge 

E ... e^ 

a*ei 

I 

I I 

na^ei 

ma^e^ 

I 

1 I 

n^a^e^ m^a^e^ 

etc. 


The action may be sum- 
Far End of Line 

aci) 
aei) 
ma^i ) 
ma^ei) 
rn^ahi ) 
m^a^ei) 

etc. 


Then by summation we have that the voltage at the end of the 
line is 

^ r o o i: 2ae, 

2 roe, -I- mo®ei -4- inraH-, . . .1 = ;; ; 

^ ^ ^ ■' 1 — mar 

The residual voltage on the condenser is 

E — e^ — a^mCi — ahn^Ci, . . . -f ahiCi -f- a^n^e^ -f- . . . 

„ , r m—n 1 

— Cj ® _ o2^) (1 — a2^)J 
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CHAPTER IX 

EDDY CURRENT EFFECTS 




1. Fundamental Considerations. Closely allied in physical con- 
ception with the diffusion of electricity along a wire is the 
penetration of flux within an iron plate and problems such as 
the electrical resistance of copper bars to currents of high 
frequency. The fundamental equations of these problems are 

based on Faraday’s and Ampere’s 
/I laws. These are derived as 
follows — 

Fabaday’s Law. For an iron 
plate as shown in Fig. 36, let the 
flux be acting downwards as 
shown. For an elementary area, 
of length 1 cm. and width dx at 
right angles to the plane of the 
flux, the total flux through it will 
be B .dx where B is the flux 
density. The e.m.f. induced in 

this circuit ABba ia — zr. B .dx. 


d 

f 

\lcm. 

♦ 

A 


Ct C I 

I ^ 

I ^ 


\dxi 


-2a 


dt 

If the current density at A is a, 

da 

then at B it will be u -4- dx. 

dx 

The potential difference round 
this element will then be 

(“ + 5^* = '>55*’ 

where p is the resistivity of the iron and the voltage drops 
across the ends AB and ab are negligible. So we have that 



Fig. 3G 


-^jiBdx) 


da 


Pd 


X 


dx 


(366) 


AMPiJEB’s Law. For an elementary area ABcd parallel with 
the plane of the paper, of magnitude 1 cm. X dx and with a 
current density of a, the total current through the area is 
a . dx. Hence the work done in taking unit magnetic pole once 
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round this area is 47ro . dx. When H is the magnetic force at 
3H 

A and H + work done is 


So that we have 


4:TTadx = 


1 3B 

^dx . . (367) 

(I dx 

where p and p, are in e.m. units. 

2. Resistance of a Round Conductor. On neglecting the field 
set up by the eddy currents, the resistance of a conductor 
may be found approximately by the following method, due to 
MouUin. The complete solution is given in Chapter XI. In a 
conductor of length I and cross-sectional area ttu®, let the 
current be /q sin wt. For low-frequency alternating current, 
the current density will be uniform and equal to Uq = 

At high frequencies the distribution is not uniform, because 
the induced e.m.f. is sufficient in magnitude to produce distor- 
tion. The density wiU be given in the form 

or = (To ± 

where a* is the density of the eddy currents. 

As the eddy currents do not appear in the external circuit, 
then S(T« = 0, i.e. the eddy currents must fiow in opposite 
directions on each side of some radius r^. Let us consider a 
section abed at right angles to the magnetic flux with be on the 
neutral radius as in Fig. 37. As the magnetic field strength 
H within the conductor increases uniformly from zero at the 
centre to {21 Ja) sin cut at the external surface, then the flux 
through the rectangle will be 

, , f’’* 2/0 . , ,, o ov^sineoi 

ad -Hdr = ad I — sm (otdr = ad{r^ — r ^)— — ^ — 

Jr ^ Jr ^ d 

So, Faraday’s Law may be written. 


^ ^ a* 


But /fl = Tra^cTo ; so 


(T* = — ro®)<7Q cos cot 


(368) 
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As 


J” 2-nra^ 

I 


= 0, so 


a 1 

-2TT^(x)r{r^ — ro^)ffo cos (otdw = 0 
0 P 


On eliminating factors independent of r we have 

^ _ ft 

4 2 ~ 



. ( 369 ) 


Tfj = al\/2 as the radius of the neutral cylinder. The 
eddy currents inside and outside this cylinder flow in opposite 
directions. 

So that from equation (368) 
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The resultant current density will be given by 


a = 


/ft • CO / r® 1 \ , 

— ^ sin ft)< -I — ( -5 — x )/ cos oit 

■na^ ' p \a^ 2 / 


. (360) 


i.e. o is the sum of a uniform density = lalTra^ together with 
a density a, which is a function of the radius : where r < 0/^/2 
the eddy current is in phase opposition to a^, and where 
r > a/\/2 it is in phase therewith. The current is crowded 
towards the outer surface, giving a skin effect. 

Resistance. The mean power loss at low frequencies is 
given by where R — plfna^. At high frequencies the 

mean loss of power is defined as \RJ^, where R, is the effective 
resistance. After some calculation we obtain that the mean 
power is 

^ + "oV Jo 

/ a* \ 

So 


and 


^ _ 1 , ^ 

R ~ 3p2 


. (361) 


3. Eddy Cuirents in Laminations. We will assume that the 
permeability is constant for the range of flux densities encoun- 
tered and that the hysteretic angle is negligible. With B — pH 
we have, for Fig. 36 — 

- dx~~ dt ^ dt 


Combining these, we get 

With H = we have 


4770" 


dj^ 

dx 


4 it^ dH 

~V "aT 


dx^ 


dm _ 

dx^ ~ 


(‘iTTplp)j(oH^^‘“^ = 


. (362) 


where = {inp(olp)j or m = ['^/(^TTpoilp)] (1 -f i)- 

This equation is of the same type as that for a transmission 
line, and is integrated as 

H — -f . . (363) 



128 


TRANSIENTS IN ELECTRIC CIRCUITS 


where A and B are constants determinate from boundary 
conditions. 

For a plate of thickness 2a the centre line may be regarded 
as an axis of symmetry. So that quantities at a: = o and 
X = — a must be identical. Thus A = B, and if is the 
amplitude of H at the surface we get 


so 


g(l+j)mo _j_ g~{l+j)ma 

_ tr / cosh (1 +i)OTa; \ 
®\cosh (1 j)ma ) 


. (363a) 


Now cosh (1 + j)0 = cosh 0 cos 6 j sinh 0 sin 6 

and on rationalizing we get 

VCKcosh 26 + cos 26)] /<}> 

where tan ^ = tanh 6 tan 6. 

Substituting these results in equation (363a), wo get 

,, „ / cosh 2mx + cos 2ma; N* , , . v 

^ ^«(cosh 2ma + cos 2ma ) 

where tan <f>i — tanh mx tan mx . . (366a) 

tan <^3 = tanh ma tan ma . . (365b) 


When the external magnetic force is given by Hq cos o)t, 
then the value of H is 

rr „ fcosh 2ma; + cos 2ma;1i , ^ ^ 

= ^» Ush2Ma + cos2ma J >''> ' 

where = — <f)^. 

At the centre line we get 


H. = H, 


•( 


cosh 2ma + cos 2ma 


)‘ 


cos {(Ot — ye) 


(807) 


i.e. along the centre line, H has its minimum value. 
The mean value of H is given by 




1 ra 

ajo‘ 


Hgdx 


^n.= 


Hq ( cosh 2ma — cos 2ma ^ 
rmy/2 \ cosh 2ma + cos 2ma ) 


(368) 
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Space Phase Angle. The angle y* will represent the dis- 
placement of H m X is varied from a to 0, or as we move from 
the outside face to the centre. Now 

tan y = tan 

and by expansion we get 

_ sinh m{a — ») sin m{a + a:) + sinh m(a -f «) sin m{a — x ) 
^ cosh m{a — x) cos m(a + + cosh m(a -f- x) cos m(a — x) 

. (369) 

At a; = a we have tan (^i — ^ 2 ) = 0, i.e. the angle of lag 
is zero. 

At a; = 0, then 

tan ye = tanh ma tan ma = tan 

Now, when (<^i — ^ 2 ) == then must be in phase with Uq. 
From equation (369) we get 

mx = ma — 2stt* 


or 


X 


a 


-.2rr 

m 


At depths given by 27 t— from the surface, or at depths of 

27T/m, ^rrfm, . . . the magnetic force is in phase with that at 
the face. 

When (^1 -f ^a) = 0, we have that at depths rr/m, drrfm, . . . 
the forces are in phase opposition. The depth 2nlm is defined 
as the wavelength, or 


A = 


2tt 


\/(2lT/i(0/p) 


Vipll^f) ™ rafe-units 


or in practical units, with = 


2TTfi(0 


X 10~®, we get 





. (370) 


Effective Penetration. This is defined by Ewing as the 
thickness of a surface layer which, at a constant flux density 

* 8 — an integer. 
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£g, would give the same flux as exists in the lamination. Let 
the thickness be d. Then 



On integrating we get, for large values of ma, that 

, 1 

{l+j)m 

Or in practical units 

d = l/mV2 = 3 670V (p/^/) • • (371) 

Thus with ju — 1 000 and p = 10“®, at 1 000 o.p.s. the 
depth is 0‘0113 cm., and at 10® c.p.s, is 0-00113 cm. 

4. Eddy Current Density and Iron Losses. From equation 
(357) we have that 

__i.^ 

^ ~ i-n dx 


On substituting for H and differentiating, we get 


^=-S[ 


.. sinh (1 4- j)mx' 


or 




a = A sinh (1 -j- j)mx . 

1 


(372) 


where ^ - l-« (l + j)™ _ 

ra 

The loss in a plate may be formulated as (p a^dx) per a 

Jo 

The mean value of J (cosh 2mx — cos 2mx)dx 


cm. 


ML” 

4ma 


(sinh 2ma — sin 2ma) 


As the r.m.s. value of a harmonic quantity is (l/V^) times 
the maximum, so, over a complete period, the effective value 
of (T* is 

1^21 


8ma 


(sinh 2ma — sin 2ma) 


and the loss becomes 

(sinh 2ma — sin 2mo) X 10“2 watt . (373) 
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Substituting for l^^j, the mean watt loss per cm.* is 

Ho^pm ( sinh 2wa — sin 2ma \ 

— vTS 1 S — X 10-T watt 

327T*a \ cosh 2ma + cos 2wo / 

In terms of mean density, we get 

£o*a ^ sinh 2ma — sin 2ma \ 


W, = B, 


X 10~^watt 


mp V cosh 2ma — cos 2ma / 

For large values of ma this expression becomes 


W. 


IGwa 


V ipf4) X 10“^ watt 


(374) 


(876) 


5. Phase Angle between Current and E.M.F. in an Exciting 
Coil of a Closed Magnetic Circuit. The power input to such a 
circuit is cos 0, where V^ax and are the 

maximum voltage and current. This expression may be 
written as 

^ Volts/Turn x Amp. Turns . cos <f> 


So that the loss per cm.® is given in r.m.s. values as 


Ft X A/cm. 
area 


COB<f> 


Now, volts per turn 
density at face {B,,) 


(aB„ X area and itTfiAfcm. = flux 


// cosh 2ma + cos 2ma\ 
— (y^ )ma ^ 2^ j 


The input power is 


oiBJ- 

‘i'trfi 


/ V // cosh 2ma -|- cos 2rno \ , , „ . 


which is equal to iron losses given in equation (374). 

co®a / sinh 2ma — sin 2ma \ 

mp ’ ■ \ cosh 2ma — cos 2ma ) 


Hence 


cos^ 


/ // cosh 2wa + COS 2ma \ 

)fi>»M*y ^ _ cos 2mo / 

1 sinh 2ma — sin 2mo 


\/2 -\/(cosh® 2ina — cos® 2ma) 
In the limit, when 2ma is large, then 

cos ^ = W2 


(377) 
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6. Effect of Hysteretic Angle. The effect of hysteresis has 
been investigated by various authors. As a full development 
would require too much space, we give a summary of the 
results given by Latour.* 

Denoting the hysteresis angle by t, -y/ll + sin t) by a and 
y'(l — sin t) by we have the following relation — 

Bmax — (Apm/ft))'v/(co8h 2max + cos 2m/Sa:) . . (378) 


where A 


PY'(cosh 2in(x.x 

ft) 


cos 2m^a) 

4:TTflI 






~ pm ^/(cosh 2mm + cos 2mpa) 
( cosh 2mm — cos 2m^a y 
\ cosh 2mao + cos 2m^a j 




4mp 


'sinh 2mm 


a 


sin 2mpa\ 

]b„ 


I cosh 2mm 

'sinh 2mm 


cos 2mpa/ 


W, 


(o^a sin T 
4mp 


a 


+ 


sin 2mfia\ 




,cosh 2mm 


W. 4 - W, = 


COB <}> = 


ft)®a / « sinh 2mxa ■ 


imp \ cosh 2mm — 
(lf\/2){ot. sinh 2mm 


cos 2m^a 

^ sin 2m^a X* 
cos 2m^a / 
- (8 sin 2m^a) 


) 


B 2 


B * 

"TO 


(cosh^ 2mm — cos* 2m/9o)i 


(379) 


(380) 


(381) 


(382) 


(383) 


(384) 


On putting a = /3 = 1 in the above expressions, we obtain 
the previous results, which neglect hysteresis. 

Experimental Data. The following test results are taken 
from Dannatt.f 


Material 

Cos ^ 

60 o.p.s. 

H.F. 

Silicon steel 0*32 mm. 

0*68 

0*78 

Silicon steel 0*052 mm. 

0*18 

0*72 

Mu metal 0*26 mm. 

0*8 

0*86 

Rho metal 0*36 mm. 

0*6 

0*84 


♦ Latour: Electrician (1919), p. 219. 

t Dannatt: “Magnetic Properties of Ferro -magnetic Laminae, “ Journ. 
l.E.E. (1936), Vol. 79, p. 667. 
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Thus at high frequencies the power factors are greater than 
the limiting value of 0*707. From the same paper the following 
values of permeability are taken — 


Material 

Permeability 

60 c.p.B. 

H.F. 

Silicon steel 0*32 mm. 

880 

230 

Silicon steel 0*062 mm. 

480 

210 

Mu metal 0*26 mm. 

9 600 

180 

Rho metal 0*36 mm. 

1 100 

160 


Again the theory does not account for all the variation in 
permeability, so that it would appear that a complete theory 
must take in account the effect 
of the surface layer, the size of 
the iron crystals, and the direction 
of magnetization with respect to 
the crystal. 

7. Eddy Currents in Armature 

Bars. We will assume that the 
slot leakage lines go straight 
across the slot and that the iron 
is of infinite permeability. Under 
these conditions the value of the 
leakage fiux is directly propor- 
tional to the width of the slot, 
and, due to the pulsation of these 
lines, an e.m.f. will be induced 
in the conductor. The presence of this e.m.f. causes distor- 
tion of the current distribution within the conductor. How- 
ever, all leakage lines do not cause distortion. Certain of 
the fiux will be linked with the whole conductor and will 
induce the same e.m.f. in all sections. The effective fiux will 
be that due to currents below the section in question. For 
example, in Fig. 38 the fine shown as 1 is effective, whereas 
that shown as 2 is ineffective so far as conductor B is con- 
cerned. The total disturbing flux in the shaded element is 

A-I I’d 1‘z IA~ 

<f> = Ud-x)l+ 47r {Wadx)~ . (386) 

8 Jx Jo 8 


-w- 


B 


I 


% 


-X 

-0 




Fig. 38 
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where If, is the current in ab-amperes in all conductors below 
the conductor considered and whose flux passes 
within the region a: to d of the conductor B, 
a is the current density in ab-amperes, 

J Wcrdx is the current in the conductor B below the section, 

i.e. in the region O to x, 

I is the armature length. 

The voltage across 1 cm. length of the conductor is 

F = n + f 


when i is the current through the element. 

As this voltage must be the same for aU elements of the 
conductor, we have 

dV ^ di d d<f> ^ 

dx ~ ^dx ^ dx dt 

But ri = pla ; so, on substituting for i and <f>, we get 


I dh 4.-nl r* 


_ IMKJ ^ V -X.lt V I ___ _ 

ply-- If. - 77 . da: = 0 

^ da; a dt « jo dt 
With currents varying harmonically, we get 


dx ^ pa 


_ Atto) C* 


Wadx = 0 


Further differentiation gives 
d^cr .47roi 
dx^ ^ pa 


j Wa = 0 

'' pa 


(386) 


(387) 


where a* = If or a = (1 -f ^ 

The solution of this equation is 

a„ — A cosh ax B sinh ax 

Substitution of (387) in (388) gives 

B = alo/W 


(388) 
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If Ii be the current in the conductor, then 




Wadx 


Hence A = a[(/i/8inh ad) — /q tanh Jad]/Tf 

So we obtain 

r cosli oloc "1 

(Tj. = ad I (Tj + (To sinh cnx — Oq tanh ^ad cosh ouc I . (389) 

For the slot arrangement of Fig. 39 the current densities 
are given in the following table, where tr^ = cti = a„. 


Level 

Density 

1 

Of^oid [coth ced 4- tanh Jati] 

2 

CnCcd coshec Jarf] 

3 

OnCed [coshec oed — tanh Jacf] 

4 

a^oid [coth oid] 

5 

a^ad [ J coshec Jad] 

6 

a^OLd [coshec ad] 


-(I) 

-( 2 ) 


-(3) 

(4) 

-(5) 

~( 6 ) 


Fig. 39 


8, Voltage Drop in a Conductor per 
Centimetre Length. The total voltage 
drop in a conductor is due to both sets of fluxes. The 
“ efiective ” flux will set up a voltage, which is the same in 
all sections, so we use the section where x — d and get 

pa 

V = pa^ == ■j[|F(A cud + /q tanh \oLd) per cm. . (390) 


The “ ineffective ” flux is due to the mutual action of 
currents in all conductors lying below the one under consider- 
ation plus the self-induction due to the current in the con- 
ductor itself, or 




47rd 


A + 



{Wadx)dx 


477 

as 


jet) 


cosh ad — 1 
\ ^ sinh ad 

— /q [tanh |ad(cosh oed — 1) — sinh ad] + (x,d{Ib — Ig) 
• ^ [(iA + /o)2ad tanh ^ad + aM^C/j - Jo)] 
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So the induced voltage is given by 

e = ^ [(R + /o)2ad tanh - h)] • . (391 ) 

A proper combination of V and e 
gives the voltage drop in the con- 
ductor. 

For a full-pitch winding with 
solid conductors (Fig. 40), and if 
of the 2n conductors in each layer 
we consider the gth conductor 
p. from the bottom of the slot, then, 
since there are (<?—!) conductors 
below the conductor, the effective 
current is Jq = 1) = A- So 

that we get 



V = ^ [/i coth ad -f /i (5 — 1 ) tanh ^ad] 

(g — 1 )e = ~ [(i + g — 1 )2ad tanh iad](g — 1 ) 

A P ^ 

ifd -M- 


(392) 


So we obtain the total impedance as, 

2n 2n 

1 1 

where M = ad coth ad = Jf , -f jiff*, 

N = 2ad tanh ^ad = Nj. + jN„. 

So we have that 

For further development, reference should be mae to the 
papers listed in the Bibliography. 




EDDY CURRENT EFFECTS 


137 


REFERENCES 

Latour; Proc, I.E.E. (1919), Vol. 7, p. 61. 

A. B. Field: “Eddy Currents in Large Slot- wound Conductors,” 
Trans. A.I.E.E. (1906), Vol. 24, p. 761. 

: “Idle Currents,” Journ. I.E.E. (1906), Vol. 37, p. 83. 

L. Fleischmann : Arch, fur Elektrotechnik (1914), Vol. 11, p. 887. 

H. W. Taylor: “Eddy Currents in Stator Windings,” Joum. I.E.E. 
(1920), Vol. 68, p. 279. 

W. V. Lyon: “Heat Losses in the Conductors of A.C. Machines,” 
Trans. A.I.E.E. (1921), Vm. 40, p. 1301. 

: “Heat Losses in Stranded Armature Conductors,” Trans. 

A.I.E.E. (1922), Vol. 41, p. 199. 

A. Press: “Skin Effect b Slot- wound Conductors,” Electrician (1924), 
Vol. 92, p. 4. 

Gilman: “Eddy Curr^ tt Losses in Armature Conductors,” Trans. 
A.I.E.E. (1924), Vol. p. 246. 

Lyon, Wayne, and H derson: “Heat Losses in D.C. Armature 
Conductors,” Trans. E.E. (1928), Vol. 47, p. 689. 

F. W. Carter: “Pa^ . Current Losses in Bar-wound D.C. Arma- 
tures,” Journ. I.E.E. (iJ26), Vol. 66, p. 1002. 


lo— (T.152) 



CHAPTER X 

FOURIER’S SERIES AND INTEGRALS 


Heaviside frequently points out the parallel between his 
methods and those of Fourier. We now use the latter to explain 
some points of the expansion theoiem. As the mathematical 
proofs can be found in standard ma'^^hematical texts, we will 
merely sketch the theory. Solutions ;f differential equations, 
where the disturbing force is expresse 1 as a series, are dealt 
with and exhibit no novel features. 

The integral is a natural developm ;nt of the series — in 
electrical language it gives the freqi icy spectrum of the 
array of impulses denoted by the seri When the impressed 
force is given by the integral the solu is readily obtained, 
provided that the resultant integral can be evaluated. An 
explanation of the unit function can be derived from this 
theory, and the function is seen to consist of terms of all 
frequencies. A very interesting use of the integral is in an 
integral equation, which in one form is known as Carson’s 
equation, and in another as jit'oinwich’s equation. From these 
equations the equivalencies oi a great many o oerators have 
been evaluated by Campbell (see Bibliography). 

1. General Theory o! Periodic Functions. If f{i, is a single- 
valued function, then in the interval — s < ^ < 6 it may be 
represented by 

■nt 27 r< 

/(O — 2-^0 “h -^1 T -^2 ~ h • • • 

o o 

I R • . R • , 

+ Bi sin 1- jfjg sin h • • • 

s s 

With the proviso that /(<) has 

(i) a limited number of discontinuities, 

(ii) a finite number of maxima and minima, and 


(iii) that 

j* exists. 


We have 

^ J fi^) cos ^ dX . 

. (397) 

and 

j J f{^) sin ^ dX 

. (398) 
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The series may be written in several forms, for example 
with 77 /s = ft> 

(i) If/(<) is an even function, i.e. f{t) = + /(— 0> fhen 

f{t) = ~ f cos f f{X) coBncoXd{coX) (399) 

^ Jo J-2)r 

1 

(ii) If/(0 is an odd function, i.e. f(t) = — /(— t), then 

00 ^ ^ 

f{t) = —y sin ncot /(A) sin ncoX . d{coX) . . (400) 

J-2(r 

1 

There are other forms of these expressions, but all converge 
to f{t) in the interval — s < t < s, and at points of discon- 
tinuity become equal to 

Lt i[/(< + u)—fit-u)] 

u-*0 


2. Determination of the Wave Form of a Function. This is 
the operation inverse to that with which engineers are familiar. 
Here the series is given, and it is required to find the wave 
shape. If 


2 'ST 1 
1 

or operationally 


. rnr 
sin — t 


. (401) 


riTT 

to — ,rt CO 

/•M - - V i ^ ^ _ i- V 2 

^ TT / f n / mrV ~~ sp / . ( nu \^’ 

V p= + (-) 


then on expansion 


m 


sp 


1 + 


\spj \sp/ 


+ • . • 


. (402) 


But we have the following partial fraction development from 
the theory of complex functions — 


sp coth sp = 1 -j- 




. (403) 
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Equating (402) and (403) we get 

/(<) = -^ (fii’cothsi?- 1) .... (404) 

bp 

Integrating the last term, and expanding the ooth terra, we 
have 

f{t) = 1 - - -I- 2e-2«> + 2e-*«’ (406) 

8 

Taylor’s theorem gives 

e-2»p/(<) =f[t- 2s) 

So, in equation (405), when t < 2s, all exponential terms 
are zero, and the value of f(t) is (1 — t/s), i.e. a straight line. 
For the period t — 2s to t = 4s the first exponential term is 
present, and we get 

fit) = [1 - t/s -f- 2f(t - 2s)] = [3 - </s] t > 2s 


We are now able to draw out the wave shape as shown in 
Fig. 41. 




(406) 
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Now, on displacing the time axis downward, we get a wave 
formulated by 

, ^ E 2E sin ncot 


(Fig. 43) . (407) 


“ 5 | 


n 




Fig. 43 

3. Application of Fourier’s Series to Circuits. Case l. Let 


a voltage represented hy E = - Eq applied to 


an 


RL circuit. Then we get 
sin ncot 




n{R + Lp) 


. (i) 


On rationalizing and discarding imaginaries, we get 

2 jj, 1 pi2 cos ncot + ncoL sin ncot g-(K/r-)» "j 

* "" n^^Z^nl R^ + n^(o^L^ + n^o>^L^\ 

= ^ ~ ~ • ( 408 ) 

where tan<^„ = and Z„ = y/{R^ + n^co^L^). 

00 

Case 2. With a voltage E — Eq-\- y E„ cos (ncot + $„) 
applied to an RLC circuit, we have 




[LCp^ + BGp -f l]5f = C[Eq + En 008 {ncot -f- f^)] 

1 

Or, on substituting a = RI2L and coq® = l/LC, we get 


(p® + 2ap + coo^)q = ^ E„ cos {ncot + ^n)VL 

1 
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With oscillatory conditions, the first term is 

pm 


+ ^ ( _ « _ a _ _y^)J (409) 


This yields, after some simplification 

e““* cos {^t — (^) 


— ^0 2 ■ 
L‘"o 


] 


(410) 


/SV(«^ + /S") 

where /? = V (coq^ — a®) and tan (f> — a//S. 

For the second term, omitting the summation sign, we get 
after rationalizing and discarding imaginaries 

I cos (ncot + — <^„) 

^2 - j^[(2anw)2 + ((oo^-n^co^)^] 

e-“‘r a sin (i„ + /?)«— (nco — ^) cos (f„ + /?)« 


+ 


£““* r a sir 

tL“" 

a sin — 


a® + (wo) — /9)® 

(ntu + |8) cos /5)< 


+ (ria» + /9)^ 

The full expression is given by 

Enf I 1 


J 


(411) 


, ^ En I COS (ncot + ^« — ^n) 

Zw ^ lV[(2ana;)2 + (V-to"r^2)2] 


n=l 


e-“Tcos (Pn<+ l«+ <9/) cos i„ + ^n") ']) , 

/S« L V[«^ + - iSn)? V[«^ + (no + m Jl ^ ^ 

, and tan 0„" 


where tan 
a 


cmoi 


(Of^ — n^(o^’ 


tan d„' = 


noi 


|8 


71(0 + /S' 

The steady-state term is given by 

n ^ 0^1 cos (^^ft;< -f |„ — .^„) 

^ L(o^^^ Z^L ^/[^2(x.no>Y -f (two^- nW)®] ' ^ ' 

1 

and on differentiation we get the steady -state current as 

00 

noiE^ sin {n(ot -f — ^„) 

L V[(2awa))2-1- (oo®-w®ft)®)®] 


. ( 414 ) 
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Case 3. Response op an Oscillograph to a Damped 
Sine Wave. The force equation for an oscillograph string is 


(F 


00 

+ + o)^^)y = ^0 + 2 cos noit . (1) 


The solution of this equation is deduced readily from the 
preceding as 


CD 

F„ e~^»* C 08 (ncot — ^„) 


+ 


“)]* “h n^(o^(2/u — a)®} 

1 

+ two other transient terms 
nco(2y—oc) 

where tan 

4. Fourier’s Integral. The series may be written as 


(416) 


2 1 . nvt r* . /nnA\,. 


or, discarding the first term 




/ njrt \ ( UttX 

cos(^— jcos(^— 


) 


+ sin(^^j8in(^)JdA . (416) 


As cos ( — 0) = cos 0, we get 
00 


(417) 


Now let the period become infinitely large. With m = nnje 
the step Trjs wifi become dm in the limit, or 

j ^00 ^00 

/(O = ^ I dm \ f{X) cos m{t — X)dX . . . (418) 
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Putting m == 0 ), and noting coa cot = we get 


These integrals are subject to the same limitations as the 
series. For an even function we get 


2 f® r® 

/(<) = - I cosm<dm /(A) cos mA . dw. 

Jo Jo 

. (420) 

If we write g{t) = mtf{m) dm 

. (421) 

then 


f{t) — j cos mtg{m) dm 

. (422) 

These formulae are known as the Fourier cosine 
The sine transforms are given by 

transforms. 

h(t) = sin mtf(m) dm 

. (423) 

f{t) = j sinmth{m)dm 

. (424) 


5. Conversion o! the Series to the Integral. As an example 
of the conversion we consider the series 


Z“ 


2E(^ sin n{'!Tls)t 


{n odd) 


(426) 


Now let the period become infinitely large. By the method of 
last section we get m = mrjs = co and dm = n/s — dco. 


So 


^3 = 

3 

77 

2Eq dm 

TTU 

TT 

ms 

TT ' m 

11 

« r 

E^ sin cot _ 



dft> 

Jo 

77 

CO 


2Eq dco 

TT CO 


(426) 


On plotting the amplitude of the component harmonics 
against to for equations (426) and (426), in the case of (426) 
we get a row of impulses at definite values of co (Fig. 44a), 
while for (426) we get a smooth graph or continuous frequency 
spectrum (Fig. 446). The contribution or response made by 
frequencies between co and {co dco) is given by the shaded 
area. 
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For a square wave given by 



Cb) 


Fiq. 44 

when the period becomes infinitely large the wave is given by 
the integral 

T, ^ fi 1 f" sin _ “1 

Carrying out the same process graphically, we get Fig. 45 (a). 
This is seen to correspond to the unit function of Chapter I. 
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So the admittance is either 

B sin (ot 


2 r 

^Jo 


d(o 


or 


2 _ 2f 

B ttJo 


B^-hX^ CO 

X cos cot 


d(o 


(437) 

(438) 


(ii) For a condenser circuit we get 

1 B . . l/coC 


Z~ B^+ 1/0)2(72 
So for unit voltage we get 

B sin cot 


+ j 


B^ -f 1/0)2(72 


A(t) 


2 r* 

Wo 


Z2 


d(o 


CO 


(439) 


The value of this integral is ^ so the current is 

i = {EIB)s-*l^^ 

(iii) When the applied voltage is alternating we get 

„ . r 1 p sin Xt I 
e = la sin o)on f + ~ I - - j- dX 

E . jE? f” 1 

== ^8ino)(,( + ^ -r[cos (cOq— X)t— cos {(Oq-\- X) t'\dX (440) 
" "TtJo ^ 


or e 


_ r, 1 r“ sin x.t ,,i 
E cos o)q( ^ “* I dA I 


E E pi 

= 2 coso)o(+ ^ [cos (o)(,— X)t + cos (o)o+^)<] • dX (441 ) 

These equations are zero for negative values of t, and are 
E sin cot and E cos cot respectively for positive values. On 
integrating equation (440), it may be shown to be equivalent to 


Er Ip 

2 + 4 


(On 


A2 


COS Xt dX 


(442) 


Now, suppose that this voltage is applied to a line B — G~0, 
and that the far end of the line is terminated through a resis- 
tance B. Then 

o-2oti- 


z 


“Ll 4- 
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and 


1 1 ri + ^£-2“'']* 

~ Z~ Zoll- 

= 7o[l + 2jue-^ + . .] . 


(443) 


Now, as Zq = \/ (L/C), then al — y/iLC ) . I — jco^T, say, so 
we may write 


Y = ro[l + 2/ie-2>.T + . . .]-= A(t) (444) 

When the applied volts are E sin oit, the current is given by 


i(t) = [1 4 - 2^£-2i'“*2' 4- . . .] sin (Off, 



Inspection of this equation reveals that the terms come in 
at different times ; thus in the intervals 

0 to 27 we get 

EYf, . , EY„r (oo 

sin co^t -f- I n n cos o)t d(o 

2 27r j 0 

27 to 47 we get the above plus 

EY EY 

sin Q)t 4- I 2jue~^^‘"‘^ — ^ — » cos cot do 
2 27r j 0 cOq*— or 

and so on. 

As e- 2 j<ur, pQQ _ gQg _ 2<x)qT), it is seen that the 
increment is retarded in phase by 2o)o7 and is also attenu- 
ated by 2fi. The form of the wave with respect to time is 
shown in Fig. 46. 


♦ See Section 3, Chapter VIII. 
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7. The Carson Integral Equation. Let us consider the 
differential equation 

(ao2)” + Oi^)"-i + . . .a„)t/=/(0 . . (446) 

or F{p).y=f{t) 


The formal solution of this equation is well known ; but we 
follow a method due to van der Pol.* For the initial conditions 

let p”y = y • . . — y = 0 


On multiplying both members of the equation by e~^ and 
integrating between 0 and oo , we have 

^00 f%CO 

f{t)e~^dt . (447) 

Jo Jo 

^00 

For the term p”y £"*’* dt, by partial integration, we get 
Jo 

^00 

{p»-^y)e-^dt . . (448) 


dt 


The first term of this expression is zero by terminal condi- 
tions, and so 

^00 ^00 

I {p'^y)e~^ dt = p\ (p”“^y)e“^‘i 

Jo Jo 

Repetition gives 

^00 ^oo 

I df = p" I ye~^* dt 

Jo Jo 


. (449) 


So the first member of the equation becomes 

/«00 ^00 

[«oP" + + • . • + ««] dt = H(p) ye-^* dt, say 

Jo Jo 

. (460) 

For the second member, we will let /(<) represent the unit 
function, and so 


r 1 

/(Oe-"‘d< = - 
Jo P 

r“ 1 

So we get H{p)j ye~^^dt = - 


or 


1 r" 

f77-T= ye-P*dt . . (461) 

P-H{p) jo ^ ' ' 


* Bibliography. 
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This equation was first given by Carson, and by it we are 
able to interpret operational expressions in terms of the 
Fourier integral, as will be shown. 

8. Development of Operator Equivalencies. We now sum- 
marize the methods used to develop the equivalencies of 
operators. In our work we have used several methods, and 
so we will refer the reader to these as examples. 

(i ) Expansion of operators by means of the binomial theorem, 
by Taylor’s theorem, or by partial fraction development. 
Examples will be found in Chapter I. 

(ii) Using known operators and deriving therefrom others 
by the use of the shifting theorem, Borel’s theorem, etc. (see 
Chapter I. 

(iii) Using the integral equation of the previous section. 
This is a particularly powerful method, as all known Fourier 
integrals can be given in operational form, and reference 
should be made to a list by Campbell (see Bibliography). 

We have already used method ii to evaluate operators in 
the form of a product (Chapter I). The idea has been extended 
by Carson and by van der Pol,* to which we must refer the 
reader for further information. 

Method iii has not been used so far in our work, and we now 
show how it is handled. From Section 7, we have the integral 
equation 

1 r® 

— U7~\ — 
p.H{p) Jo ^ 


Let y — t". Then by the well-known integral, we have 



dt = 


1.2.3 . 

^n+l 


n 


So 


In! 

= 3^+1 or H{p).y = p^—^=^ 1) 


pH[p) p” 


and 


1 


V' 


-1 


p'*~^ (n — 1)1 
As a further example, we have 


etc. 


I 




P 




n I 


♦ loc. cit. 
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SO that + X), which is the equivalent of 

sin Xt. In Appendix I we have collected some of the more 
important results for ready reference. 
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CHAPTER XI 

CIRCUITS WITH VARIABLE PARAMETERS 


General. The assumption of constant parameters does not in 
general agree with conditions obtaining in practice. For ex- 
ample, it is known that resistance varies with temperature and 
that the inductance of an iron-cored coil varies with the current. 
There are numerous other such effects. In subsequent dis- 
cussion we segregate parameters as follows — 

1. Those varying with respect to x, the independent variable, 
where integration is with respect to x. Such equations lead to 
solutions in a power series, of which the best known is the 
Bessel. The differential equation is of ordinary linear type 
with polynomial coefficients. 

2. Those depending on the unknown y, the dependent 
variable. Such equations lead to a functional integral equation. 
As the solutions of these are given in the form of successive 
approximations, we will discuss them and pass to graphical 
methods. 

3. Those depending on the unknown y, but the connecting 
function is complicated. Solutions can only be obtained 
graphically or by use of a machine integraph (see Bibliography, 
Chapter VI, for reference to the latter). In general it will be 
found that there is no great difficulty in setting up the equation, 
but the numerical solution of it is quite another matter. So, 
for engineering purposes at least, resort is had to graphical 
methods. These lead to no general result, as each problem 
must be treated anew. 

1. Power Series, Bessel Functions. To introduce the first 
group we consider a line whose characteristic constants are a 
function of the line length, e.g. let Rx = RJx and Cx = Cffc. 
With G = L — 0, we have 


ax X 


. (462) 
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or 


dx^ X dx 


B, 


kW=:0 . 

V = Al^ikx) + BKf^ikx) 


( 454 ) 


where k^ — — Cfpcp — A and B are integration constants. 

oc 

The solution of this equation is given by a Bessel series. 

2. A Summaxy of Mathematical Results. For an equation 
such as 

^ 1 ^ 

dx^ X dx 

where m is integral or nonintegral, it is shown in mathematical 
texts that the solution is given by 

y = AJ^ikx) + BY^{xk) (456) 




(456) 


where A and B are constants determinate from boundary 
conditions. Here Jm(kx) denotes the series of the Bessel 
function of the first kind and mth order, and Ym{kx) that of the 
second kind and mth order. For the zeroth order we have the 
following expressions — 


Jo (kx) = 


(kx)^ (kx)* (kx)^ 

■!" 2H^ ~ 2 * 4 ^ ■ 


■ (467) 


when the argument is small. When it is large, then 

Jo = y( (*® - •^/4) • 

For the second kind we get the following expressions — 

Fo {kx) = I {<r.— \gn,kx)jQ{kx)-\-^^~— (1 + i) 


(458) 


(459) 


where a = Ign* 2 — y = 0*11593 for a small argument. When 
it is large, we have 

Fo S {kx- 7114.) . . (460) 

Generating Function. The form of the expression for the 
various orders may be developed from an expression 

_ g — !t/2< _ _ ^ (461) 

On expansion we have 

+ ? + + • • •][^“l + ^l(^)) •] 



CIRCUITS WITH VARIABLE PARAMETERS 


165 


Multiplying out and putting f = 0 gives the expansion for Jo(*) 


\2) ^ • • • 


Collecting the terms in t, we get the expansion 

T / X ^ 

— ' ■ ■ 

Or we have that 

_|. tJi{x) + • ■ • (462) 

On substitution of — t, the generating function becomes 
giarsin 6 _ J^[x) -f- 2[J^ COS 26 -|- J 4 COS 40 . . .) 

+ 2j(JiBin 0 . . .) 

Expanding and equating reals and imaginaries, we have 
COS {x sin 0) = Jo{x) + 2(J2 cos 20 + J4 cos 40 . . .) , (463) 
sin (a: sin 0) = 2 {J^ sin 0 -|- 3 sin 30 . . ,) . , (464) 

And it may be shown 

1 r®" 1 

•^o(^) = T- ® d0 = :r- e^*'“*®d0 . (466) 

Ztt Jo Jo 

or generally 

A-m rir 

«^m(«)='^ e^*®“®cosm0d0 . . . (466) 

aSTT jo 


gjx 000 9 QQg ffiQ 


(466) 


Differektiation. Differentiating Jff{kx) with respect . to 
{kx) gives 

kx {ikxf _ , . 

2 *1” 2*4 • * • {kx) 

If the differentiation is carried out with respect to x, we get 


k^x 

+ 2 ^ • • • 


kJi {kx) 


So we get the Recurrence equations 

xJ'm{x) = mJM — xJm+x{x) 

= — mJ„{x) + a:J„_i(a:) 
Modified Forms. For an equation in the form 

dx^ xdx ■ 


(467) 


. (408) 
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with ^ = 1, the solution is 

y = ^Jo ijx) 4- BYq {jx) 
or, as usually written — 

y = AIq (x) + BKq (x) . . (469) 

where /q is the modified form of Jq of Fq. 

For small arguments 

^0 (^) = 1 + ^ + 2^ • • • = 2^; Jo 

and for large arguments 

gx r j 2232 j *1 

"" \/{2jtx) + si ^ • (sip ■!■•••] 

For series of the second kind, with small argument we get 

Ko {x) = + lgn|Jio (a:) 

+ ^ + (^ + i)2^ + ■ • • • (472) 

and for large arguments we have 

'"> = '■V(£)[‘ - + ^‘ ■ (li? + • ■ •] <«») 

The recurrence formula gives 

I -mix) = Imix) : h'ix) = hix) 
K-mix)^K^{x):K^'ix)=^-K^ix) . . (474) 


Oscillating Functions. For an equation such as 


dx^ 



the solution is given by 

y = AJq ijikx) BYo ij^kx) . 

The generating function is now given by 

gim-m = [1 + yixt +...][{- ^fx/t . . . ] 
So 




or, as usually written, ber x j bei x . 


(476) 

(476) 

(477) 
(47S) 
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For functions of the second kind, we have 

^0 = ker a: + j kei a: . . ( 479 ) 

where ^4 

kera:= [a — lgna:]bera:-f (7r/4)beia:— . . ( 480 ) 

kei a; = (a — Ign a;) bei x — (n/i) 

ber a: + ^ — (1 + i i ‘ ' ‘ ^ 

Differentiating, we get 
d 

JJT-. Jo = ber'(Aa;) + j bei' (kx) 


so — = k (ber' kx + j bei' kx) . , (485 

Approximations. When x > 10, we may use the following- 


ber 

X = 

l^COS 

1 

> 

tt/S) -f- 

1 . 
^sin 
8a: 

{x /\/2 + 

,/8)] 

bei 

X — 

j^sin 

(W2- 

77/8)- 

1 

— cos 
8a; 

{xl \/2 + 

„/8)] 

ber' 

X = 

l^cos 

{xJ^Z 2 4 - 

77/8)- 

3 

(W2-~ 

»/8)] 

bei' 

X = 

sin 

{xJ \/2 -f 

77/8)- 

3 . 

{x|^/ 2 - 

„/8)] 

ker 

X = 

cos 

(a;/-y/2 -j- 

77/8)- 


L {x/^ 2 -\- 

.7/8) j 

kei 

X = 

f— sin {a;/\/2 - 

L. 

f ^/8) - 

-^cos(x/^/ 2 -^ 7 tI 8 ) 

ker' 

X = 

— COS (xl\/ 2 - 

-^/8) + ^8in{a:/V2- 


kei' 

X = 

sin 

(a:/ V2 - ■ 

77/8) + 

3 

(X/V 2 - 

,/8)] 


, )/ 

V(2va:) ® 




When the arguments are very large, then the first term is a 
sufficient approximation. 

Vector Diagram. An oscillating function may be repre- 
sented by a right-angled triangle with the hypotenuse or 
modulus 

Mf) l(f>Q = (ber* x -f bei* a:) 
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at an angle <f>Q to the base line, such that 
tan (f)Q = (bei a;)/(ber x) 

The second kind may be similarly represented 

iVo /0o = V (ker* x + kei^ x) tan 6 — (kei a;)/(ker x) 

The following facts are useful in respect to differentiation — 

^ . Jo ijh) = her' x+j bei' x . . (484) 

^^••^oO‘^*)=i"*[ber'«+jbei'a:] = -Ji(j«a:) . (486) 
So we may write 

ber' a: + bei' a; = — (Jlf^gM) . . (486) 

As — we have that 

Y'(ber'* x + bei'^ x) = l(<f>i — tt/4) . (487) 

Similarly we get that 

Y' (ker'* X + kei'* x) = /{d^ — tt/4:) . (488) 

Thus we get readily that 

ber' X = Ml cos {(f>i — n/i), bei' x = sin {<f>i — 7r/4) 
and the following combinations may be developed — 
ber X . bei' a: — bei a: ber' x — sin {<f>i — ~ 

ber X ber' a: + bei a; bei' x —M^Mi cos {<f>i — <f>o— 7r/4) 

bei' X ker' x — ber' x kei' x = M^Ni sin (<f>i — 0^) (489) 

ber' X ker' x + bei' x kei' x = M^Ni cos (^^ — B^) 

Tables. Tables of Bessel functions will be found in several 

of the references to this chapter. We give below the roots of 
where m is the order of the function and 5 that of the root. 


! 

a 

1 

m = 0 

m = 1 

m = 2 

m =« 3 

m 


1 

2-404 

3-832 

6-136 

6-379 

7-686 

8-780 

2 

6-620 

7-016 

8-417 



12-339 

3 

8-664 

10-173 

11-620 


14-373 

16-700 

4 

11-792 

13-323 

14-796 

16-224 

17-616 

18-982 

5 

14-931 

16-470 


19-410 


22-220 

6 

18-072 

19-616 

21-117 

22-683 


26-431 

7 

21-212 

22-760 

24-270 

26-749 


28-628 

8 

24-363 

26-903 

27-421 



31-813 

9 

27-494 

29-047 



33-612 

34*983 
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Mental picture. A general idea of the form of these functions 
is as follows — 

(i) with real argument is of the same shape as e", 
with the exception that at a: = 0, /^(a:) = 0 ; and at x — oo, 
I^{x) = oo. 

(ii) Kf^{x) behaves as with the restriction that at a; = 0 
it is infinite, falhng to zero at a; = oo. Functions of this type 
are frequently eliminated by boundary conditions. Recalling 
the definitions of the incident and reflected waves, we see that 
lQ(a:) is the reflected and Kq{x) the incident wave. 

(iii) Jo(a:) approximates to a cosine wave. 

(iv) Fo(a:) approximates to a sine wave. 

From the table of roots it is seen that the periods of these two 
functions are both when the value of x is above 10. 

Operational Expressions. As we treat the problems of this 
chapter by formal mathematical methods, we will now show 
some operational expressions which lead to the Bessel series. 
For example, we have on expansion as a convergent series the 
following expression — ' 

A 1 la* 1-3 o* 

V(A*- a*) “ V(1 - aVA*) “ ^ + 2 A* 2* ’ A* ‘ ’ 

= 1+15 a*:®* + 2^ “*** • • • 

= /„ {ax)* 


On expanding the expression- 


' A 
A + 2« 






+ ('AY 

2a ^ 2 I • V 2a / 


^ _J L_ 

V 2a ■ 's/(2a) ’ 

•• (a + 2a) ^ v'(2^aa;)[^ j 


_i_ 

* O l/O^A«' 


‘ \A+2a/ * y'(27raa;)L Sao: 2!(8aa:) 

_ g-d* j^(aa;) , , ,-j- 

/ A /A — aV A 

A+to j = i ATi j 

expression at the start. 

♦ See equation (470). 


• ■ 


t See equation (471). 
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If now we substitute 


^2 \ 

^2 j for a^, then 


Vl>2 - aMl - xyvH^)] »V L / 

An expression such as 

2A ^ 2A 

A — a A 

9e-a* 

“ V(2aA- A®) ' A-a 
By expansion and integration, this becomes 


T 2 Hi r 1 ,12 . 32 1 ^ ^ , 

L7raa;J ‘L 8aa; 2 !(8aa:)2 • • • J“" 2 0 (®^)- 

Further expressions are given in Appendix I. 
"s. 3. Skin Effect in Isolated Tube. Revert- 

S/\ ing to the problem discussed in Chapter IX, 
^ \ we will here take into account the field set 
\ j up by eddy currents. For a tubular con- 

I y ductor as shown (Fig. 47), the voltage 

/ across the ends of an elementary tube of 
radius x and thickness dx is 


Fig. 47 


But Ri 


pi 

2‘nx dx 


a (27 tx dx) — pa per cm. length. 
d<f> 


where p is the resistivity and a the current density. The 

r* 21 

disturbing flux is ^ = j — - dx 

where /* is the current enclosed between radii a and x. So 
dj, 24 2 r* 

~ I 2iTax dx 

dx X xja 

On differentiating with respect to x, since the p.d.s across all 
tubes must be the same, we get 

da d 4aT C* 

P1Z~ Tf~:^ ' \ ox dx=z0 

dec Cbv 00 j ^ 


ax dx — 0 
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Removing the integral sign by a second differentiation, we have 

d^a 1 da da 
dx^ ' X dx~ p dt 

When a varies harmonically with time, we have 


(490) 


\ d ( da\ 47r . , 

•*» U'tC/y lA/tL J p 


(491) 


where w* == 4aim/p. 

The solution of this equation is 

a = AJq ipmx) -f BKq (j*mx) . 

Usually, this equation is written in the form of oscillating 
functions, thus — 

a — il[ber mx + j bei Tax'] + B [ker mx + j kei tux'] . (492) 

A and B are constants determinate from boundary conditions. 

Determination of A and B. For an isolated tube, at a; = a, 
then Ix = 0 and da/dx — 0. This condition gives 

.4[ber' ma -}- j bei' ma] + B [ker' ma j kei' ma] = 0 
B ber' ma + j bei' ma 
A ~ ker' ma -f j kei' ma 

r 

Also the total current in the tube is I == 2n\ xa dx. 

I = 2i 7- j* x[AJ^ + BKf\dx 

jx ber mx . dx = (x/m) bei' mx 

\x bei mx .dx — — (xjm) ber' mx 


(493) 


Now 


So I 


= 2it ^ 




A (bei' mx — j ber' mx) 


X *]^ 

-f- — -B (kei' ww; — j ker' ma;)J . (494) 


From equations (493) and (494) we get 
ml 

A — 2^ (ker' ma + j kei' ma)IA {mb) 
Tthl 

— B= (ber' ma + j bei' ma)l A {mb)j 


(496) 
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where 

A {mb) — (bei' mb— j ber' mb) (ker' rm -f j kei' ma) 

— (kei' mb— j ker' mb) (ber' ma j bei' ma). 

For any numerical problem the values of A and B may be 
evaluated with the help of tables. 

Resistance and Inductance. When the current is uni- 
formly distributed, the resistance per centimetre length will be 
/t)/77-(6®— a*). At a: = 6 the disturbing flux is zero, and the 
voltage drop is given by Due to the voltage drops being 
equal, we have 

ZoCTo = = IZ' (say ) 

So = (496) 

Substituting, we get 

I -4 (ber mb -f j bei mb) + B (ker mb -f j kei mb) | 
” 2n |[Aa; (bei' mx — j ber' mx)-j- Bx (kei' mx—j ker' ma:)]*/ 

Let X — B/A. 

Then we have 

_ •P*” r (ber mb -f j bei mb) + % (ker mb -\-j kei mfe)! 

~ ^ 27rb L (ber' mb -\-j bei' mb) + x kei' mb ) J 

or 

Z' jmb 6* — a® 

bZ^T~W 

t (ber mb + j bei mb) -j- x (ker mb -\-j kei mb) 1 
(ber' mb -|- j bei' mb) x (ker' mb -f j kei' mb ) J ’ ' 

The real part of this expression gives the resistance, and the 
imaginary part the inductance. For further development we 
note from tables that, when the argument is in the neighbour- 
hood of 8, the “ker” and “kei” terms are small and may be 
neglected. So the expression for the impedance becomes 


Z' _ 

.mb 

r, 

1 Jo .mb r a^~\ Mo i , 



^~2 



B 

mb/ 


M 


^DO 


pj 

cos (^0 — <^1 + 37r/4) 

. (498) 

and 

X 

mb/ 

an 

Mo 


Rdo 



j^8in(^o-<^i -f- 3rr/4) 

. (499) 
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And for very large values of mb we get, 
R mb / a®\ 

X mb 


R. 


2V2 




4. Solid Conductor. Inserting a = 0, we see that B must be 
zero, and the appropriate equation is 

<r* = AJq {mxji) 

If Uq is the density at a: = 6, we have 
(To = AJo {mbf) 

Jo (fnx) 




® Jo {mb) 
or in terms of the moduli, we get 

= o'o ~ 

The current density is times the density at the surface, 

and the angle of lag is $g — Sj. The ratio ^ is known as 
the penetration coefficient. ® 

The value of the current flowing in the wire is 27 r J ax dx. 


/ = 


2rr(rn 


Jo {'•nb) 

The average density is 


rb 

I Jq (mx) xdx = 
Jo 


27r<To6 Ji{mb) 
mf Jo(m6j 


(601) 





So 


or 


= or. 


= O’, 


7r62 

2(ro / bei' m6 — j her' mJb \ 
mb \ her mb + j bei mb j 

mb Jfo {mx) , , , , , , , , „ . 

T ' Mj^{mb) 1*^0 M) + 37r/4 

w6 / her wx + j bei mx \ 

2 \ bei ' mb — j ber' mb) 


(602) 


Resistance and Inductance. As for the previous example, 
the resistance and inductance may be estimated by equating 
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the real and imaginary parts af the voltage drop. Thus in any 

d(f)x 

elementary tube the voltage drop is e* = pa^ + and for 

harmonically varying currents we get e = pa^ + At the 

outside layer the disturbing flux is zero, so 

e = p(Tq = I [i2« + 3^i) 

where R, and X, are the effective resistance and inductance. 
As I is equal to the average density times the area, we get 

D , P mb ( her mb + j hei mh\ 




'/ her 
\ bei' 


(603) 


* ' •' * ~ vft^oTat) ” Trb'^ 2 \ bei' mb — j ber' mb / ' ' ' 

On equating the reals and imaginaries in this expression, we get 

^ p mb f ber mb bei' mb — bei mb ber' mb 

■“* ~ Ia2 • 'oTI j ) • (504) 


(606) 


nb^ ‘ 2 \ ber'^ mb + bei'* mb 

p mb f ber mb ber' mb -f- bei mb bei' mb 
7r6* ■ 2 V ber'* mb -f- bei'* mb 


The latter expression for is due to internal linkages only. 
These expressions may be further simplified.* It can be shown 
that for small values of the argument 

„ . „ 1 / mb \* 1 / mb \® 




,(tJ 


/^Y 
i(2 J ■■■ 

1 

;m 262 • • -J 


(606) 


L,a-21og,(./6) 2ri 3 3 1 

Lj,o — 2 log (s/6) m6 l^ 2 8'\/2m*6* 8m*6* ‘ ‘ "J ‘ ' 

Losses in Wire. The eddy current losses in a wire are given 
by /*J?„ and on substituting we get 

TTOr^pb M-, (mb) , , . „ , . 

TT 7 W » A ' ' 1 I /CAOV 


m ilf Q (mb ) 


cos ((^0 — (^1 + 377'/4) . (608) 


5. BgssgI Cable. Heavisidef considers a cable whose constants 
are a function of the length. Thus 1/i?* = and = G^fc. 
The equation is 

— f — — ^ = m*F 

a; dx\R dx J 


AIq (mx) + BKq (mx) 


. ( 609 ) 


* Russbli.: Phil. Mag. (1909), Vol. 17, p. 624. 
t O, Heaviside: Electro ^magnetic Theory^ Vol. II, p. 212. 
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On charging an empty cable we get 
V = A/q (mx) . 


R dx 


. (510) 


— A-^If^'(mx). . (611) 


To determine the integration constant, let the voltage eXx — I 
be ^o> then 

'' = *0^ .... (M2) 


and I = ^ [/q' {mx)]A . . . (613) 

-Ko 

By the recurrence formula we get 

T _ ma; (mx) 

- R, /o (ml) ^0 

Denoting mxj^ by u for harmonic quantities, we have p — jco. 
With large arguments we have 


J„ {u) 


2 \ ( 2n + 

— ) cos u : — 

rru j \ 4 


mx f( l\ cos {Ux — S-n/4) 

■^0 V \*/ (“« “ ) 

mx // l\ cos tt, — sin „ 

-FT-./ - • ; — : . 


At X = I, /, 


- Ro^/\xJ cos Ui + sin Ui ' ® 

ml / 1 — sin 2ui N 
Rn \ cos 2m, j O' ' 


(616) 


Let Ofl = ^/{RffCo). Then m = jonoP^. On substituting, we get 
^ ^.oo^p* [ I- j sinh 2jao?p* \ ^ ^ , 




J?o \ cosh 2joL^lp^ 


(616) 


Dealing with real terms only, and noting that tanh approaches 
unity for large angles, we have 


jp _ E^a4 1 

i?„ - i?„ ■ V(^0 ■ 


• (617) 


Expansion by Operational Processes. With a voltage 
equation we have that 

_ la {mx) 

^ - h i.ml) 
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The determinantal equation is Jq {ml) — 0. Applying the 
expansion theorem, we get 



L V” ‘"■'U 

. (618) 

Let — 

— = BoCop. 


Then 



So that 

d» s 


Again 



So 

"L + w) J ■ 

. (619) 


The value of the roots of the equation Jq {slj^) = 0 must now 
be obtained. Analytically the work is complicated, but for 
numerical work the roots are readily obtained from tables. 

It is interesting to note that submarine cables are taper 
loaded. Thus the Newfoundland — Azores cable is loaded by 
mumetal wound directly on the conductors in a close con- 
tinuous wrapping. An inductance of 0-066 to 0-230 henry per 
mile is obtained. It is claimed that tapering (i) reduces the 
reflection losses from the non-loaded part, and (ii) permits the 
use of lower-permeability materials for loading. 

For further discussion reference should be made to Starr.* 
6. Filter Circuits. This problem has been discussed in 
Chapter VII, but as Bessel functions give a more compact 
solution for calculations we will now develop the theory. For 
a steady voltage E, the response is given by 

• _ sin/g,< 

nL Z-i cos STT 

1-1 

when the effect of resistance is neglected. 


• Stabb: Phil. Mag. (1934), Vol. 17, p. 83, 
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We will now write it as 

2E f ‘ cos 

0 cos stt 




dt 


(620) 


Stt 


where /S, = 2v sin ^ and the step in the summation is 7r/2n. 


STT TT , „ 

2^' = S' 


Let 0 

So with proper limits we get 

4:E r* cos {2vt sin 6) cos 2nd 


^E n f" 

" ttL Jo Jo 


dt do 


(621) 


cos' Stt 

Now, COS Stt = (— 1)*; then cos^ stt = 1. 

The Bessel generating function shows that 

cos { 2 vtsind)=jQ( 2 vt)-{- 2 [J 2 (2vf)cos26+t/4 (2«;<)cos40. . .] (622) 


All the terms vanish at the limits of integration with the 
exception of 

^itj2 

2 J 2 „ (2vt) cos^ 2n6dd = ^ . Jjn (2*^<) 

46 


/: 


So that we get 

2E 
^n- L 


n 2E r‘ sr(- 

j^J2„(2vf)dt=-^ 


(- iy.(vtf”+^* 
r(2n + s-i-l) 


dt 


2E 2n ”1“ ^ "1“ 1 
vL 2n 2s I 

2E 


'2n + 1 


{2vt). 


(623) 


= :^-G,.J 2 n + i{^vt). say. 


(624) 


For the current in the rth mesh we get 

, 2E cosh {n — r)cc sin ^,t 
”” nL / f jSg cos STT 

4E n 

= ^ J J cos {2vt sin 0) cos 2n6 cos (n — r) 20 dO dt 

With the same discussion as above we note that 

cos2n0 . cos (» — r)20 = cos 2n0(co82n0cos 2r0 + sin 2n0sin2r0) 
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But COS* 2nd = cos* stt = 1, and sin 4n6 = sin 2s7r = 0. 
So we get the effective integral as 



{2vt) cos* 2rd = ^ Jjr (2^0 


Hence 


2E f* 

if = •^2r (2^0 dt 


2E 


(vO . Cr 


. (626) 


So to find the current in any mesh all that is necessary is to 
alter the order of the Bessel function. Reference should be 
made to a paper by van der Pol,* wherein these results are 
obtained otherwise. 

7. Impedance of an Earth Return. In order to calculate the 
short-circuit fault current of a line to earth, it is necessary to 
know the inductance of the earth return. This problem may 
be solved approximately by considering all the earth current 
to return within a circular trench, the radius of the trench to 
be equal to the height of the conductor above ground. We 
consider that the flow of the earth current is parallel everywhere 
to the axis of the trench and that the magnetic lines are centred 
on that axis. 


If or is the density, then I 


rh 

I nax dx. Along the trench we 
Jo 


have pa + 


dt 


e. 


But e is independent of r, and so 


^ r _ n 
^ dr ^ dr dt 


Moreover ^ 
dr 


2 p 

r Jo 


■nra dr. 
3a 

PFr 


_ A !!i r 

r Jo 


ra dr 


or 


d^a da 3 
+ P^ = Si 


ar* 

a*cr 


ar 

1 da 


dr^ ^ r dr 


dt 

271 

P 


= — .j(oa 


when harmonic quantities are considered. 

* Van deb Pol: Journ. I.E.E. (1937), Vol. 81, p. 387. 
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Hence a — A . (mrj^) -f- BK^ (mrj^) 

At r = 00 , or = 0. So A = 0. 

This gives a = BKq (mrj^) 

Again, the total earth current 

r* 

I = I irradr 
Jo 

== ttB J* r Kq {mrj^) dr 

= W)] 

So we have 

_ I mji P Aq (mrj*) '1 
7th L^i(wAj*)J 

When we consider the voltage at r = A, then E — pa. 

So ® = 
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CHAPTER XII 


OIRCinTS WITH VARIABLE PARAMETERS (contd.) 


1. Methods of Successive Approximations. General. We have 
already seen that when the circuit parameters depend on the 
independent variable the solution may be obtained as a power 
series. We now investigate the matter further, for the ease 
where the operational impedance can be considered as con- 
sisting of two components, one which is constant with respect 
to time, viz. [Z{p)^ (i.e. t is the independent variable), and one 
which varies with time Z{p)t. Thus we have the equation 

[Z(p)o -f = E 


Uf\ - ^ 

- Z{p), Zip), 


(626) 


The first term of the right-hand side of equation (626) is 
the current due to an impedance which is invariable with 
respect to time, while the second term is due to the variable 
impedance. It is observed that the resultant current is less 
than the constant term. On writing 

EjZip), = i,{t) 


we may rewrite the equation as 


m 


1 + 


Z{p)t 


Zil>)o 

By expansion this gives 

i(<) = (1 — a -}- 

Zip)t 




where a — 
So we have 

where 


Z{P), 


{(t) — io{t) iiit) -f- *2(0 • • • 

= zlffl ’ 

*2^'^ "" zjfio ■ 


(627) 


. (628) 


. ( 629 ) 
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Now, Z{p)t . toil) is the voltage across the variable com- 
ponent with a current iQ(t) flowing. We will denote this by 
v(i) ; lfZ(p)o = A(t) is the admittance. 

From the discussion on the Duhamel integral it is seen that 
ii(i) may be written as 

d r‘ 

i^(t) = ^ J v(A) A)dA . . . (630) 

If, for example, resistance be the variable component of 
Z{p)f, then we have 

H(t) = I • £ \ri^) ■ hW] .A{t-x).d?. 


init) = £ [^(A) . /n-i(A)] .A{t- X)dX . (631) 


This method is not adapted to numerical work; but, in 
spite of the wearisome mathematics, it does give valuable 
information respecting the character of the oscillations. 

2. A Variable-capacity Circuit. To 
demonstrate the above method we will 
consider an LC circuit wherein the ff 
capacity of the condenser is varied by 
means of a cam as shown in Fig. 48. 

Then the equation is 

H^hQ+hiSinodt 

Lpi — E Fio- 48 

Here ^ = K[hf^ + \ sin cut] = Sq-\- By sin cot = -j- 

E = Lpi^[8^ + am‘l 

E-8,{t).{ilp) 

Lp + {So/p) 


So that 

or 



i{t) — iQ{t) 


s,{‘) ■ (</P) 


and 


. (632) 
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where io(<) is the current flowing when the condenser plates 
are fixed at A = iio- This equation may be written as 

m I, s.(‘> 1 , pi(«) IT I.,,, 

“ , , Ml . i “ > 2(1-) ■? + >J ■ ■ T“‘ * 

Z{p) p 

i.e. i{t) == io(0 — *i(0 "h * 2(0 • • • • • • • (633) 

where ii{t) = ^ . io(<) = ^ ^i(^) • ^ - A) . dA 


We will now develop a few terms of this approximation. For 
the first term we have 


Ep E 


where ^ 


We may consider /Si(A) . - — ~ as a voltage v{X), in order to 

bring our expressions into line with the previous discussion of 
the Duhamel theorem. Here, then — 


Vi(A) = 81 sin (wA J io(A)dA 


= Si sin coA . . (1 — cos j8A) 

Mi) = ^.sinj3< 


As .4(0) = 0, we may write 


{* v{X) .A' {t- X)dX 
Jo 

SiE r* 

~ J <wA . (1 — cos jSA) sin p{t — A) . dA 


_M r • , 

~ L 8/32 ■ Jo 
From which we get 

L_ 

hW- 22,2^2 L 


2BW + 2 / 38 ) . , 

sin (Ot - ^ 2 )(ft, 2 _ 4 ^ 2 ) 
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Further development is cumbersome, but inspection of 
results reveals that i{t) is multiperiodic. 

3. Helmholtz’s Method and Taylor’s Series. We now give 
two methods which are used frequently to give solutions in 
the form of successive approximations. The first of these is 
the method of Helmholtz, and is best illustrated by an example, 
as follows. Let 

+ h,x^ = a cos (at . . (634) 

with the boundary conditions t — 0, x = c, and px — 0. 

Let us now seek a solution in the form 

* = + c^iit) (636) 

Put a = a^c. 

On substituting in equation (534), and equating like powers 
of c in both equations, we have 

cos (Ot 

+ ^^^2 + Hi == 0 

etc 

By integration of these auxiliary equations, we get for the 
first that 

(l\ 

rf, = .4 cos nt H 5 ; cos b)t 

The second auxiliary equation becomes 
PY2 + — — hf>-^ 

But is now known from the previous solution. The repetition 
of this process gives the values of <f> 2 , ^3, etc., with 

a: = + cVa + (637) 

This method will be used to obtain a solution for a variable 
inductive circuit. 

Taylor’s Series. In radio work non-linear circuits are the 
rule, and in dealing with these recourse is had to a mathe- 
matical expansion known as Taylor’s series. This states that 
i{f{x) is a function of x such that the first n derivatives exist, 
and are finite and continuous in the interval from x — a to 
x — b, then throughout the interval f{x) may be expressed as 

fix) =f{x)o+ {x-Xo)f'ix)o+ ^? - ^pI /2(a;)o . , . 

^ (n-1)! 



r-^x)^ + B„ . (638) 



174 


TRANSIENTS IN BLBOTBIC CIRCUITS 


where x, x^ lie within the interval a, b, where all the derivatives 
are evaluated at a: = Xq, and where J2„ is a remainder, after 
n terms, which tends to zero in many (the useful) cases. 

Thus, consider a rectifying device such that its El character- 
istic is as shown in Fig. 49. Let the device operate close to 
E'q, /(,. Then, irnder steady operating conditions, a small 
variation of voltage to E will alter the current to 

I = I,-\-K[E-Eo] 


A closer approximation would be given by several terms — 

/ = /o -f K,[E - E^] -f K^[E -E^f+... (639) 

To determine the coeflEicients Ki, 
Ag . . . we differentiate and put 
E = Eq. For the first differentiation 
we have 

^ = K, + 2K,{E-B,) . . 



Fig. 49 


and so we have 


“ IdE} 


The second yields 


at Ea 


K, 


ir"! 

2! L<i«J 


at En 


which, in ejfect, is Taylor’s series. 

With an applied voltage such as 

E = Ef^-\- E^ sin o)t -f E^ sin Zoit 

we get 


/ = Jo •+• K^Ei sin o)t + E^ sin Scof] -f- K^E^^ sin (ot -f E^ sin Scof]* 
by taking two terms in the series (539), i.e. 

/ = /fl -b K^E^ sin wt -f E^ sin 3a><] -f- K^E^— E^ cos 2 tat 
+ 2EiEg cos 2o>f— 2EiEg cos 4<uf -f E^^— E^ cos 4ft)f] , . . (641 ) 


This expression shows that in this case the current response 
consists of components 
(i) due to Eq acting alone. 
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(ii) due to terms of the frequency of the applied voltage, 

(iii) due to terms whose frequency is the sum or difference 
of the applied frequencies. 

Graphical methods of obtaining the response will be seen in 
texts dealing with radio engineering problems. 

4. Some Formal Mathematical Methods. As illustrations of 
some of the mathematical methods that may be used in 
solving variable circuits we now consider simple circuits 
such as 

(i) BL circuit with R constant and with inductance varying 
as Zr = Lq(\ + e cos (ot ) ; 

(ii) RL circuit with resistance given by i? = ^ 0(1 -f ® cos cot) 
and L constant ; 

(iii) LC circuit with L constant and with capacity given by 
(1 + « cos cot). 

0 Oq 

These elementary circuits will suffice to demonstrate the 
difficulties that the solutions entail. 

(i) Variable-inductive Circuit. The fundamental equa- 
tion is now written as 


Ri -}- pLq{\ -f- e cos cot)i = E 


Denoting cot by x, we have codt = dx and with — ^ = a, 
e = E/coLq we get ^ ® 


r d 1 

a + ^ (1 + e cos a;)J z = 


(642) 


Assuming a solution, such as 

* = »o + -f • • • 


n«=i00 



««0 


(643) 


substituting it in equation (642), and equating like powers of 
e, we get 

e 

oci-o or Iq 


Dealing with the first power of e, we have 
di-i d 


Hence 
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Denoting ^ by D, we have, in operational notation — 


1 


D 


“ a ‘ D -f a ' D® + 1 


The solution for the steady-state condition is given by 

c 


a(a®+l) 


(a sin X — cos x) 


For the second power of e we have 

aij + -b D(ii cos a;) = 0 
c 


or 


(a + D)t 2 = 





a(a® + 1) 

— r 

a + D 


[a cos 2x + sin 2x] 


«D® 


7 + 


2D 


a(a®+ l) a + DVD® + 4^Da-f4, 

The steady-state solution is given by 

•• “ ir ( » » +r)(«‘ ~+ l ) • t<«* - 2) 2* - 3“ 2a:] 


In a similar fashion the values of * 3 , *4 . . . may be evaluated. 
So the full solution is given by 


e{ I , . 

^ a a2 -j- ' 1 ^ 


+ 


[(«® — 2) cos 2x — 3a sin 2a:] + 


a® + 4 

jEi r, , « • 3a . ^ 


a® + 1 


cos (Ot + 


0.“ 


a® + 4 


cos 2vDt ■ 




(644) 


Thus it is seen that an inductance of this nature would give 
rise to several frequencies in the response. 

(ii) Variable-resistance Circuit. The fundamental equa- 
tion is now 

Dpi -f i2(l -f a cos mt)i — E 

From the theory of differential equations the solution is 
obtained as follows. 
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Using the transformation w = gi, we have 


^ (1 + a cos 0)0 - ^ . (646) 


We choose g such that 

\ dg R 

gdtL^ ' 

i.e. g = 

R Ra 

where « = ^ ^ 

By substituting in (646) we have dwfdt = EgJL. 
or w — ^ jgdt — ^ j A 

where A is an integration constant. 

So i = 4- 2 e-** 

From the theory of Bessel functions we have 

00 

. « V f cos n(ot {n even)! 


-[ 


(646) 


g/>8ln<»J _ J 


mi) + 221 


j sin nwt (n odd) 


■^.m> 


where it is noted that the Bessel is independent of time. So 
that 


= wm •7 + 22/'"- {/d 

1 

f , (cos no)t\ J 
Now e“M . ,1 dt 

J (sin noit) 


cos ntat 
\ sin nu)t 


JnWW 


Ism noit 


[gitMt (_ i)n e~^^]—jn(o [£5"“* + ( — 1)" 


a* -j- n^(o^ 

When n is even, this reduces to 

, r « cos ntot 4- nco sin n(of 
2e“‘ . J- . — 

L + 7\ror 


and when n is odd, to 


^ . .fa sin no)t — no) cos mof 

_ 
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On denoting n by 2q and -f- 1, we have 

00 

. . «“* X -/>..> , fa cos + 2g(w sin 2flwn 

+ Wj) ~~~'J 

1 

. O.- V r „Jasin(2^+l)a)i-(2g+l)ft)COs(2g+l)ft)q 

+ 2j 2, L «*+(2« +!)><.> J 


3=1 

So the value of i becomes 


^ , -M. « I y 008 2 ,»( + 2 JO, sin 2 j<^) 


U-^aH^g' 


3»» 00 


2 *^2a+l iMi) 

+ (2g+l)2(«a • (2g + 1 ) ft)< — (2g + 1 )c!) cos (2g + 1 )(o<] 

(647) 


9=1 




From the expansion of the Bessel functions it will be seen that 
all terms are real, whether n is even or odd. The introduction 
of the variable resistance brings in a whole array of frequencies. 
From (i) and (ii) it is seen that the circuit deviates considerably 
from the invariable state. 

(iii) Varying Capacity. We have already used this circuit 
as an illustration of the method of approximations. We revert 
to it here as giving us a form of Mathieu’s equation. Here we 
have 

Lpi^r^ = B 

But ^ ^ (1 + « cos coi) = (0i suy 

where cc is the angular velocity of the cam shown in Fig. 48. 
Thus we have 

LpH + i[So-\- S {t)] — 0 . . (648) 

This is a form of Mathieu’s equation of which the expUcit 
solution is unknown. However, from an electrical point of view 
we are able to discuss the circuit, and so the equation. If the 
variable term is absent, 8 (t) = 0, and it is known that the 
response is periodic with coq^ = IJLC. When the variable 
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portion of the capacity is inserted, then we have ripples super- 
imposed on the original periodic response. This combination 
gives a “beat” phenomenon. Let the period of recurrence be 
T, and suppose that at ^ = 0, and pq — If,. For the 

next period we may assume that q = mQ^ and pq = ttiIq, where 
m is some multiplier. If m is less than unity, the oscillations 
are decreasing and the circuit is stable ; but if m is greater than 
unity, the oscillations are increasing and the circuit is unstable. 

If we imagine the cam to give a rectangular ripple so that 
/{?(<)= 1, then with cdq^ = l/L(7o we may rewrite the 

equation as / A.K\ 

-f ji = 0 . . (649) 


This will hold for an interval, say 0 < (ot < n. 
During the next interval n < (ot < 2n we have 


pH + 




These equations may be integrated as 

ii — Cl sin qit -f- cos qit 

where + ^) 

and sin q^^t + cos gji 

, AK\ 

L ' 


where 




. (660) 
. (661) 

. (662) 


At cot = TT these solutions should fit as regards both amplitude 
and slope; at the end of a complete period they must be m 
times as large as at the beginning. So 


(H) cofosfr 

^ (H) (Ot’^ir 


di^ 

dt 

dt 

(^' 2 ) a>e-2iT 

= m (ii) 

di^ 

di^ 

d't a)t‘=27r 

= m-T- 

<u<-=0 


On writing these equations in full we have, with w/o) = d, that 
Cl sin qid -j- Cg cos g^S — Cg sin gg^ — cos q^d = 0 
Ciqi cos qid — C^qi sin g^^ — C^q^ cos q^d -1“ C^q^ sin q^O = 0 
0 -f- mC/g — Gg sin 2q^d — C^ cos 2gg0 = 0 
— mGiqi . + 0 . — C'ggg cos 2ga0 + C^q 2 sin 2ga6 = 0 
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These equations will give the solutions of the G, if the deter- 
minant formed by omitting the C, from the above is zero. 
This gives 


w®— 2 m ^ cosg^i 0 . cosg' 26 ' 




8ing'i0sing20 1 + 1 = 0 . (663) 


Thus m = + ± 1 ), where + is the term within the 

brackets. 

When + > 1 , then m > 1 , so that the system is unstable. 

If A = 1, the system is seen to be stable. 

When + is between + 1 and — 1, m is a complex number 
and does not agree with the original assumption. The real 
part is less than unity, so that a stable system could be 
expected. 


With A less than — 1 , one of the values of m will be smaller 
than — 1 , i.e. the amplitude is reversed. The other will be 
greater than — 1 . So, after consecutive periods, we have 


1 st 

2 nd 

3rd 

4th . 

A ' 

. . periods 

m 

m* 

m® 

m* 

. . amplitude multiplier 

+ 

— + 

. . sign of amplitude 


It is seen that the amplitude of the third period will be in the 
same direction as the first, but will be larger. The criterion for 
stability is thus that 


cos q^d . cos q^O — 


+ 92^ 
2 ? 1?2 


sin q^d . sin q^d 


(664) 


is numerically less than unity.* 

6 . Linear Integral Types of Equation. The form of the 
equation given in Section 1 is known as the linear integral 
type or Volterra type. To introduce equations of these types 
we may consider an equation such as 


y„ + ai{x) . Vn-i + . y „_2 . . . + a„{x) . y == ^{x) ( 666 ) 


Let = 



say. 


* For a fuller discussion of this problem, see Van deb Pol and Strutt : 
Phil. Mag. (1928), Vol. 6, p. 18. 
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Then 

yn 


_i = ju(x) dx + Cl 

Vn-i = JJ 

j ■ ■■ J ’dx-^Gi 


n timee 


+ ('i 


yfl — 2 


(»- 2 )! 


■ ■ -h G„ 


On substituting in equation (656), we get 

u{x) + a-^{x) J u(x) dx + (ii{x ) . j* J u{x) dx ... — <f>{x) + SC'(x) 

= f{x), say 

where SC? (a:) = a^C?! + + • • • 

fx r^(x ty*~^ 

Now noting that J m(<) . if" = J — lyj- u{t) . dt, we have 

«(«) + (%(*) + a2(«) . {x—t)+a^(x) . • . . . y{t)dt = f(x) 


or 


u{x) = f{x) — I [%(») . . .] u{t) dt 
Jo 

= /(*) ~ J K{Xit) . u(t) . dt 


where K{x^) = [a^^ix) a^ix) . {x — t) . . .]; an equation 
which is of the same form as given under Section 1. 

The specific forms of this type of equation are named 
after the mathematicians Abel, Poisson, Fredholm, and 
Volterra. Equations of this class are solved usually by methods 
of successive substitutions on the lines of Caque’s method of 
Chapter I. As the question of convergence arises, a full dis- 
cussion is beyond our scope. For numerical work these solu- 
tions are practically useless in view of the slow convergency, 
and so step-by-step methods as discussed by Whittaker and 
Robinson* are used. There is one development of interest 
based on operational methods due to Koizumi, f 

* Whittaker and Robinson: Calculus of Obaervationa. Loc. cit. 
t Koizumi: Phil. Mag. (1931), Vol. 11, p. 432. 
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6. Solutions of Equations having a Kernel K (x — ^) by 
Operational Methods. The Carson integral equation (Chapter 
X (7)) is given by 


pHif) 

or in operational notation 


1 p 

W ~ Jo 


h,{t ) . dt 


H(jp) ^ 

If H(p) can be factorized in the form Hi{p) . H^{p), and 
if the auxiliary operational equations are 

WJp) ^ H^) ^ 

then, by Borel’s theorem — 


)o ^ pH{p) ~ pH,{p) . pH,[p) 


= j^[^j\w.hAt-^)d}^y^dt 

On dealing with an equation of the Poisson type, such as 

f(x) = <l>{x) - a .K(x-X)dX . (666) 

we multiply throughout by and integrate with respect to x 
between 0 and oo, obtaining 

j* f{x).e~^^dx = j <f>{x).s~^^.dx—(x.j |^J ^(A),X(a;— 

(667) 

Operationally we have 

= = ^(D), say 

5^(5) “ = «(^) 
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We have the following operational equation — 




. (668) 

So 

A(D) 

^ ~ B{D) 

1-a 


Now let 

n(D\ 

~ R(D) , 

“ D 


Then 

^(i))=-C(C)[«^-l] 




. (669) 

This may 

be converted to 



— « f f{^) -L (x — 

Jo 

. (660) 

where L{x) = R{D). 

Thus, suppose we have an equation such as 



f(x) = (^(x) ~ “ J (* ~ ^) • d}. 

. (661) 




184 


TRANSIENTS IN BLBOTRIO CIRCUITS 


When f{x) is known, the value of <f>{x) is obtained provided 
the integral is tractable. With an equation such as 


^(0=1+ f 

Jo 

K{t) = -t = -^=B{D) 




R{p) 


P 


P^+l 


== sin X 


So 


<f>{t) = 1 J . sin (a: — i)dS = cos t 


Again, if<f>{x) = J*^ P < 

Then, by Campbell and Foster’s Table (No. 601), we have 

1 


(663) 




T 


,ot— n~l 


r(a— w) 

So that = r(a— n) . 

With « = 1 and n — p, this reduces to 

x-" = r(l — p)Di‘ = B{P) 

KeveAiD)-^- ^ 

> ~ B{D) ~~ r(i - p ) . z)"-! ~ r(i - p) T{p) 


D 


But 


1 


sin ptT 

r(i — p) . T{p) ~ p 

sin /ITT 


m = 


, and we now get 
dX 




(a; - A)^-'' 


7. Graphical Methods. Possibly for practical work, where 
the solution of a definite problem is required, graphical methods 
of integration are most used. From the results obtained no 
generalization can be made, and every new set of conditions 
requires a fresh start on the problem. We now illustrate the 
methods by considering such simple cases as the rise of voltage 
and of current in an iron-cored coil, in the field circuit of a 
dynamo, and in a circuit now known as the ferro-resonant 
circuit. 
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Rise of Voltage and Current in a Coil. We consider next 
how the current in a coil builds up. The final value of the 
current at steady state will be VjR, where V is the applied 
voltage and R is the resistance. On neglecting leakage flux, 
we have that 

F = + . . (664) 


where N is the number of turns on the coil. 
Now, V = RI ; so we have 

We mav write this as 

d<f> 

R. 

at 

Solving this equation, we get 

r1.m 


( 666 ) 

(666) 


where is the initial flux at f = 0, while <f>i, is the flux at any 
given time. 

We now introduce a quantity A(f>, which is a linear function of 
Ai, so that 


A(f> 4’> 
M~ I 


. (667) 


where <f>, and I are the steady-state values. 
Using these values, we have 


t 


_ f"’ 


A<f> 


(668) 


Now, N<f>, 
at 

Then 


L,I, where L, is the self-inductance of the coil 

L, 


t 


L, f" 

R' 1. 


A<f> 


(669) 


This equation may be solved by a graphical method due to 
Riidenberg.* 

For any coil we may obtain the flux current relationship as 
shown in Fig. 50, where are represented the quantities involved 


♦ Rudenbbro : Elektrische Schaltvorange (Springer), 
13 — (T.152) 20 pp. 
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in the differential equation. We can now plot {<f>, 1/A^) direct. 
The value of the integral from = 0 to ^ is given by the 
area ABCD. On denoting L^jR by T, the circuit time constant, 
we now have that the area ABCD — tjT. Thus we know the 
value of tJT necessary to give a flux (f)^. Repetition of this 
process gives the curve {<f>, tjT), and by reference to the magnet- 
izing curve we are able to obtain the curve (i, tjT). 

Rise oe Current in the Field Circuit oe a Dynamo. 
If <f> be the air gap flux, then the generated voltage is 


e = <j>Z (W-/60) (p/a) X 10~® = K(f>, say. 



The field flux is vcf), where (v — 1) represents a leakage factor, 
and so the voltage induced in the field windings by this flux 
is given by 


vN 


# 

dt~ K ‘ dt 


(670) 


so that the equation for a self-excited machine is 

vN de „ . 

^ = e 


or 


K 
vN de 
Tdi 


= e- 


Ri = Ae 


(671) 


where Ae represents the difference between the generated 
volts and the resistance drop in the field circuit. The solution 
of the equation is 

The time constant of the winding is given by 

vN<l>g vN 


LJR 


(672) 
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SO that the simplest form of this equation is 

tJT = \\\|^e)de . . . (673) 

Jeg 

The magnetization curve (e, i) may be plotted from test 
results, and is as shown in Fig. 51 (a). All quantities in the 
above discussion are shown on this figure. The plot of (1/Ae) 
against e is shown in Fig. 51 (6). The area ABCD gives the 
right-hand side of equation (573). Thus we know the value of 



tjT necessary for the voltage to rise to e^. Repetition gives the 
characteristic (e, tJT), and the (*, tjT) function may be obtained 
by reference to the magnetization curve. 

The Febeo-besonant Cibouit: Steady State. When an 
iron-cored coil is in series with a condenser, the voltage across 
the true inductance is given by v{t) = L{i) . i{t). The form of 
this curve would be that of the magnetizing curve. Thus the 
expression for v{t) would introduce insuperable difficulties if 
we were to proceed as outlined in Section 1. We will suppose 
that the iron-cored coil and the condenser are both available, 
and that test results connecting V and i are also available. 
To introduce the method we consider initially an LG circuit 
having constant parameters. The equation is 

I/pi -+- {llpC)i — E sin (ot 

Here the potency of p is jw, and so Ifp is l/jo. Since the 
equation must hold for all values of time, we may rewrite it as 

j{o}LI — IJcoC) = Eq . . (574) 

or APPLIED VOLTS = INDUCTIVE DBOP — INCAPACITY DBOP 

Here, the relationship between the inductive drop and the 
current is a straight line through the origin at an angle <f) to 
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the current axis, so that tan (f> = cdL, as shown in Fig. 62. 
Now we select a point A on the voltage axis so that OA = Eq, 
and draw a Une AB at an angle of 6 degrees to the horizontal 
such that tan 0 = edl — IJoiC. It is seen now that OB^ 
satisfies the equation Vj,— Va = E. So, for an appHed voltage 



OA, the current will be given by OB^ Repetition of this 
construction will give the relationship between the current and 
the applied voltage. 

When the frequency is varied, it is evident that (f> wiU 
increase, while 0 will decrease, with increase of frequency. 



Fig. 53 


So that at some frequency the fines AB and OB will be parallel 
with each other, viz. when w* = IJLC. 

If we plot the relationship (i, to) for different values of the 
applied voltage, a set of curves such as those shown in Fig. 63 
result for different values of E^. Here values of E^ are shown, 
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to some arbitrary scale, as 1, 2, 3 ... . The curve marked 
0 is the natural frequency curve. The current/voltage relation- 
ship at constant frequency is obtained if we select any fre- 
quency, such as xy, and plot the current against voltage. Then, 
should the voltage vary with time, the time and voltage are 
correlated and the current time relationship so obtained. 



Now, when we consider the iron-cored coil, we get the 
saturation characteristic of Fig. 54, instead of the straight line 
OB of Fig. 62. It will be observed that the capacity fine, in- 
serted at the correct applied voltage, cuts the inductor charac- 
teristic at the three points C, B, D. In the same way as for 



Fio. 55 


the linear circuit, the current/frequency curve can be obtained 
for different values of the applied voltage (Fig. 55). The 
curve 0 here represents the natural frequency curve. It will be 
observed that in Fig. 54 the locus of the point D gives the curve 
to the left of (y = cOg, while B gives the lower and C gives 
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the upper part of the right-hand half of the curve. From Fig. 66 
it will be observed that a decreasing voltage at constant 
frequency gives an increased current, which is contrary to 
established principles. Thus the upper part of the right-hand 
curve is unstable. All the others are stable. When the voltage 
apphed to this circuit is gradually increased, it will be observed 
that, at a certain voltage, the current suddenly increases from 
a very low value corresponding to the point B of Fig. 64 to 
a high value corresponding to D when partial resonance condi- 
tions exist. 

8. Empirical Formulae. When an empirical formula is 
available, the circuit equation may be solved mathematically. 
However, any empirical formula, at best, is an approximation, 
and this fact must not be forgotten. 

Building-up of Voltage op Shunt Dynamo. To illustrate 
this method we will consider how the voltage builds up on a 
shunt dynamo, running at constant speed. We consider the 
field to have a constant inductance L, and we neglect the 
leakage fiux. The e.m.f. across the brushes is 

E = Ri pLi — Ri np^ X 10“^ . (676) 


where the fiux is expressed in megalines and n is the number 
of turns in the shunt field. Froelich’s equation may be written 
as i — (f>l{a — b<f>), and the voltage across the brushes is given 
by F = <}>M or 


<I>M 


R<f> 


a 


b<f> 


d4 


-2 


SO that 


100 a-b<f> 

n <f>[{aM-R)-bM<f>'\^ 

\ a bR 1 

~\j>{aM-R) + {aM-R) {Ma- R-bM4>)\^ 


Integrating both sides, we have that 
~nr ^ \_aM-R ' M(aM-R) ' 


At t = 0let<f> — <f>n, the residual flux. Making this substitu- 
tion, we have 

loot a ^ R /aM—R—bM<f> \ 

IT "" aM-R • M{aM-R) ' ^^^*\aM-R-bM<^ J 
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Denoting M<l> by E and M<f>g by Eg, we get 
lOOf a E B faM—R—E N 

IT "" aM- R ' M{aM~R) ' E, ) ‘ 


(676) 


This expression gives the relationship between E and t. 

9. Maintenance of Oscillations. In radio work, oscillations 
are maintained in a tuned BLC circuit by the addition of a 
thermionic valve. The explanation that the valve adds a 
negative resistance to the circuit means that the total circuit 
resistance must be zero precisely. However, such a delicate 
condition is not required in practice, and so we conclude that 
the linear theory is inadequate. Rayleigh met the same 
difficulty.* For a linear RLC circuit we have 

p^x + apx -f- = 0 


and for a non-linear circuit we get an equation, such as 
p^x -f- cof(x) . px -j- co^x = 0 


If y = px, then 


pH = 


dy dx . dy 
dx dt ^ dx 


dy 

So we get ^ ^ ^/(*) • y + ^ 


Let Q = — “(“ 

ft) 

The equation then becomes 

dqjx = dxl[F{x) — q] 


F{x), where F{x) 


-x 

f{x)dx. 

Jo 


Poincardf and LienardJ have studied this form of equation, 
and it seems, from their work, that the explicit solution will 
probably remain unknown. 

Before considering the graphical method we will consider 
an equation of this type from a mathematical point of view. 
To illustrate, we will consider a simple oscillator circuit, which 
in its simplest form can be represented as in Fig. 66. This 

♦ Rayleigh : Theory of Sound, Vol. I, Section 68a. 

t PonsrcAB^J: J. des Maths (1882), Vol. 8, p. 251. 

X Lienabd: R.Q.E, (1928), Vol. 23, pp. 901 cuid 946. 
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treatment is based on a paper by Van der Pol.* It wUl be 
assumed that the characteristic of the unstable circuit element 
is defined by a single-valued function of the independent 
variable, and that this characteristic is a symmetrical cubic 



Fig. 56 


function. So, for the circuit, we have the following relation- 
ships — 

e = Ri^ = Lpi^ = (l/p^)*3 
H = -y{e- ^e*) 

Also i-^ -j- = 0 . . (677) 

Hence the general equation is given as 

— y{e— Oe ^ jedt = 0 . (678) 

On differentiating again, we have 

Cp^e y{l - Zpe^)]pe eJL = 0 . (679) 

Using X — tj y {LG) = (o^t, a = {y — 0) ^/{LJG), and y 
= e\/[Z^yl{y — <?)], we obtain that 

+ = 0 . (680) 

If <r = 0, then y is sinusoidal, with a period = 2ttJ(Oq. As a 
increases from 0, the oscillations deviate from the simple 
harmonic type. If a is negative, the oscillations will be posi- 
tively damped ; but if positive, the amplitude of the oscillations 
will increase. 

Let a and \/[3/5y/(y — (r)] be small, so that we may assume 
that the oscillations are approximately sinusoidal and are 
given hy y = u sin x, where % is a slowly varying function of 
X. By differentiation we have 

dy . du 

— = u cos X -}- sin X . j- 
dx dx 

* Van due Pol: Proc. I.R.E. (1934), Vol. 22, p. 1067. 
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du 


Due to assumptions respecting a, we will assume that is 
negligibly small. * 


A • 

Ag»m gjs 


-23- COS X 

ax 


u sin X 


Also 


2 ^ ^ 1 . 

^ dx 3 dx 


= — COS X 

4 

So that equation (580) may be written 

diJL 

2u-^ au^(l — u^/4:) = 0 


(681) 


Using the transformation wv?‘ = 1, we have 
dwjdx + ci{w — i) = 0 
The solution of this equation is given as 

M, == J[1 + £-<'«-«.)] . (682) 





Fig. 67 


where is a constant of integration. Then 

tt = 2[1 + 

y = u sin x 

^ l(v—G\ sin cga^ 

■ ■ ^ v( 3/Sy ) ^ [1 + 

After a very long time the amplitude of e will increase to 
2v'[(y-^)/3/5y] = 2^(173/3) 

If we cannot make the saving assumption that a is small, we 
must proceed by graphical means. To demonstrate, we will 
consider the case of a series dynamo supplying power to a 
separately excited shunt motor (Fig. 67). The armature of 
13 *— (T.152) 
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the latter merely oscillates. The voltage equation becomes 
V = JjBq L^pi -}~ R^t -f" Lj^pi e 

and for increasing speed we have ^ = K<l>i and e = a^o), 
where a and K are constants. 

Let Rq-{- Bi = B and Li = L. Then we have 

clK C 

V — iB Lpi + chfxo = iB + Lpi + -j- <t>^ I idt (684) 

For a series dynamo the connection between the voltage 
and current is 

V — pi — Jyi® 

where p and y are constants. 

The equation becomes 

ClK C 

Bi + Lpi + -j- <f)^ j idt — {p— iyi^)i = 0 
clK 

or Lp^q + Bpq + -j- <f>'^ q - [p — )’“]?? = 0 . (686) 

Let t — xV LG. In other words, we are measuring time 
in (periods/27r). Then 

d ( dq dx\ d^q / dxV 1 d^q 


dt\dx dtj~dx^ 

___ . . aK 1 
Writmg -J- <f>^ = we get 


/ dx^ 
\dt) 


LCdx^ 


i I P ^ 4. 1 

O' dx^^ V{LC) ' dx~^^'^ 


1 Y fdq 


LG \dx 


( 686 ) 


With q — y^/[{R — p)LGIy'], we have 

y—a{y— \y^) + ?/ = 0 
where a = ■\/{G{p— R)IL'] 

dZ 

If we write Z — py, pZ = Z we get 


Z^-a{Z-W) + y = 0 


( 687 ) 
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The graphical solution is obtained in the following way. 
First, we construct the curve y — a\Z — as in Fig. 68. 
At any point M, MPQ is drawn parallel with the y axis. At 
the point P, where MPQ cuts the curve, PN is drawn perpen- 
dicular to Ox, and NM is drawn. At the selected point M we 


y 


Fra. 68 

draw the ordinate Mm. Now the subnormal through M, or 
— is Nm, so that equation (687) will now read 

Om = ON -f Nm 

A line RR^ at right angles to MN will be a tangent to the 
required curve. Repetition of this process will give the locus 
of If as a family of spirals. Eventually, the volutes curl asym- 
potically against a certain curve. This curve is a cyclogram, 
and by usual methods can be converted to a wave. Using the 
conversion factors, we revert to an {i, t) wave. This method is 
of the “hit-or-miss” variety: a judicious selection of the 
point M will assist in shortening the work. Hak* has given 
guides to the selection of initial values, based on the following 
development. Equation (686) may be written as 

LpH + {\IC)i + {R- p + iyi^)pi = 0 

On using dimensionless time, we have 

y(z)- “ “ 

* Hak: R.O.E. (1936), Vol. 26, p. 896. 
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With i — y\/[5LG{R — p)/y], we get 

p%y _ (;(1 _ y2)py y = 0 . . (688) 

where a = (p— R)^/{C|L). 

At t = 0, y = 0, and then dyfdt = G^: at any other time 
{y„, tn) the tangent is given by 

(^), ^ “ Jo ^ 

As the equation is vaUd for all values of time, a curve ful- 
filling this equation may be plotted by the usual step-by-step 
methods and the integral evaluated. Hak gives the following 
approximations for the initial values — 

Slope; (7i = 2 + 0 - 16 cr 2 (for u < 1) . . . (590) 

= 1-026 + 0-654a + 0-489/(T (for u > 1 ) (691 ) 

Period: T — tt (for(r<l) . . . (692) 

= 1-989 + 0-745(7 -1- 0-686/(7 (for (7 > 1 ) (693) 

Max. value at t^ax '• 

5-728 

” 3-365 -f- (7 ~ ■ 


( 694 ) 
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EQUIVALENCIES OP OPERATORS 


Operator for 

Developed for y 

p/(p + P) 


pKp + /3r 

[<«-l/(n- l)l]e-^« 

mp + p) 

e-“< — 1 

With coq^ > = 

— a^; tan <l> ~ co/ol 

P*/(P^ + 2a2) + ft)o*) 

— (coq/co) sin (co/ — <f>) 

pKp^ + 2(xp + ft>o*) 

£-(xt (sin o)t)lcjD 

l/(p<“ + 2ctp + Wo**) 

(l/wo®) [1 - (Wo/w) £-“‘sin (wl + ^)] 

With > (jt)^ (roots — m & — n) 

pV(P* + 2(xp + coo*) 

[^e~ne _ rne~^^^y{n — m) 

p/(p* + 2ap + coo*) 

[£-m< — — m) 

l/(p* + 2ap + o>o*) 

(l/wo®) - [(e-»«Vw) 

— (e~’^/n)]l(n — m) 

With a* = Wo’ 

^ (roots equal) 

p 7(P* + 2ap + Wo®) 

e-»« (1 - a<) 

p/(p* + 2ap + Wo®) 

<e-“‘ 

l/(p® + 2ap + Wo®) 

[1 - £-»« (1 + a«)]/wo* 

p7(p + P) 

)i(0) - 

pKp + P)^ 

te-pt 

pKp - P)^ 

— 

pHP + a) (P + P) 

{e~P^ — e~a^)/(a — 

p®/{p + a) (p + P) 

(ae”«^ — jSe~^0/(oc — /3) 

M/(p - a) (p + /3)] - 1 

(pe-^ — (xe~P^)l{(x — j3) 

P 

(y ■— j5)e~oc<+ (a— y)e“"^^+ (j3— a)e“*v« 

(p + a) (p + P) (p + y) 

(a-/S) (iS-y) (y-a) 

p* 

(P + a) (P + /3) (P + y) 

a(/? — y)B-<^ + p(y — a)£“^< 

+ y(a - ^)e-y< 

{« - P)(P— y) (y - a) 

ap/(p* + a») 

sin 
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Operator for 

Developed for y 

+ «*) 

cos a< 

a*p/p + a^) 

2 sin* (^a<) 

pKp + + ia) (P + ^ - ia) 

(sin od . 

pKp — /3 + i«) (P — /S — j(x) 

— (sin cd . bP^)Iol 

p 7(P + + ja) (p + /3 — ja) 

[cos a< — (pjoi) sin a<] 

p7(p - P + }«){p - p~ y«) 

— [cos cd + (/5/a) sin a<] 

pV(P^ — co^) 

cosh (ot 

coplip^ — CO*) 

sinh cot 


(cosh \/t .)/\/(^0 

[p/V(p + /^)] 

[e~^^ \/(^0] cosh 2 a/(/c)<) 

[p/V(P + )?)*] ep/(i>+^) 

[e~^^i/(tr/)] sinh 2 V(pO 

H-* 

II 

1 

i”/nl 

pfi— a+1 = p~(a— n— 1) 

i<*-n-i/T(a - n) 

p— a+1 

tf^-yVia) 

p^ 

i-*ITii) = 1/V(’r<) 

pW) 

1 d! j-i /(A)d!A 

Vtt dJo - A) 

p-\ 

<Vr(i!) = 2 v(<^) 

pi 


pi 

fir-K-e 

pn—l+l 

1,3.5. . .(2n- 1) 

(H* '-W-"-* 

p— n—i+1 

1 . 3 . 5. ..{2n- 1) 

p(p + /?)* 

- 

pip + iS)~* 

e-PflVi^t) 

pip — /3)-* 

ePijVi'^t) (if t < 0) 

P{P + /3)-5 

2e-« ^/itfn) (if <> 0) 

P(P + j8)-“ 

g--j27ra(t)4-M;)e-^« i«~l/r(a) 

p{p— /S)-« 

— . e'"i27ra(t?+f4>) £/9t^-l^r(a) 

P[(P 4- - (P + y)^~“] 

^-2 (e-Pf E-Yt)ir(0L -- 1) if <> 0 
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Operator for H(p) 


Developed for y 




(—^) e-«» 

\g + a/ 


P* 

P + y 


j>3 

p + y 

p 

(p + y) V(P + 

P ViP + P) 

p + y 

pi 

1 + V(M 

VjS _ 1 
V(p + P) 

VHpIP) - 1 ] - 1 

p 

1 + Vipp) 


With p — h^q* 


J/iWl I 

Vv,r/Lrt 2 31 ^‘^ 226 I<S‘’'‘ • 'J 

[\lh v'(»rf)] e-“V4A*i 

^ (2AV<) “ 


V(— y) 

_JL_r 

y 'i/iirt) L 


e~y^ erf [ VC— yO] < > 0 


1 1.3 

1 + irn + 


2yt ^ (2y0* 

, 1 - 3.6 1 

(2y<)» +-**J 


-TT-Jv + Vi— y) e-y* erf [ -/(— yi)] 

Vi’rt) ^ 0 


_ ^ f-L ^ • 

“ Vi^t) L2y< (2y 


ViP-y) 

1 


3 

(2y<T» 

, 1 - 3-6 . 

+ (2y0» • 

e-yt erf ^/i(P - y)<] 


V(ir<) 


+ v'o? — y) fi-y* erf 


pvi^t)'^ PVP 


V[(P-y)n 
^IPeif ViUP) 


Vi'nPi) 


erfc. V(P^) 

erfc. ViP^) 


1 2 t— p 

- g A( 0 ) + ^g^ ^(„pt) 


^ «<//» erfc Vit/P) 
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TRANSIENTS IN ELECTRIC CIRCUITS 


Operator for 

I Developed for y 

i 

i 

pi 

'vp "" V(w0 t ^ 

1 + vm 

.3(|)‘-3.(f, )%...] 

pIV(P* + a*) 

«^0 (“<) 

P/V(j>®— a*) 

Jo 0«<) 

plViP^-p’^) 

(1/^) iiTo (/9<) 

V 

InW 

ip+ Vip^- 1)]" V(p*- i) 

P 

Jf^it) 

[p + V(P* + 1)]" Vip* + 1) 

pI VHP + «■) {p + P)] 

e-i(«+W /<, [J (a - P)l] 

VHP + a)/p] 

e-i“* [a</i (ia<)+ (1 + a^-^o (i“0] 

VCp + «) j 

Vp + V(p + a) 

- i [/i (i«<) + h (i«0] 

p(p® + /3®)““ 

V"-/ < V-i , 

T{oc)\2p) 

Hp + /?)“ + r*]-“ 

Vw/ < \“-i , 

^)(rr) 

pe-^lP 

- [Ji (2 Vt)lvt] 

e-ilP 

Jo (2 V<) 

(1/p) e-^lP 

Ve jj (2 Vt) 

pe-«V(J>’+a*) 


pe-aV(P'-P') 

a]8i:, [p V«* + <*)] 

V(0l} + <«j 

Vioc* — p*) 

V(a» + <»)] 

TT 



APPENDIX II 

MATHEMATICAL FORMULAE 

(I + x)^ I + mx + ^ «* + ... 

(1 + == 1 — a; -f ic* — , . . 

(1 — aj)“^ ^ I + X x^ + ... 

(1 + «)* = 1 + ix - - 

qA 

S±* = 1 i * + - ± _ + . . . 

c^® = cos X + j sin X 
6® = cosh X + sinh x 
sin (« ib !/) = ^ cos V db cos x sin y 

cos (a; ± 1/) = cos x cos y ^ sin x sin y 

sinh (a; ± 2 /) = sinh x cosh y ± cosh x sinh y 

cosh (aj ± y) == cosh x cosh y ± sinh x sinh y 

sinh (a? ± jv) = sinh x cos y cosh a^ sill 
=a V (sinh* X + sin* y) . 
cosh {x ± jy) = cosh x cos y ± i sinh x sin y 
= V (cosh* X — sin* y ) . /<^2 
cos* X + sin* a; = 1 
cosh* X — sinh* a; = 1 

5 = 1 ? — I ? , 

sinh X 1 + (Tr/a?)* ^ 1 + (27r/aj)* 

* coth * = 1 + 1 -+ + 1 + (2./i)i + 

r(l n) — n r(n) = n(n 1) T(n — 1) = n 1 
r(n) . r(l — n) = Tr/sinnTT. 
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Amperes Law, 124 
Approximations, successive, 170 
Armature conductors, eddy current, 
133 

Artificial line, 94 
Attenuation factor, 87 

Ballistic Galvanometer, 62 
Bessel cable, 164 

functions, 154 

operators, 159 

Borel’s Theorem, 11, 56 
Boundary conditions, 47 
Breaking circuit, 48 

CAQuii’s method, 4 
Carson’s equation, 150 
Compensation theorem, 2 
Conductor, eddy current, 125, 160, 
163 

Critical resistance, 23 

Damped line wave, 60 
Damping, 19 

of instruments, 30 

Decrement, 22 
Difference equation, 93 
Duhamel integral, 56 

Empirical formul®, 190 

Faraday’s law, 124 
Ferro -resonant circuit, 187 
Field discharge resistance, 49 
Filter circuit, 96, 166 
Fourier’s integral, 143 

, application, 146 

series application, 141 

Graphical methods, 185 
Guillemin, 12 

Heaviside’s Expansion Theorem, 7 

Notation, 2 

Helmholtz Method, 173 

Impulse generator, 36 

Integral operational evaluation, 63 

Line ; Solution, Characteristic Vibra- 
tions, 90 


Line Solution, Operational, 105 

Steady State, 83 

wave, 109 

Maintained oscillations, 191 
Mathieu, 178 

Motional capacity machines, 65 
Mutually coupled circuits, 39, 42 

Natural period, 23 

Operators equivalencies, 151 
Oscillographs, 30 
increased range, 32 

Plates, eddy current, 127 

, losses, 130 

, penetration, 129 

Reciprocity theorem, 2 
Reflection, 112 
Restriking volts, 51 
Roots of Z(p), 10 
, discussion of, 12 

Shifting theorem, 55 
Starting conditions D.C. motors, 67 
Superposition theorem, 2 
Synchronous machines, 70 

motors, 71 

, pulling into step, 75 

Generators in Parallel, 72 

Taylor’s series, 173 
Test transformers, 42 
Th^venins theorem, 2 
Time constant, 17 
Traffic signals, 49 

Unit Function, 6, 145 

Voltage distribution, 100 

in insulator strings, 98 

in transformer windings, 
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Volterra equations, 180 

Wave form determination, 139 
length factor, 87 
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